










































Aath. Mech. 
2 April 1943 
ie, 


Hilbert 
a und der 
en tiefsten 
s gegeben 
durch die 
war groß 


n Königs- 
abilitierte 
durch den 
tblickende 
zu Seinem 
stillen für 
; schenkte 
er Mathe- 
> er Lehr- 


1eorie der 
ürdig zur 
ingschen 


; über die 
lagen der 
wohl zum 
aupt ganz 
logischen 
jedeutung. 
n nomina- 
g welche, 
rt hat es 
Kantzitat: 
ind endigt 
bert, der 
; Konven- 


Analysis, 
er unver- 
o schönen 
ie ebenso 
Werkzeug 
;h ist, wie 
hätte. Er 


ıtet, nicht 
den. Die 
angspunkt 


‚fassenden 
1 anderer 
\rithmetik 
oxien der 
ue hinein 
n? Diese 

dadurch 
ı Symbole 
nzahl von 
direkt zu 
zu zeigen, 
ng dieses 
»r Schüler 
en, das er 
mel. 477 





termany.— 
am Harz. 


ur a u ar rn ea 7 
propaganda purposes or in the popular press. Quotation of titles and abstracts ) 
in publications of purely technical or scientific nature is not prohibited. 


ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 








Band 23 Juni 1943 Heft 3 
Inhalt: 
Seite ® Seite 
Hau pta ufsätze. F. Schubert: Zur Theorie des Kleine Mitteilungen. H.Görtler: Ergänzende 
stationären Verdichtungsstoßes . - -.. - -» - » 19 Bemerkung zu der Aı beit von D. Küchemann: 
D. T. Dumitreseu: Strömung an einer Luftblase „Störungsbewegungen in einer Gasströmung mit 
im senkrechten Rohr . . 9 a te re 


K.KlotterundG. Kotowski: Über die Stabilität 





der Lösungen Hillscher Differentialgleichungen Buchbesprechungen. E. Eckert: Die Berech 
mit drei unabhängigen Parametern. 1. Mitteilung 149 ei, 2 des w ärmeübergangs in der laminaren 
F. Reutter: Der starre Kreiszylinder im isotropen Grenzsebicht UmBRTOmIET Körper. M Die Tätig- 
istiechen Melle. . 0 0. ah keit der Baltischen Geodätischen Kommission in 
E.-A. Deuker: Zur Stabilität der elastischen den Jahren 1988 bis 1. . . 0... 18 
re Er Zuschriften an den Herausgeber .. . 18 





HAUPTAUFSÄTZE 


Zur Theorie des stationären Verdichtungsstoßes. 
Von Feodor Schubert. 


Aus den dimensionslos gemachten Grundgleichungen werden einige für die 
Berechnung von Verdichtungsstößen wichtige Verhältnisgrößen hergeleitet. 
Die sich dabei ergebende einfache Abhängigkeit von den Stoßbedingungen 
erlaubt es, den allgemeinen stationären Verdichtungsstoß im wesentlichen 
wie einen ihm äquivalenten geraden zu behandeln und die Änderungen 
der 'Zustandsgrößen in leicht zu erzeugenden Diagrammen darzustellen. 


1. Einleitung. 


In Gasströmungen mit Überschallgeschwindigkeit können Verdichtungsstöße auftreten, 
die von B. Riemann 1860 theoretisch vorausgesagt wurden. Beim Durchgang der Strömung 
durch die Front eines solchen Verdichtungsstoßes erleiden Geschwindigkeit, Druck, Dichte, 
Temperatur und Entropie einen Sprung, während die Energie konstant bleibt. Es ist das 
Verdienst von Hugoniot, 1887 diese Verhältnisse im. einfachsten Falle des stationären 
geraden Verdichtungsstoßes einwandfrei aufgeklärt zu haben. Beim gerade Verdichtungsstoß 
erfährt die Geschwindigkeit keinen Richtungswechsel, während die Stoßfront auf der Ge- 
schwindigkeit senkrecht steht. Für den allgemeineren Fall des schiefen Verdichtungsstoßes 
hat Th. Meyer‘) 1907 die Theorie entwickelt. Eine sehr anschauliche graphische Darstellung 
hat A. Busemann?) 1929 durch die Einführung der Stoßpolaren erzielt. Eine solche Stoß- 
polare beschreibt jeweils für eine bestimmte Machsche Zahl der Anströmung vollständig 
die Kinematik der Strömung beim Durchgang durch die Stoßfront. Für die dynamischen 
Verhältnisse, z. B. den Druck- und Entropiesprung, sind zusätzliche Diagramme notwendig, 
die auch von Busemann angegeben wurden. Wenn es nun gilt, die Stoßfront bei der 
Umströmung eines Körpers aufzubauen, was rechnerisch auf eine Kombination der Ver- 
dichtungsstoßgleichungen mit den charakteristischen Differentialgleichungen der Überschall- 
strömungen hinausläuft, so erweisen sich diese graphischen Hilfsmittel als nicht genau genug 
und in der Handhabung mühsam. Deshalb sieht man sich genötigt, eine für die Rechenpraxis 
zweckmäßige und möglichst einfache Darstellung der Verdichtungsstoßtheorie zu entwickeln. 
Eine solche wird im folgenden hergeleitet, wobei die Zusammenhänge aller physikalisch 
bedeutsamen Größen an der Stoßfront durch Vermittlung eines Parameters wiedergegeben 
werden. Es zeigt sich dabei, daß der schiefe - Verdichtungsstoß in weitem Maße auf den 
geraden Verdichtungsstoß zurückgeführt werden kann. 





1) Th. Meyer: Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgeschwindigkeit 
strömt. Forsch.-Arb. Ing.-Wes., Heft 62, Berlin 1908. 

2) A. Busemann: Artikel „Gasdynamik‘“ im Handbuch uer Experimentalphysik, Bd. 4, 1. Teil, Leipzig 1931. 
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2. Bezeichnungen. 


Mit p werde der Druck bezeichnet, mit o die Dichte, mit T die Temperatur, mit ®w die 
Geschwindigkeit, mit a die Schallgeschwindigkeit; k sei das Verhältnis der spezifischen 
Wärmen. Die Größen ünmittelbar vor dem Stoß sollen 
den Index 1 tragen, die hinter dem Stoß den Index 2. 
nn, Durch das betrachtete Flächenelement der Stoßfront wird 
diejenige Normalebene gelegt, die den Vektor der Anström- 
geschwindigkeit w, enthält. Die Spur des Stoßflächen- 
elementes in dieser Ebene bildet mit der Anströmge.- 
schwindigkeit w, den „Stoßfrontwinkel“ y (Bild 1). 
7\ Ks] Der Stoßfrontwinkel y ist der Parameter 
in der neuen Darstellung der Verdichtungs- 


73 








Bild 1. Änderung der Geschwindigkeit w i Pat x . 7 

im Verdiehtungsstoß. stoßtheorie. Sein Bereich erstreckt sich von — (ge. 
rader Verdichtungsstoß) bis zum Machschen Winkel «, der Anströmung (Verdichtungsstoß 
mit verschwindender Stärke): 


n 
sys 77 . 


Mit w,„ und w; werden die Geschwindigkeitskomponenten normal und tangential zum Stoß. 
flächenelement bezeichnet; a* sei die kritische Schallgeschwindigkeit. 


3. Verhältnis der Drucke, Dichten, Temperaturen und Normalgeschwindigkeiten 
vor und hinter dem Stoß. 


Bei Vernachlässigung von Wärmeleitung und innerer Reibung liefern die Erhaltung der 
Masse, der Energie und des Normalimpulses die folgenden grundlegenden Gleichungen für 
den Verdichtungsstoß : 

0, Wnı = 0, Win 
: na, ,. kp» 
TE Fi A a hr > 
am +tp = Win t+P; 


Der tangential zum Stoß angesetzte Impulssatz sagt zusam- 
men mit der Kontinuitätsgleichung aus, daß die Tangential- 
geschwindigkeit w; vor und hinter dem Stoß dieselbe ist 
(Bild 2). In dem Energiesatz (2) werden von beiden Steiten 


2 
. w 2 } 
der Gleichung 5‘ subtrahiert, so daß in dem resultierenden 





Gleichungssystem nur noch die Normalgeschwindigkeiten 

Bild : a auftreten. Ferner werden die Gleichungen zweckmäßig 
ild 2. Zerlegung der Geschwindig- 2 r s R 

keiten w, und w2. mittels o, und a, dimensionslos gemacht, indem man die 

erste Gleichung durch o,a,, die zweite durch a}, die dritte 

Pı 
mn 


durch o,a} dividiertt. Da a? ist, wird bemerkenswerterweise der eine auftretende 


. . q 
dimensionslose Parameter _Pı _ —. Man erhält 
30a & 
Yn _ Bı 9m 
0, 4, 


k „| Mg 
7 Rt 0, u 0, ey u Mm . . . . . . . . . (2a), 


2a! 
—. 2 n LP 5: B 
0,0 os re ee ee (3a). 
Dies sind drei Gleichungen für die Dimensionslosen », n, er} Als einziger Para 
1 1 1 1 


enthalten. Da w„,=m,siny ist, kann man 


(la), 


meter ist in dem Gleichungssystem - ku 
a, 


m — Ma,siny setzen, wobei Ma,="" die Machsche Zahl der Anströmung ist. Indem 
1 1 


"ns nach (1a) in (2a) und (3a) eingesetzt wird, erhält man 
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1 FRLT Bu; 
k—1l 2u k- 
wi 1 


‚2 
= al ige De BEN RER 1 LT RE ERF BHrn 06 ur 3 2 7 5 Fe 


(2b), 


Mininnti Pı ; 
Durch Elimination von folgt: 


0; 
vd 1 Z=0 
k—i 20 k—-lo, ki-io, kik-la 5 20 


za 


A.” lı o&5 __k wu k mn,0; Im, 


wofür man auch schreiben kann: 


( 1 0 k DR (0; +3 (2 1)=0 
k—1lo, k-lar/\e, su 


Nach Absonderung der trivialen Lösung o,=o, erhält man 


| 1 n 2) 0, k+1l mn, 
k—1 


"as 
2a} 


. SG-DeE'’ 


C 


k+1lm k+1 
1 ? sın? „ 
er 5 Mai sin?’ y 
2 0 # 2 
k—I mn, 


‚81 
1+ Pr} 14755 Ma}! sin? y 


Aus (3b) und (6) ergibt sich 


Ps _ „Mai sin’y — 1,1 
Tu k+1 k 


1 
Dividiert man durch Pı ;=7, so erhält man 
o,a k 
k—1,. 
Pr ZT 
k+1 


p, k+1 


Ma? sin? 


Weiter ist nach (1) und (6) 


k—1 
1-+ 5 Mai sin? y 
£ Ei 


o # 


®? sin? „ 
Ma; sin? } 


: RE i 
Da nach der Gasgleichung 1-2. ist, so folgt 
1 ıC2 


k—1\/(k—1 i 
(% Ma} sin? y — 5 -) " 5 Ma; sin’ y+ 1) 


- — u ea. Me 
5) Ma} sin? y 


In diesen Beziehungen (6), (8), (9), (10) ist die Theorie des schiefen Verdichtungsstoßes auf 
die des geraden zurückgeführt, wobei an Stelle der Gesamtgeschwindigkeiten beim geraden 
Verdichtungsstoß hier die Normalkomponenten der Geschwindigkeiten auftreten. Die Mach- 
sche Zahl der Anströmung ist durch die Machsche Zahl der Normalanströmung Ma, sin y zu 
ersetzen. Man muß sich aber vor dem Trugschluß hüten, daß man jede beliebige Beziehung 
aus der Theorie des geraden Verdichtungsstoßes formal auf den schiefen Stoß übertragen 
könnte, indem man an Stelle der vollen Geschwindigkeiten einfach die Normalgeschwindigkeiten 
tretendäßt. Beispielsweise gilt für den geraden Verdichtungsstoß die Prandtlsche Gleichung 


we, 


während für den schiefen Verdichtungsstoß 
3) Wie nachträglich bemerkt wurde, hat L. Crocce in L’Aeroteenica Bd. XVII (1937), 519—536, (Singolaritä 
della eorrente gassosa iperacustica nell’ intorao di una prora a diedro, GI. (10’)) Formeln angegeben, die den 
Gleichungen (6) und (8) entsprechen. 
9*® 
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wu, N, mat? — w; 
ist, wie schon Th. Meyer a. a. O. gezeigt hat. 

Eine Beziehung für die vollen Geschwindigkeiten beim schiefen Verdichtungsstoß, die 
ebenfalls nur den Parameter Ma, sin y enthält, ergibt sich aus (2) unter Verwendung von (6) 
und (7): 

wi; — m; 4 (k Ma; sin’ y+1) (Ma; sin’ y — 1) in 

" 5) Ma? sin? y a 

5) Mai sin? ) 

Für die praktische Berechnung von w, dürfte allerdings diese Formel weniger in Frage 
kommen. 


4. Drosselfaktor und Entropiesprung. 


Um Druck und Dichte in der Strömung vor dem Stoß bei gegebenem Geschwindigkeits- 
feld auszurechnen, hat man die Energiegleichung 


up, ww. 1841 , k » 
ZU u 3 


uf de (12a) 
0 
und die Poissonsche Gleichung 
p 


es 133 
ok ( 1) 


zu verknüpfen. Dabei bezeichnet p, den Kessel- oder Ruhedruck, o, die Ruhedichte. Daraus 
ergibt sich 
k—1 


)\ k o \k-1 k—1/w\? 
(2) (2) —=1- =) N Dr a 18 
Pu 0, k+ 1 \“ 
Bezeichnet man mit p) den Ruhedruck hinter dem Verdichtungsstoß, mit o, die betreffende 
Ruhedichte, so hat man entsprechend in der Strömung hinter dem Stoß als Energiegleichung 


2. # 18:+1 k » 


FD *2 
k—1lo 2 


(12b) 


„er "EIz 
und als Poissonsche Gleichung 
(13b). 


Druck und Dichte ergeben sich aus 
k—1 


pi) * _fe? ee 
(4) =(2) =1 k+1 =) er a 2 Pre 
Po 


AFTER » Me oe act 
Die beiden Verhältnisse = und — sind nach dem zweiten Teil der Gleichungen (12a) und 


{) 00 
(12b) einander gleich und sollen als „Drosselfaktor“ K bezeichnet werden. Wendet man 
die Poissonsche Gleichung in der Form (13a) und (13b) auf die Drucke und Dichten un- 
mittelbar an der Stoßfront an und dividiert diese Beziehungen durcheinander, so erhält man 


P: (©) 2 Br ) 
v:\o) »i\0) ’ 


was nach der Definition des Drosselfaktors gleich K*-! wird. Folglich ist 


k 1 
K=(&)' (2)' . 
\e Ps 


Durch Einsetzen von (6) und (8) bekommt man schließlich 


(16). 


5) 


k+l k 
K=-)" a Ma} sin? y ' k—1 1 


k—1 \ ’ 
1-++ 5 Ma} sin? y k Ma; sin’ y— ——5 


_ 


Durch den Drosselfaktor ist die Berechnung von Druck und Dichte in der Strömung hinter 
dem Verdichtungestoß aus dem Geschwindigkeitsfeld ermöglicht. 
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Der Entropiesprung an der Stoßfront läßt sich leicht mittels des Drosselfaktors ausdrücken. 
In mechanischen Einheiten ist die Entropie der Masseneinheit mit R als Gaskonstante 


IB, die 
‘on (6) also ist der Entropiesprung 


(11). 
Mithin ist nach (15) 
RlnK. 


Frage 


5. Der Ablenkungswinkel. 


Zur Berechnung des Ablenkungswinkels ö der Strömung durch den Verdichtungsstoß 

keit (siche Bild 2) geht man am besten von den Beziehungen 
rkeits- 
Wu, ir + ti 


MW, 
(12a) 
Dureh Division ergibt 


(13a) 
woraus nach (9) 
araus Ba 
Ma? sin? y 
0) u ’ 
ALRETAE 
(14a) >— Ma; sin’ y 
Tende folet. Diese Gleichung löst man nach Ö auf, indem ein trigonometrisches Additionstheorem 
‘hung auf tg (y - Ö) angewandt wird: 


k+l 
120) 5 Ma} 
cote d Ma! sın’y 1 l/tgy 


In dieser Formel treten Ma, und sin y getrennt als Parameter auf. 


13b). 


6. Hinweise für den praktischen Gebrauch. 
Abschließend werden die gewonnenen Formeln für die beim Verdiehtungsstoß maßgeb- 
14b). lichen Größen zusammengestellt: 


3 J 
und = Ma? sin? y 
man 
| un- 
man 


Ma? sin’ y+1 


a 
+1’ a sm®y— 77 


k—1\ /k—1 
k Ma? sin? y = | 3— Mai sin? y 


e + 


.) 


\ ara: sin? y 


Ma; sin? y 1 


(16). "2 5 Ma} sin? y k Ma? sin? y — 


nter 
cotg h) 
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Zahlentafel i. 





PR 0% T; 


Die Dimensionslosen j 
u MM 1 


und K in Abhängigkeit von Ma, siny für k= 1,405. 


Pa T, 7 Pa 


Ma, sin y 


Ma, sin | 
1 7 Pı | 7, Pı 





1,00 | 1,000 | 1,000 | 1,000 10,35 | 3,834 | 2,69 | 0,831 
1,05 1,084 | 1,033 | 1,000 10,70 | 3879 | 2759 | 0317 
1,10 | 1,169 | 1,066 | 0,999 11,06 3,923 2,820 0,303 
1,15 | 1,254 1.098 0,997 3% 11,43 3.965 2,881 0,291 
1,20 | 1,341 1,129 0,993 3, | 11,80 4,006 2.944 0.278 
1,25 | 1427 | 1,161 0,987 12,17 | 4046 | 3008 | 0867 
1,30 | 1514 | 1198. |) 0,979 12,56 | 4,086 | 3,073 | 0,255 
1,35 | 1601 | 125 | 0,9m 1294 | 4124° | 3139 | 0,2447 
1,40 | 1687 | 1,257 | 0,958 1334 | 4,161 | 3206 | 02344 
1,45 | 173 | 1290 | 0985 1374 | 4197 | 3273 | 02246 
1,50 | 1859 | 1394 |) 0980 | 1414 | 4232 | 3392 | 02151 
1,55 | 193 | 1358 | 0913 3,5: 14,56 4266 | 3,412 | 0,2061 
1,60 | 2,027 1,392 | 0,895 | 1497 | 4,300 | 3,482 0,1975 
1,65 | 2,110 | 1,428 0,876 15,40 | 4332 | 3554 | 0,1802 
1,70 | 28 1,463 0,856 15,83 4,364 | 3,627 0,1813 
1,75 | 2,273 1500 | 0,835 16,26 | 4,395 3,700 0,1737 
1,80 | 2353 | 1538 | 0813 16,70 | 4485 | 3775 | 0.166 
1,85 2431 | 156 | 01 || 17,15 | 4454 | 3850 | 0,1596 
1,90 2,508 | 1615 | 068 |ı | 1760 | 4,483 | 3,927 | 0,1530 
1,95 2,583 1,655 | 0,745 | 18,06 4,511 4,004 0,1467 
2,00 2657 | 1665 | 0722 |) 1853 | 4538 | 4,083 | 0,1407 
2,05 | 2730 | 1237 | 069 | 19,00 | 4,564 | 4,162 |. 0,1350 
2,10 2,801 1,79 | 065 || | 1947 | 45900 | 4242 | 0129 
2,15 2811 182 | 062 | 1995 | 4615 | 434 | 0,122 
2,20 2,9399 | 1867 | 0,69 20,44 | 4639 | 4,406 | 0,1192 
2,25 | 3,006 | 1912 | 0607 || 2094 | 4.663 | 4,489 | 0,1144 
2,30 301 108 0 21,44 | 4,687 | 4,574 | 0,1098 
2,35 3,135 | 2005 | 0563 || 2198 | 4700 | 4,659 | 0,1054 
2,40 3,197 | 2,052 | 0,54 2245 | 4,731 | 4,745 | 0,1012 
2,45 32358 | 2101 | 03 || 297 | 1758 | 4833 | 0.0972 
2,50 3317 | 23151 | 0501 | | 2349 | 4,774 | 4921 | 0,0934 
2,55 3375 | 2201 0,481 24,02 | 4795 | 5010 | 0,0897 
2,60 3,132 | 233 | 0402 || | 2455 | 4815 | 5,100 | 0,0862 
2,65 3a87 | 2305 | 048 |! 3,10 | 4,834 | 5,191 | 0,0828 
2,70 3500 | 2358 | 06 || | 35,64 | 4853 | 5,283 | 0,0796 
2,75 3392| 2113 | 0408 || 26,19 | 4,872 | 5,376 | 0,0766 
2,80 3,643 | 2468 | 0,892 |) 26,75 | 4890 | 5470 | 0,0736 
2,85 3,098 | 254 | 085 || 27,32 | 4,908 | 5,566 | 0,0708 
2,90 3741 | 2581 |) 0360 | | aus8 | 4925 | 5662 | 00681 
2,95 3788 | 269 | 035 | 2846 | 492 | 5,759 | 0,0656 
29,04 | 4959 | 5857 | 0,0881 
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Zahlentafel 2 (1. Teil). Ablenkungswinkel ö in Abhängigkeit vom Stoßfrontwinkel y 


bei verschiedenen Machschen Zahlen Ma, für k == 1,405. 














Ma, = a) 
1.10 1,20 1.30 1,40 1.50 1.60 1.70 1.80 1.90 2.00 













0) 0 V 0 0 0 














































73 2.05 2,37 2,69 3,00 3.30 
86 0,65 1.34 2,02 2,71 3.38 4,04 4,69 5,31 5,92 6,50 
8 0,94 1,95 2,97 3,98 4,97 5.94 6.89 7.80 8.68 9,52 
82 1,18 2,51 3,84 5.16 6.45 773 8.093 10.09 11.21 12.28 
80 | 1,37 3,00 4,62 6,23 7,79 9,31 10,76 12,16 13,48 14,73 
78 1,48 3,39 5,29 7.16 8,97 10,71 12,37 13,96 15,45 16,86 
76 1,51 3,69 5,84 7,94 9,96 11,90 13,74 15,48 17,11 18,64 
74 1,45 3,87 6,25 8,56 10,77 12.88 14,86 16.73 18.47 20.09 
12 1,29 3,93 6,52 9,01 11,38 13,63 15,74 17,70 19,52" 21,21 
70 1.03 3,87 6,64 929 11,80 14,16 16.37 18.41 20.29 92.03 
68 0.66 3.68 6,61 9,40 12.03 14,49 16,77 18,87 20,80 22,56 
66 0,17 3,36 6,43 9,34 12,07 14,60 16,94 19,09 21.05 22,84 
64 2,91 6,10 9,11 11,93 14,53 16,91 19,09 21,08 22,88 
62 2,32 5,62 8,73 11,61 14,26 16,69 18,90 20,90 22,71 
60 1,60 5,00 8,18 11.13 13,83 16,28 18.51 20.52 22.34 
58 0,75 4,24 7,49 10,48 13,22 15,70 17.95 19.97 21,79 
56 3.33 6,64 9,69 12,46 14,96 17.22 19.26 21,08 
54 2,29 5,66 8,74 11,54 14,07 16,35 18,39 20,22 
52 1,12 4,53 7,66 10,49 13,04 15.33 17.38 19.22 
50 3,27 6,43 9,29 11,87 14.17 16.24 18.08 
48 1.87 5,07 7.96 10,56 12,89 14,96 16,82 
46 0,34 3,57 6,50 y12 11,47 13.57 15.44 
44 1,94 1,90 7,56 9,93 12,05 13.94 
42 0,17 3,17 5,86 8,27 10,41 12,33 
40 | 1,30 4,03 6.48 8.65 10,60 
38 2,07 1,55 6,77 8,75 
36 2,49 4,75 6,77 
34 0,29 2,60 4,66 
32 0,29 2,41 





0,00 







Der Zusammenhang der physikalischen Größen wird 
in dieser Darstellung fast ausschließlich durch den 
zusammengesetzten Parameter Ma, sin y vermittelt. 
Nur in dem Ausdruck für ö treten Ma, und y als 
zwei getrennte Parameter auf. Dieses Gleichungs- 
system hat sich bei allgemeinen Betrachtungen in 
der formelmäßigen Rechnung praktisch bewährt, 
für zahlenmäßige Auswertungen bei speziellen Bei- 
spielen erwies es sich als besonders brauchbar. 
Nach den Formeln (I) bis (V) sind die Zahlentafeln 
1 und 2 berechnet. Die Bilder 3 und 4 bringen 
graphische Darstellungen dieser Formeln. In Bild 4 
sind die Maxima der einzelnen Funktionen durch 
eine gestrichelte Kurve verbunden worden, deren 
Gleichung angebbar ist. 
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2° Bild 3 (links). 
je Die Dimensionslosen n, 82 nr und K in Abhängigkeit von 
1 1 


01 















| 
0 Ma, siny für k= 1,415. 
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Zahlentafel 2 (2. Teil). 





220 2,60 2,80 3.00 3.20 





0 0 0 0 0 0 0 0 0 
3,86 1,37 1.83 5.25 9,62 5.96 6.26 6,54 6,78 7.00 
1.59 8,57 9,47 10,27 10,98 11.63 12,20 12,71 13,17 13.59 
11,08 12,49 13,75 14,87 15,87 16,76 17,55 18.25 18,88 19,44 
14,25 16,00 17,56 18,94 20,15 23.28 22,17 23.01 23.75 24.41 
17.03 19,06 20,85 22,41 23,77 24,97 26.01 26,93 27,7 28.46 
19,42 21.65 23,59 25,28 26,74 28,00 29,11 30,08 30,92 31,66 
21.40 23.77 25.82 27,58 29,09 30,39 31,52 32,51 33.36 34.11 
22,98 25,45 27,56 29,36 30,89 32,21 33.34 34,33 35,18 35.92 
420 26,73 28,87 3).08 32.22 33.53 34,65 35.62 36.46 37,20 
25,08 27,64 29,79 | 31.60 38 13 34,42 35.53 36.48 37,31 38,02 
25.65 28,23 30,37 T 32,16 33,67 34,95 36,04 36,97 37.77 33.47 
25,95 28,52 30,65 32,43 33,92 35,17 36,24 37,15 37.93 38,61 
26,00 28,56 30,68 32,43 33.90 35.13 36,17 37,06 37,83 38,49 
25.83 28,38 30.47 32.20 33.65 34,86 35,88 36,76 37.50 38,15 
25.45 27.99 30,07 31,78 33,20 34,40 35,40 36.26 36,99 37.62 
24,90 27,43 29,49 31.18 32,59 33,76 34,75 35,59 36.32 36,93 
24.19 26,70 28,75 30,43 31,82 32,98 33,96 34,79 35,50 36,11 
23.33 25,84 27.87 29,54 30,92 32.07 33.04 33,86 34,56 35, 
22,33 24,84 26,87 28.53 29,90 31,05 32,01 32,82 33,52 34,12 
21,21 23,72 25,75 27,41 28,78 29,92 30,87 31,69 32,38 32° 
19,96 22.48 24,52 26,18 27,55 28,69 29,65 30,46 31,15 81.7 
18.61 21,14 23.19 24,86 26,24 237.38 28.34 29,15 29,85 
17,14 19,70 21,76 23.45 24,84 25,99 26.96 27,77 28,47 
15.57 18,16 20,25 21,95 23.35 24,52 25,50 26,32 27.03 
13,88 16,51 18,64 20,37 21,79 22,97 23,97 24,80 25,52 
12,09 14 77 16,94 | 18,70 20,15 21.36 22,37 23,22 23,95 
10,19 12,93 15,14 | 16,94 18,43 19,66 20.70 |! 2157 
8,16 10.97 13,25 15,10 16,63 17,89 18,96 19,86 20,62 
6,01 8,90 11.24 13,16 14.73 16,04 17,14 18,07 18,86 
3,71 6,70 913 | 11,11 12.74 14,10 15,25 16,21 17,04 
1,25 4,35 6.88 8,94 10,65 12,07 13,26 14,27 15,13 
1,84 1.48 6,64 8.43 9,92 11,17 12,23 13,14 
1.90 4,18 ı 6.07 7,64 8,97 10.09 11.05 
152 | 853 5,20 6,62 7,81 8,84 
0,77 2,57 4,09 5,38 6.48 
1.33 2,73 3.93 
1,12 


00 0 
X © 8 
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er er er 


0° 














Bild 4 (rec. ‘s). 

Ablenkungsw. kel d in 
Abhängigkeit v.m Stoß- 
frontwinkel y mit aa: als 
Parameter für k= 1,405. 
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Zahlentafel 


” 








7,46 
14,44 
20,60 
25.76 
29,92 
33,17 
35.62 
37,41 
38,65 
39,44 
39,85 
39.95 
39,80 
39,42 
38,87 
38,16 
37.31 
36,35 
35,29 
34,14 
32,91 
31,61 
30,24 
28,82 
27,34 
25,80 

24.21 
22,57 
20,88 
19,13 
17,32 
15,44 
13.49 
11.45 

9,29 

7,00 

4,51 

1.76 


0 
7.82 
15.11 
21,49 
26,79 
31,02 
34,29 
36,75 
38,52% 
39,73 
40.48 
40,86 
10.93 
1,75 
10,35 
39,77 
39.04 
38.18 
37,21 
36,14 
34,98 
33,75 
32,45 
31,09 
29,67 
28.20 
26.68 
25,11 
23,50 
21,84 
20,13 
18.36 
16,54 
14,66 
12,70 
10,64 
8,47 
6,14 
3,59 


0.73 


38,27 
40,00 
41,17 
11,87 
12,21 
42,23 
12.01 
11.57 
10,96 
10,21 
39.33 
38.34 
37.26 
36.09 
34.86 
33.56 
32,20 
30,80 
29,34 
27,84 
26,30 
24,72 
23,10 
21,44 
19,74 
17,99 
16,19 
14,34 
12,41 
10,41 

8,29 

6,01 

350 

0.64 


39.82 
11.53 
42,64 
43,29 
13.57 
43.55 
43,28 
42,81 
42,17 
11.38 
40,48 
39,47 
38.: 


37,2 


0 
9,18 
17,58 
24.72 
30.47 
34,90 
38.20 
40.59 
42,26 
13,34 
43,96 
14,21 
14,17 
43.88 
43,39 
42.73 
41,93 
41,02 
40,00 
38,90 
37,72 
36,48 
35,19 
33,84 
32.44 
31.01 
29,54 
28.04 
26,51 
24,95 
23.36 
21,75 
20,11 
18,44 
16,74 
15,02 
13,25 
11,44 
9,57 
‚61 


5 
D) 
‚4 


0 
9,7: 
18,5: 
25.9 
31,85 
36.32 
39.61 
11.96 
43.58 
44,60 
45,17 
45,37 
45.28 
14.95 
44,43 
43,74 
42,92 
41,98 
40,95 
39,84 
38.65 
37,41 
36,11 
34,76 
33.38 
31,96 
30.50 
29,02 
27,52 
25.99 
24,45 
22,88 
21,31 
19,71 
18,11 
16,49 
14,87 
13,24 
11,60 
9,95 
8,30 
6,64 
4,99 
3,33 
1,66 
0 
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ein und benutzt man die Tatsache, 
ı 


an . A _ 
Führt man noch den Machschen Winkel «,=are sin Ma 
IM 


daß für Fa = const. die in (I) bis (IV) links stehenden Dimensionslosen unverändert bleiben, 
2 Ya | 

so kann man leicht zu einfachen Diagrammen gelangen, welche in der Praxis, besonders bei 

überschlagsmäßigen Rechnungen, gut verwendbar sind. Man braucht nur jeweils in einem 

sin «,, siny-System eine Anzahl ausgewählter Geraden durch den Koordinatenursprung zu 

ziehen und an sie die entsprec henden Zahlenwerte zu schreiben. Bild 5 zeigt als Beispiel 

ein derart angelegtes Diagramm. 


20 
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r - " a Da r - - , 
Bild 5. Geraden gleicher ““ -Werte im sin «4, sin y-Feld. 
[#) 
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Für den rechnerischen Aufbau der Stoßfront ist es sehr wichtig, den Zusammenhang 
zwischen w, und ö zu haben, um die Verknüpfung zwischen Verdichtungsstoßgleichungen 
und den charakteristischen Differentialgleichungen für Überschallströmungen direkt vornehmen 
zu können. Zu diesem Zweck wird ö bei der betreffenden Machschen Zahl Ma, der An- 
strömung als Funktion von y nach Formel (V) in dem in Betracht kommenden Teilbereich 
berechnet, der sich meist leicht abschätzen läßt. Um w, ebenfalls als Funktion von y zu 
berechnen, geht man etwa von der Beziehung 

n;—= m, C0Sy= m, cos (y — Ö) 
aus, die sich aus Bild 2 sofort ergibt. Danach ist 
cos y 
u (y— 6)’ 
womit auch w, über das bereits berechnete ö leicht als Funktion von y ermittelt werden kann. 
Das angedeutete Vorgehen beim Aufbau der Stoßfront erwies sich in der rechnerischen 
Einfachheit und in der Genauigkeit früheren Verfahren als überlegen. 

Zum Schluß möchte ich Herrn Prof. Tollmien für die freundliche Beratung bei der 

Abfassung dieser Arbeit aufrichtig danken. 449 
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Strömung an einer Luftblase im senkrechten Rohr. 
Von D. T. Dumitrescu in Bukarest. 


Die Aufgabe: Form und Geschwindigkeit einer in einem vertikalen Rohr 
aufsteigenden, unendlich langen Luftblase zu finden, führt unter Vernach- 
lässigung der Kapillar- und Zähigkeitskräfte auf ein Eigenwertproblem 
der Differentialgleichung Ap=0, das in erster Annäherung unter der der 
Wirklichkeit gut entsprechenden Annahme gelöst wird, daß der Meridian- 
schnilt der Blase in der Umgebung der Kuppe durch den Krümmungskreis 
erselzt werden kann. Der Eigenwert, der sich als Funktion dieses 
Krümmungsradius angeben läßt, wird dadurch bestimmt, daß die asym- 
ptotische Lösung mit stetiger Tangente in den Krümmungskreis übergehen 
soll. Berechneter und experimentell gefundener Wert stimmen gut überein. 


I. Einleitung. 


Aufgabenstellung und Formulierung der Randwertaufgabe. Die Strömungen schwerer 
Flüssigkeiten mit freier Oberfläche sind trotz ihres häufigen Auftretens verhältnismäßig 
wenig untersucht. Von den bekannten Arbeiten nennen wir diejenigen von Blasius'), 
welcher sich mit einem ebenen Problem beschäftigt hat, und eine von Förster?), welcher 
eine Reihenentwicklung im Raum für rotationssymetrische Fälle angegeben hat. 


Die folgende Arbeit’) handelt von der Form und der Aufsteigegeschwindigkeit einer 
unendlich langen Luftblase, die in einem vertikalen, flüssigkeitsgefüllten Kreisrohr aufsteigt. 
Die Form der Blase ist aus Bild 1 zu ersehen*). Die dreidimensionale, rotations- 
symmetrische Strömung wird als stationär und wirbelfrei betrachtet, und zur 
Lösung des Problems wird die Potentialtheorie herangezogen. 
Der Einfluß der Kapillarkräfte darf, besonders bei größerem Durchmesser 
des Rohres, vernachlässigt werden. 
Um die Randbedingungen auf der freien Oberfläche zu erhalten, betrachten 
wir allgemeiner eine Blase von einer bestimmten Dichte ö,, die in einer anderen 
Flüssigkeit, deren Dichte ö, ist, aufsteigt. Beziehen wir den Vorgang auf ein 
relativ zur Blase ruhendes Koordinatensystem, so handelt es sich um eine 
stationäre Strömung. 
Ist dann Oz die Symmetrieachse des Rohres, p, der Druck, w die Ge- 
en schwindigkeit eines Teilchens und 9=-+gz das Potential der äußeren Kräfte 
mit Wasser ge- in der Flüssigkeit ö, (z nach oben positiv), so lautet dieBernoullische Gleichung 


füllten Rohr auf der Trennungsfläche 
aufsteigende 


2 
Luftblase. mt .2 
Ö ur 7 92 == 0 

; 2 


. . . . . 7 
und wenn wir den Druck im Staupunkt (2=0) mit p, bezeichnen, so ist ec, = z ‚ also 
1 
.2 


p, =Pmp+g95,2— m e; 


Andererseits besteht in der ruhenden Blase die Gleichung 


wo (,— a8 ist, wie sich leicht ergibt, wenn wir verlangen, daß p,=p, für 2=0. Daraus folgt 


P:=P +9 6,2 


und weil auf der Trennungsfläche der Druck beiderseits derselbe sein muß (p,=p,), folgt 


w=/252l1 3) 


1) H. Blasius: Funktionentheoretische Methoden in der Hydrodynamik. Z. Math. u. Phys. Bd. 58 (1910), S. 90. 
2) Förster: Über Flüssigkeitsstrahlen, deren Formen Drehungskörper sind. Z. Math. u. Phys. Bd. 62 (1912), S. 319. 
3) Für die Anregung zu dieser Arbeit und für die wertvollen Ratschläge bei ihrer Durchführung bin ich meinem 
hochverehrten Lehrer Herrn Professor Dr. L. Prandtl zu besonderem Dank verpflichtet. 

4) Photographie, vgl. S. 148, 





. u. , ' r Z. angew. Math. Mech. 
140 Dumitresceu, Strömung an einer Luftblase im senkrechten Rohr Bd.?3 Nr.3 Juni 194 





oder wegen d,<6, folgt als Randbedingung für die Geschwindigkeit auf der Blase einfacher 
w—=yY2gz. Die mathematische Formulierung unseres Problems lautet jetzt, wenn wir mit 
p das Geschwindigkeitspotential bezeichnen, eine Lösung der Differentialgleichung 197 =0 zu 
finden, so daß 1. auf einer noch nicht bekannten Stromlinie w = y29z wird, und 2. die Wände 
des Rohres selbst Stromlinien sind. 

Es handelt sich also um kein eigentliches Randwertproblem, und die mathematischen 
Schwieriekeiten bestehen in der Tatsache, daß die erste Randbedingung auf einer noch nicht 
vorzegebenen Kurve zu erfüllen ist, deren Form gleichzeitig mit der Lösung des Problems zu 
bestimmen ist. Es wird sich zeigen, daß es sich um ein Eigenwertproblem handelt. 

Das im folgenden durchgeführte Verfahren zur Lösung des Problems beruht auf zwei 
Entwicklungen, einmal einer Reihenentwieklung vom Staupunkt (2=0) aus, dann auf einer 
asymptotischen Entwiel lung von z= »> herein. Diese beiden Entwicklungen werden stetig 
und mit stetiger erster Ableitung "aneinander angeschlossen. 


II. Theoretischer Teil. 


1. Dimensionsbetrachtungen. Die Aufsteigegeschwindigkeit der Blase läßt sich bis auf 
eine Konstante, den Eigenwert unseres Problems, durch Dimensionsbetrachtungen bestimmen, 
Zu diesem Zwecke wollen wir das Z/-Theorem von Buekingham’) in der von Bridgman‘) 


angegebenen Form anwenden. 

Sind a, ß, y, ... meßbare Größen und Dimensionskonstanten einer Naturerscheinung, dann 
lautet das //-Theorem folgendermaßen: 

Wenn die Gleichung ® (a, P,y,...)=® in den »n Veränderlichen a, ß, y, ö,... eine „voll- 
ständige Gleichung“ sein soll (d. h. eine Gleichung, die ohne Änderung ihrer Form richtig 
bleibt, wenn die Größen der Grundeinheiten irgendwie geändert werden, denn nur auf solche 
Gleichungen ist natürlich die Dimensionsanalysis anwendbar), dann hat die Lösung die Form 


v0 


FiR,2,,8 rs, 


a1 + ++ 


wo die 7; die na — m unabhängigen Produkte der Argumente a, f, y, d, ... usw. darstellen, 
die dimensionslos in den Grundeinheiten sind und = die Anzahl der Grundeinheiten ist. Man 


kann die Lösung nach einem Produkt auflösen, und man bekommit 
Sy u u A 5 SE Re 5 


wo die p, gr, ... so beschaffen sind, daß a ' pP” y" 0"... dimensionslos ist. Diese Form ist 
der mathematische Ausdruck der Homogenität in den Dimensionen. 

: Die charakteristischen Größen unseres Problems sind die Aufsteigegeschwindigkeit 
w[LT"] der Blase, die Diehte [ML ®] der Flüssigkeit, die Erdbeschleunigung y|1,T ?]|-und 
der Radius R[Z] des Rohres. Wir haben »n—=4 Veränderliche und m =3 Grundeinheiten 
und können infolgedessen dureh das //-Theorem nur ein einziges (n - m=1) dimensionsloses 
Produkt haben, das gleich einer Konstante ist. 

Weil wir die Aufsteigegeschwindigkeit der Blase bestimmen wollen, setzen wir den 
Exponent dieser Größe Il und bekommen als dimensionsloses Produkt wo’ g°* R7T; setzen 
wir das Produkt gleich einer Konstante 4 und lösen nach w, so bekommen wir 


vr. 


Die rechte Seite muß die Dimension einer Geschwindigkeit haben. Setzen wir für Ö, g, R 
ihre Dimensionsformel ein, so erhalten wir 


[1 T-]= [MP [LT pr. 


Durch Vergleich der Exponenten von Z, M, T auf der linken und rechten Seite ergibt sich 
ein System von drei linearen Gleichungen für die «drei Unbekannten p, q, r. Die Lösung 


l 


5) 


lautet p=0, = Infolgedessen ist 


w=AyYyR. 


Also ist die Aufsteigegeschwindigkeit «der Blase proportional der Wurzel aus dem Radius des 
' . .. we-:, . n r 
Rohres und der Schwerebeschleunizung. 4° gt die sogenannte Froudesche Zahl und 
4 i 


gibt das Verhältnis der Trägheitskräfte zu den Schwerekräften. 
5) Buckingham: Phys. Rev. Bd. 4 (1914), S. 34. 


6, Bridgman: Theorie der physikalischen Dimensionen. Deutsche Ausgabe herausgegeben von H HIoell. 
Leipzig und Berlin 1932, 
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Die Aufsteigegeschwindigkeit der Blase hängt, wie zu erwarten, nicht von der Dichte ö 


acher 
ab; denn ist z. B. die Dichte doppelt so groß, so wird auch die auf jedes Element wirkende 


Fe Schwerkraft doppelt so groß, und weil die doppelte Kraft durch die doppelte Masse ihres 
Tände Elementes ausgeglichen wird, bleibt die Beschleunigung und infolgedessen auch die Ge- 
schwindigkeit unverändert. 

schen 2. Reihenansatz für das Potential. Weil unsere Strömung rotationssymmetrisch ist, so 
nicht müssen das Geschwindigkeitspotential und die Stokessche Stromfunktion die folgenden 
NS zu Gleichungen in Zylinderkoordinaten befriedigen: 

zwei IF 2a 

- ei Fr iz a dr 

einer —=yerfy 
stetie 


Die der Grenzbedingung für r—= R entsprechenden Lösungen dieser Glei- 
chungen lauten’): 








s auf ’ : 
fr ü 
men, 6GJ, (Bi 2) u: 7, 
ıan") \ 
dann ni 2 
r\-3 ; 
y=r Yy BEAT; 5) RLT, | 
„voll- zus 4 
h . e- | 
chtig BEN “ b Re 
olche wo ß; die Nullstellen der Besselschen Funktion erster Art und erster TE 
Form Ordnung urd e; willkürliche Konstanten bedeuten und wobei /’ eine | 
willkürliche Konstante ist, die auch =0 gesetzt werden kann. (Es ist | 
J,(B) = für d, =3,83171, A, = 7,01559 usw.) Wenn wir noch eine Parallel- == 
ellen strömung, deren Geschwindigkeit w, ist, überlagern, so bekommen wir RER ie 
h D oOort 1atenwanl,. 
Man N a 
a, rt  , 
Th; n)e . 
! 
n 
u 
nm ist Wo 2 F Y2. 2 
‚> y "+r Y cJ, (B: 5) 
uf | —— \ L, 
ekeit _. 
| und Wir machen die Formeln dimensionslos mit w,=AygR und R, wo wıe oben A eine dimen- 
eiten sionslose Konstante, g die Erdbeschleunigung und R den Radius des Rohres bedeuten. Die 
loses dimensionslosen Größen sind 
| o"= , w"- 7 u" = - "— = 
uno / AygR’’ AygR®’'’ iygR’ iygR’ 
etzen 


Rn ER 0 
R [) 5; i VgR®' 

Dann lauten das Geschwindigkeitspotential und die Stromfunktion, dimensionslos geschrieben, 
folgendermaßen: 


n 

1 = u _ 9 

g* Par ] + F Y k; Js (Pis)e ss 
i=1 


sich n 
I)SUNg E ” Ber: 
5 * u 7 
v7 1 \ k;J, (Pd) e ’ 
A — 


3 





Die Komponenten «* und ®* der Geschwindigkeit sind durch die Formel 
© 
n 


UTRRCOCL 


In 


3;n 


k; Pi J. (P; &) e er 
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gegeben. Wenn wir verlangen, daß &=n=0 ein Staupunkt sein soll («u=v=0), so folgt 
aus (2) wegen J,(O)=1 






n 
A => - 2 k; Pi " . . > x - . 5 s ° ® . ; . . (3). 


= 





Die Gl. (3) bestimmt die dimensionslose Größe A unseres Problems als Funktion der Kon- 
stanten k;. 


In der Einleitung wurde gezeigt, daß auf der Trennungsfläche: 


w”=u+v=2g3 





oder 
12 Bien. ® 
worllonl + ion) |=20R5. 
also 
> 
BEN En nn ee 


sein muß. Wenn wir in der Gl. (4) «* und v* durch ihre Ausdrücke (2) ersetzen unter Be- 
n 


2 k,ßi, bekommen wir nach Division der Stromfunktion durch & 
61 


=0 ist auch eine Stromlinie) die folgenden fundamentalen Gleichungen unseres Problems: 


n n 
"u = BR 1 3 Bi £ 
; Sr ) k; J, (Bi; &) e Bin we {. ) k; B;=0 ER re ac 


i=1 sel 


rücksiehtigung, daß = — 


n 2 . 
Fund, DL +2 Ih Bit Ft BiI (Bi) e A 89...» 
Die Gl. (5) gibt die Form der Blase, und die Gl. (6) bringt die noch zu erfüllende Rand- 
bedingung (u + v?—=2gz) zum Ausdruck. (Die weitere Randbedingung, wonach die Wände 
des Rohres Stromlinien sein sollen, ist bereits erfüllt.) Es handelt sich jetzt darum, die will- 
kürlichen Konstanten k; so zu bestimmen, daß die entsprechenden Werte &, n simultan die 
Gl. (5) und (6) befriedigen. 

Sind diese Konstanten bestimmt, so ergibt sich die Aufsteigegeschwindigkeit und die 
Form der Blase aus Gl. (3) bzw. (5). 





3. Bestimmung der Koeffizienten k;. Zur Bestimmung der Koeffizienten k; machen wir 
von den bekannten Reihenentwicklungen Gebrauch. 


991er et 
m #9=(5): i male) e+z1[2)® sale) + la) 
er 12a ler ale)" +)" 


Durch Einsetzung und Entwicklung nach Potenzen von 7 bekommen wir unter Einführung von 


2= 2 (2). 2= I u(d). 2= (6). 2-2 n(5) Mn) 


n n 
ER £ @. B. ;- Ri 
I BER IL n={ Sie taıgı "arg: +) 
i=1 i=| 
2 a 
Ti EI et 7 8), 
2? B.:, Bi . u’; 
+57 1B,— 121° a1381°°- 5141: + Nr 
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Nunsasem=e (5). (Bd) e Fish 








t=1 i—=1 
IE TE 
- Bi ta) 
9 
92 a B, Be B, L \ (9), 
Target rassıe 79 
Ba u 
+a7)Bı 8 191° 131° er In’ 
n n n ß n B 
N BI lid)en Pit b; kßi= 2 3 SAUER 2) 4) 
i=1 (= i=1 i=1 
WIE FRE 5 
— 2\ = "H+ayt rt an‘ + 
{Br Bıpy Berg Bear..), Pr 
ilı ey ey j" 
DE Te 
+alı Pete ayet [7 
Ei ee na KA re 
n=— TE 


die Gleichung für die Kontur der Blase mit &< 1 und positiven a,, also |n|= La, Pr 


Dann gehen wir mit diesem Ausdruck für n in die Gl. (8), (9), (10) ein, quadrieren und 
addieren darauf die beiden letzteren. Ordnen wir nun die beiden erhaltenen Beziehungen 
nach Potenzen von £, was wegen der absoluten Konvergenz der Potenzreihen erlaubt ist, so 
erhalten die fundamentalen Gl. (5) und (6) die folgende Gestalt: 


1 1 
D7 4 B,aa—B})&#®+ BE {B,—12B,aa+24B,a+24B,a,}&*+--:-=0. . (12), 
.. . ; 1 . j 
B#®+1-BB+4Ba+B3ö+  =-zta®+a,&+-.). .. 0. 0..(13). 
Die Bestimmung der Zahlen B,, B,, B, ... erfolgt durch Vergleich der. Koeffizienten von £ 


in den beiden Gl. (12) und (13). 
Ist a, beliebig gegeben, so folgt aus (13) 


B= 4 di; 
Mit diesem B, läßt sich B, aus (12) zu 
B,=4B,a, 
bestimmen. 
Sind B,, B,, a, bekannt, so folgt aus (13) mit beliebigen a, ' 
B.=- Bi+4Bial — 4B,B,a,+B -> 7 
T 2 


Haben wir auf diese Weise B,, B,... B„,. berechnet, so bestimmen sich die k; aus dem 
folgenden System von n linearen Gleichungen mit » Unbekannten (siehe (7)). 


x, (&) +4, (&) +.) —B, 
x, (&) +3, (6) +.) —B, 











. Mr u u . Z. angew. Math. Me 
144 Dumitresceu, Strömung an einer Luftblase im senkrechten Rohr Ba.2ı Nr. we 





Die. Determinante dieses Systems lautet bis auf einen konstanten Faktor: 
PB ße RX PB) Br BI —B)--- (Bi Pu) 

(B2 — 3) (Pr — Pr): - - (Ba — Bu) 

BB): (Pa — Pu) 


(Pr —1 Pn ); 


die wegen ;=ß; für alle i=#*%k von Null verschieden ist. Das System (14) besitzt also eine 
eindeutige Lösung. Die Integrationskonstanten k; sind also Funktionen der willkürlielen 
Koeffizienten a, und lassen sich eindeutig durch (14) bestimmen, wenn die Form der Blase 
vorgeschrieben ist, allerdings nur in einer genügend kleinen Umgebung des Nullpunktes. 


Nach dem Grenzübergang n >» lıandelt es sich um ein System von unendlich vielen 
Gleichungen mit wnendlich vielen Unbekannten, dessen eindeutige Auflösung uns möglich er- 
scheint. Jedenfalls wäre mit diesem Verfahren das Randwertproblem gelöst, eine Lösung der 
Differentialgleichung Ap=0 zu finden, so daß die Wände des Rolıres Stromlinien sind, 
Ö = 2 


; „| ist, wobei noch «a,, a, usw,, 


a ee 
und auf einer vorgegebenen Stromlinie 7: \+( A 
s / 


d. h. die Gestalt der Blase, vorgegeben werden kann. 


/ 


Es scheint also auf den ersten Blick, als ob das Problem unbestimmt wäre. Dies ist 
jedoch nicht der Fall, denn von allen nach dieser Betrachtung möglichen Lösungen wird nur 
eine das von uns vorausgesetzte asymptotische Verhalten für „> —» zeigen. Es gilt also 
eine asymptotische Lösung zu suchen und sie an die Potentiallösung anzuschließen. 

Der Verfasser war leider durch seine Einberufung zum Heer infolge des Kriegseintritts 
des rumänischen Staates gezwungen, an Stelle schärferer Rechnungen sich hierbei einer sehr 
rch“n Methode zu bedienen, die schnell arbeitet, dabei aber doch einen brauchbaren Zahlwert 
für den Eigenwert 4 liefert. 


. . v ” . ” . e 

4. Theoretische Bestimmung des Eigenwertes. Zum Zwecke der genäherten Ermittlung 

des Eigenwerts A sei der Krümmungskreis der Blase im Nullpunkt und die niedrigste asymp- 
totische Lösung als die erste Approximationskurve der Blase betrachtet. 


Sei o dieser Krümmungsradius, so läßt sich die Blasenkontur durch den Kreis 
e-n’ +0 


oder, da n für £==0 durchweg negativ ist, durch 


ı) s<o<I 


in der Nähe des Nullpunktes approximieren. Die photographischen Bilder (siehe den unten 
folgenden experimentellen Teil) werden zeigen, daß sich in der Tat die Blasenkuppe außer- 
ordentlich gut durch einen Kreis approximieren läßt, und hierauf beruht der Erfolg dieser 
ersten Näherung. 


‚ 2 
Nach der Entwicklung von V \ \ nach Potenzen von — hat die folgende Gestalt: 


PR TE hy 3(&) Bier, 3 5 
= eislc)+3 + el +1 . . . . . . . . (15), 


/ 


ist also von derselben Form wie der Ansatz (11), so daß sich die unbestimmten Koeffizienten dx 
als Funktionen des Krümmungsradius o durch die einfachen Bezieliungen 


1 1 


a, = r ad, = Erz 
2 8 0° N) 3 160° ’ 


(16) 


darstellen lassen. 


Sind die a; festgelegt, dann lassen sich die B;, (siehe Formeln 7), sodann die k; und schließ- 
lich A durch das System der linearen Gl. (14) bzw. durch die Formel (3) bestimmen, wie wir 
in dem vorstehenden Kapitel gezeigt haben. 
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06 Zahlentafel 1 gik. die Werte von a, Bx, A für ver- 

u | | schiedene Krümmungsradien o. In Bild 3 ist die Kurve Ä=f(o) 
| | | | aufgetragen. 

Te RE 





| | | | Zahlentafel 1. Berechnung des Eigenwertes A als Funktion des 
oc KARA Aa" Krümmungsradius im Scheitel 0. 





oe| a, Q, B, u: BB k, | ka | k; | 





eine | | | | | 05 1 1. 0,707 | 2,828 113,435 | — 0,0940 |+ 0,01678| — 0,02197) 0,460 


chen ; 
vr 0,6 | 0,833 | 0,579 | 0,645 2,151 | 8,516 | — 0,0832 | — 0,00839 — 0,00915| 0,471 


ktes. ig 0,7 | 0,714 | 0,364 | 0,597 | 1,707 | 5,793 | — 0,0962 | — 0,01097|— 0,00423| 0,488 


Bild 3. Der Eigenwert 4 als 


ielen Funktion des Krümmungsradius 0,8 | 0,625 0,244 0,559 | 1,397 | 4,148 — 0,1140 | — 0,00613)— 0,00251| 0,505 


or im Scheitel ®. 

; der Um andererseits eine erste Näherung einer asymptotischen Darstellung der Kurve zu be- 
a D’w 
4 
BEW,, Blase für jeden Querschnitt dieselbe ist. Macht man weiter die Annahme, daß die Ge- 
schwindigkeit im unteren Teil des Rohres nahezu senkrecht gerichtet und für jeden festen 
Wert von z genähert den konstanten Betrag w=y2gz hat, dann erhält man als Ausdruck 

für die Kontinuitätsgleichung: 


sind, 


ao 


kommen, verlangen wir, daß die Durchflußmenge @ — zwischen Rohrwandung und 


s ist 
nur au - 
also TU2RY”— er? y2gz=T (2RriAygR 
mit R=Radius des Rohres und w,—=4/ygR oder dimensionslos gemacht 
ji? 
u er 


tritts 
sehr 
wert 


(17). 


Weil ä durch das oben dargestellte Verfahren eine Funktion des Krümmungsradius,o ist, so 
entspricht jedem Krümmungskreis vom Radius o, d. h. jeder Lösung in der Umgebung des 
Nullpunktes, eine asymptotische Lösung. Der Krümmungsradius o bzw. die Dimensionslose A 
ist also der Parameter einer einparametrigen Schar von Kurvenpaaren. Zur eindeutigen Be- 
stimmung dieses Parameters A fordern wir, daß die Kurven (15) und (17) in ihrem Schnittpunkt 


lung 
ymp- 





a ] Bild 4a (links). „ 
| | Asymptotische Ge- 

+ — stalt der Blase für 
den Krümmungs- 
radius im Scheitel 
0=0,7. 

Bild 4b (rechts). 

Asymptotische Ge- / 

stalt der Blase für 7 





2 / den Krümmungs- 
inten | |  radius im Scheitel 
ußer- k ) 1 } 0e—=0,8. 37) 05 — 1 15 


ieser 














dieselben Tangenten besitzen. Die Lösung wurde für 
oe=0,; 0,6; 0,7 und 0,8 und die zugehörigen A-Werte 
versucht. Die Bilder 4a, 4b zeigen die verschiedenen 
Lagen von Krümmungskreisen und entsprechenden 
asymptotischen Lösungen. 

Wie aus den Bildern 4a und 4b zu ersehen ist, 
muß offenbar die Lösung zwischen o=0,7 und o=0,8 
gelegen sein. Die Tangente der beiden Kurven ist durch 

dn 


dn_tnE dp 
u“ i-P 2” 





stalt: 

















(15), 





n dk 








(16) m 05 —) 
Bild 4e. Asymptotische Gestalt der Blase für r i f ä 
den Krümmungsradius im Scheitel @=0,75. Der Krümmungskreis o = 0,75 schneidet sich mit 
seiner entsprechenden asymptotischen Lösung im Punkte 
£=0,71, n= 0,50, und die beiden Kurven besitzen in diesem Punkte eine gemeinsame Tan- 
gente (Bild 4c). Sie bestimmen somit eindeutig die Lösung unseres vereinfachten Problems. 
Der entsprechende Eigenwert ist A = 0,4%. 


gegeben. 


10 
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III. Experimenteller Teil. 


1. Versuchsanordnung. Zur Durchführung der Versuche*®) wurden vier Glasröhren von 
0,99 em, 2 em, 3,76 em und 7,0 cm Durchmesser und etwa 120 cm Länge verwendet. 


Die Glasröhren d (s. Bild 5), welche ie > 
unten glatt abgeschnitten waren, besaßen gr 
am oberen Ende als Abschluß einen Drei- 
wegehahn. Die Röhren wurden im Be- 
obachtungsfeld von einer Glasküvette c AR, DE 
mit planparallelen Wänden umgeben, um >) 
verzerrungsfreie Aufnahmen machen zu 
können. Die Küvette, welche ein Blick- Bi 
feld von 16%X 42 cm ermöglichte, war 
mit Wasser gefüllt. 

Der Unterteil der Röhren ragte u; 
etwa 4 cm in einen mit Wasser gefüllten Tr Te TI TTTTE TEE 
Behälter e. Die zur Bildung der Luftblase uns) Vemuniumssrisehe- 
benötigte Luft wurde unter Zwischen- 
schaltung eines Vorratsbehälters f aus einer Preßluftleitung entnommen und durch ein U- 
förmig gebogenes Glasrohr in das Versuchsrohr von unten eingeleitet. Auf diese Weise war 
es möglich, Luftblasen von konstanter Länge zu erzeugen. 


> 


—— Druckluft 














Durchführung der Messungen. Die Aufsteigegeschwindigkeit und die Form der 
Blase wurden auf photographischem Wege bestimmt. Zur Beleuchtung diente eine Bogen- 
lampe mit Kondensor, deren Licht auf eine an der Glasküvette befestigte Mattscheibe g fiel. 
Die Mattscheibe war nötig, um ein gleichmäßiges diffuses Licht zu erhalten. Außerdem 
wurden Versuche zur Belichtung mit elektrischen Funken gemacht, die jedoch zu keinen be- 
friedigenden Ergebnissen führten. Als Aufnahmeapparat diente die Leica IIIA, welche Be- 
liehtungszeiten bis "/ıooo Sekunde und rasche Bildfolge zuließ. 


Zur Zeitregistrierung wurde ein Synchronmotor h verwendet, auf dessen Achse ein 
Zeiger befestigt war, der in 0,6 Sekunden eine volle Umdrehung machte. Der Zeiger wurde 
mit einer zweiten Bogenlampe angestrahlt. Im Inneren der Küvette nahe dem Versuchsrohr 
waren zwei Marken im Abstand von 10 cm als Maßstab angebracht und außerdem zwei Paß- 
marken außen an der Küvette. Zur Durchführung eines Versuches wurde das Versuchsrohr, 
nachdem es genau senkrecht ausgerichtet war, mit Hilfe einer Wasserstrahlpumpe aus dem 
Behälter luftfrei mit Wasser gefüllt. Anschließend wurde Preßluft in den Vorratsbehälter 
eingeleitet und dieser abgesperrt. 

Nachdem der Synchronmotor angeworfen war, wurde der Absperrhahn rasch geöffnet, 
so daß die Preßluft durch das U-Rohr in das Versuchsrohr einströmen konnte und hier die 
Blase erzeugte. Sobald die Blase im Blickfeld der Kamera erschien, wurde die erste und 
kurz vor Verlassen des Blickfeldes die zweite Aufnahme mit je !Jıooo Sekunde Belichtungszeit 
gemacht. Für die Bestimmung der Form der Blase wurde der Apparat nähergebracht, um 
eine größere Abbildung zu bekommen. 


Auswertung. Ausgewertet wurden die Aufnahmen in einem Vergrößerungsapparat. 
Es wurden auf einem Blatt Papier die Kuppe der Blase, die Zeigerstellung der Zeituhr und 
die Paßmarken eingezeichnet. Dann wurden die Paßmarken der zweiten Aufnahme mit denen 
der ersten zur Deckung gebracht und hier ebenfalls Blasenkuppe und Zeigerstellungen ein- 
gezeichnet. Außerdem wurde der Maßstab markiert. Nun wurde die Strecke zwischen der 
Blasenkuppe der ersten und zweiten Aufnahme gemessen und durch die Länge des Maßstabes 
dividiert und so die wirkliche Länge erhalten. Ebenso wurde der Winkel zwischen der ersten 


und der zweiten Zeigerstellung der Uhr gemessen und auf diese Weise die Zeit festgestellt, - 


die zwischen den beiden Aufnahmen verstrichen war. 


Um einen Irrtum um ganze Zeigerumdrehungen zu vermeiden, wurde die Zeit außerdem 
mit einer Stoppuhr festgehalten. Aus dem zurückgelegten Weg der Blase und der ver- 
strichenen Zeit wurde die Aufsteigegeschwindigkeit der Blase bestimmt. 


2. Ergebnisse und Vergleich mit der Theorie. Die Versuche haben bestätigt, daß die Auf- 
steigegeschwindigkeit der Blase bet den weiten Rohren proportional der Wurzel aus dem 
Radius R des Rohres und die Schwerebeschleunigung g ist, wie sich im ersten Teil dieser 


8) Durch die Kriegslage war ich genötigt, die Versuche im Physikalischen Institut Göttingen auszuführen. Herrn 
Professor Dr. R. W. Pohl, der mir die Möglichkeit hierzu gegeben hat, bin ich zu großem Dank verpflichtet. 











Nr 


[——— 


in Aa a AM A Me a Di 


u 


u 


. Math. Mech. . a; a ’ 
ee 3 Juni 1943 Dufitrescu, Strömung an einer Luftblase im senkrechten Rohr 





Zahlentafel 2a. Zahlentafel 2b. 


Geschwindigkeitsmessungen im Rohr mit D=0,99cm. Geschwindigkeitsmessungen im Rohr mit D=2,0 cm. 





D= 0,99 cm Wm = 6,316 cm/sec D==2,0 cm wm = 14,85 em/sec 


IE En. a ee. 4% u a" 


2 
sec sec ] em ‚ em’sec sec | sec } em | em/s2c 


“ 
pe 





| 226 : | 0,331 | 0,198 | 4,398 | 27,39 | 6,227 

| 236,2 | 0,609 | 0,365 | 4,562 | 28,63 | 6,271 
238,7 | 0,675 | 0.405 | 4,605 | 28,93 | 6,282 
245,5 | 0,872 | 0,523 | 4,723 | 29,75 | 6,298 
247,0 | 0,917 | 0,550 | 4,755 | 29,94 | 6,296 
251,4 | 0,037 | 0,0222) 4,822 | 30,47 | 6,318 
252,3 | 0,063 | 0,0378) 4,837 | 30,58 | 6,322 
253,4 | 0,085 | 0,051 | 4,851 | 30,71 | 6,330 

9 | 268,5 | 0,370 | 0,222 | 5,022 | 31,94 | 6,360 
10 | 263,2 | 0,360 | 0,216 | 5,016 | 31,90 | 6,360 
11 | 266,4 | 0,459 | 0,275 | 5,075 | 32,29 | 6,362 
12 | 267,2 | 0,471 | 0,283 | 5,083 | 32,38 | 6,370 


0,396 | 0,237 | 1,437 | 214,30 | 14,91 
' 0,514 , 0,308 | 1,508 | 224,48 | 14,88 
0,67 | 0,402 | 1,602 | 236,96 | 14,79 
| 0,67 | 0,402 | 1,602 |237,57 | 14,83 
0,78 | 0,468 | 1,668 |248,48| 14,89 
| 0,927 | 0,556 | 1,756 | 259,15 | 14,76 
| 0,018 | 0,0108| 1,8108 269,33 | 14,87 
' 0,081 | 0,0486| 1,8486| 274,90 | 14,87 
' 0,123 | 0,0738| 1,873 |278,90 | 14,88 
| 0,274 | 0,164 | 1,964 | 291,63 | 14,85 
| 0,319 | 0,191 | 1,991 296,0 | 14,86 


& 1 mn Da wm DD 


a  n 
0m 
u 


” 
s bedeuten: 


E 

D Durchmesser des Rohres (cm), 

l Scheitelabstand der beiden Blasenaufnahmen in der Vergrößerung (mm). Die wahre Länge 
für 100 mm ist 82,5=m, 

t, Zeigerstellung der Uhr (Bruchteile der vollen Umdrehung), 

t,—=0$61, Zeitintervall, das ?, entspricht (sec), 

T,=t+n0,6 Zeit zwischen den beiden Aufnahmen (sec) (n, Anzahl der vollen Umdrehungen des 
Uhrzeigers), 

L=I!/m von der Blase in der Zeit T zurückgelegte Strecke (cm), 

v=L/T Steiggeschwindigkeit der Blase (cm/sec), 

Om Mittelwert der Aufsteiggeschwindigkeit der Blase. 


Zahlentafel 2c. Zahlentafel 2d. 
Geschwindigkeitsmessungen im Rohr mit D=3,76cm. Geschwindigkeitsmessungen im Rohr mit D=1,0 cm. 





D=7,0cm Um = 28,64 cm/sec 


=. BERKER” 


sec x | sec 


158,2 | 0,447 | 1,048 | 21,17 179,8 | 0,271 | 0,162 | 0,762 | 21,793 | 28,59 
163,7 0,490 | 1,090 | 22,959 | 21,06 213,5 | 0,500 | 0,300 | 0,900 | 25,878 | 28,75 
177,5 0,576 | 1,176 | 24,894 | 21,17 214,3 | 0,522 | 0,318 | 0,913 | 25,975 | 28,45 
179,0 0,584 | 1,184 | 25,105 | 21,20 216,7 | 0,535 | 0,321 | 0,921 | 26,266 | 28,52 
179,8 0,584 | 1,184 |25,147| 21,24 223,7 | 0,585 | 0,348 | 0,948 | 27,115 | 28,60 
181,2 0,0 | 1,200 |25,413| 21,18 224,2 | 0,584 | 0,350 | 0,950 | 27,175 | 28,60 
187,0 0,030 | 1,230 | 26,227 | 21,32 225,7 | 0,585 | 0,351 | 0,951 | 27,357 | 28,76 
192,2 0,072 | 1,272 | 26,956 | 21,19 226,2 | 0,581 | 0,348 | 0,948 | 27,418 | 28,92 
205,8 0,162 | 1,362 | 28,863 | 21,19 226,5 | 0,602 | 0,361 | 0,961 | 27,154 | 28,57 
215,0 | 0,217 | 1,417 | 30,154 | 21,28 ı0 | 228,3 | 0,622 | 0,373 | 0,973 |27,672| 28,44 
218,5 0,231 | 1,431 | 30,645 | 21,41 11 | 241,1 | 0,690 | 0,414 | 1,014 | 27,224 | 28,82 

12 | 246,4 | 0,738 | 0,439 | 1,039 | 29,866 | 28,74 

13 | 256,7 | 0,812 | 0,487 | 1,087 |31,115| 38,62 

14 | 256,8 | 0,816 | 0,489 | 1,089 | 31,127 | 28,58 


sec cm , em/sec 





ww 


X 1 Dr DD 





oo 
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Zahlentafel 3 (links). Mittelwerte 
der dimensionslosen Abszissen & und 
ÖOrdinaten n nach den Messungen. 











Zahlentafel4 (rechts). DieWerte des 
gemessenen Eigenwertes A als Funktion 
des Durchmessers des Rohres und der 


dimensionslosen Zahl P=D 2 
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Arbeit durch Dimensionsbetrachtung 

ergeben hat. Die Länge der Blase 

war ohne Einfluß. Bei den engeren 

Rohren bewirken Kapillarkräfte eine 

Verringung der Aufstieggeschwindig- 

keit, Die Zahlentafeln 2a bis d geben : 
die direkten Messungen wieder sowie 

ihre Auswertung nach dem im zweiten 

Kapitel angegebenen Verfahren. Die 

Form der Blase ist aus den Bildern 6 

und 7 zu ersehen. Die Bilder 6a, b, c,d 

stellen die Röhre in halber Größe dar; 

die Bilder 7a, b, c, d dieselben Röhren | 


auf gleichen Maßstab gebracht. 

Die Form der Blase ist fast ge- 
nau dieselbe für alle vier verschie- 
denen Röhren. Für das Rohr mit 
0,399 «m Durchmesser hat die Blase 
eine spitzere Kuppe, was vermutlich 
dem hier bereits merklichen Einfluß 
der Kapillarkräfte zuzuschreiben ist?). 
Die Zahlentafel 3 gibt die Mittel- 
werte der dimensionslosen Abszissen » | 
und Ordinaten £ der Blasen a und b. | | ' 








Bild 6. Blasen, in halber Größe wiedergegeben. 

















F 


Die Werte des gemessenen Eigen- | 


wertes A als Funktion der dem 
Durchmesser des ‚Rohres D Br Bild 7. Dieselben Blasen wie in Bild 6, 
tionalen dimensionslosen Zahl '°) jedoch auf gleichen Maßstab gebracht. 











—_.__ berechnete Kurve 
— erperimentelle ” 























se 05 Ze 15 2 


Bild 8. Die Werte des gemessenen Eigen- Bild 9. Berechnete und experimentell gefundene Blasenform. 
wertes 4 als Funktion der dimensionslosen 


77) 
Zahl P=D| - 


P=D v2 (D= Durchmesser des Rohres, g = 981 cm/see Erdbeschleunigung, ö-1gr/cm? 


Dichte des Wassers uind C= 70 dyn/cm Kapillarkonstante) sind in Bild 8 aufgetragen (Zahlen- 
tafel 4). Der Einfluß der Kapillar- und Zähigkeitskräfte macht sich bis zu einem Durch- 
messer D>=3cm bemerkbar. Von hier an strebt A asymptotisch gegen den Wert Ä=0,5 in 
guter Übereinstimmung mit unserem theoretisch gefundenen Wert A=0,496, den wir ohne 
Berücksichtigung der Kapillar- und Zähigkeitskräfte gefunden haben. 

Die experimentelle und die gerechnete Näherung der theoretischen Form der Blase ist 
in Bild 9 ausgezogen bzw. gestrichelt dargestellt. Der Berührtängspunkt des Krümmungs- 
kreises und der asymptotischen Lösung liegt bei £=0,71, 7 = 0,50. 

Wie bereits oben erwähnt, beruht der Erfolg dieser Methode, welche uns nur die erste 
Näherung des Eigenwertes A lieferte, auf der Tatsache, daß sich die Blasenkuppe in der Um- 
gebung des Nullpunktes sehr gut durch einen Kreis approximieren läßt (für alle &< 0,6, < 0,3). 

Wegen der geringen Abweichungen der experimentellen von der theoretischen Kurve lolınt 
sich der große Aufwand nicht, die erste Näherung durch eine zweite weiter zu verbessern. 


®) Möglicherweise ist auch die Wirkung des Glasrohrs als Zylinderlinse mit im Spiel, die bei dem engsten Rohr 
sich am stärksten bemerklich macht. 
10) Vgl. z. B. L. Prandtl: Führer durch die Strömungslehre, 111. Aufl. Braunschweig 1942, I. Abschnitt, $ 11. 
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IV. Zusammenfassung. 


Die Aufgabe, die Form und die Aufsteiggeschwindigk@it einer unendlich langen Luft- 
blase, die in einem vertikalen flüssigkeitsgefüllten Kreisrohr aufsteigt, zu finden, führt unter 
Vernachlässigung der Kapillar- und Zähigkeitskräfte zur Lösung der linearen Differential- 
gleichung zweiter Ordnung Ap=0 für das Geschwindigkeitspotential ünter einer besonderen 
Grenzbedingung. 

Zunächst wurde durch Dimensionsbetrachtungen gezeigt, daß die Aufsteiggeschwindigkeit 
der Blase w=AygR ist, wo g die Erdbeschleunigung, $ der Radius des Rohres und Ä eine 
dimensionslose Konstante bedeutet. 

Die Lösung dieser achsensymmetrischen Strömung ist mit £ und n als mit R dimensions- 
los gemachten Koordinaten durch die Ausdrücke 


n n 
1 > g —ßs;n . E E” u u 
ym— 147 > kIBHEe" und y=7 +7 2 KI(He fi" 


i=1 6m 1 


gegeben, welche durch die Bedingung J, (#;)=0 die erste Randbedingung, wonach die Wände 
des Rohres Stromlinien sein sollen, von selbst erfüllen. Die zweite Randbedingung, daß auf 


. a u 0; - 
einer noch nicht bekannten Stromlinie 52) + Tr -=7|n| sein soll, wurde in folgender 
Weise gelöst. Die unbekannte Form der Blase wurde in der Nähe des Nullpunktes durch 
ihren entsprechenden Krümmungskreis dargestellt. Die Reihenentwicklung dieses Kreises 
n 


2 ; N op\’ ... ’ 
n= 2 a; &* wurde in die Stromfunktion sowie in den Ausdruck + pn -;|n| ein- 
1 a 08 07 ' Eh 


gesetzt und auf Grund einiger Rechnungen wurden. die Integrationskonstanten k; als Funktion 


n 
des Kreisradius bestimmt. Der Eigenwert A läßt sich durch die Formel = —  k; ß; auch 
i=1 
unmittelbar als Funktion des Kreisradius darstellen. 

Die dimensionslose Konstante A (Eigenwert des Problems) wurde dann durch die 
Forderung bestimmt, daß die asymptotische Lösung mit stetiger Tangente in den ihr zu- 
geordneten Krümmungskreis übergehen muß. 

Der theoretisch bestimmte Wert 4=0,4% stimmt mit dem durch die Versuche ge- 
fundenen Wert A=0,49, ebenso die theoretische und die experimentelle Kurve der Blase gut 
überein. 437 


Über die Stabilität der Lösungen Hillscher Differential- 
gleichungen mit drei unabhängigen Parametern. 
Erste Mitteilung: 

Über die Gleihung y"+ (A+y,cosz+7,cos2z)y=0. 


Von Karl Klotter und Gertrud Kotowski in Berlin. 


(Mitteilung aus dem Vierjahresplan-Institut für Schwingungsforschung, Berlin.) 


Das Stabilitätsverhalten der Lösungen der. genannten Differentialgleichung 
wird dadurch untersucht, daß der Stabilitätskörper im A-y,-y,- Raum 
bestimmt wird. Zur Darstellung dienen Schnitte y,—=const durch diesen 
Körper. Der Schnitt y,=0 stellt die „Struttische Karte“ für die 
Mathieusche Differentialgleichung dar. 


1. Einleitung. 


Seit den Differentialgleichungen mit periodischen Koeffizienten eine immer größere Be- 
deutung auch bei den technischen Anwendungen zukommt [1]'), ist die Gewinnung syste- 
matischer numerischer Ergebnisse (z. B. die Tabulierung zusammengehöriger Werte für die 
Parameter der Differentialgleichung, die die Grenzen der stabilen und instabilen Gebiete der 
Lösungen bezeichnen) zu einem dringenden Bedürfnis geworden. Denn die Durchführung der 
Rechnungen ist, selbst wenn, wie bei der Hillschen oder Mathieuschen Differentialgleichung, 
die Methode zur Bestimmung der Lösungen längst klargestellt ist, im allgemeinen recht zeit- 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. 
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raubend, so daß für eine Benutzung Zalılenmaterial bereitgestellt werden muß. Das vor- 
handene numerische Material ist aber leider noch sehr lückenhaft. Soweit uns bekannt ist, 
sind numerische Untersuchungef systematischer Art bis jetzt nur für folgende zweipara- 
metrige Differentialgleichungen durchgeführt worden: 


a) die Mathieusche Differentialgleichung [2], [4] 
yY"+(i+yeose)y=0, 

b) die Meißnersche Differentialgleichung [3] 
yY"+lkty)y=0, 

c) die spezielle Hillsche Differentialgleichung [4] 


yY'(A+y c0s2)+ 71 =6, 


Die Rechnungen in [2] sind dabei die ausführlichsten; sie enthalten außer den Tabellen 
für die Wertepaare A und y der Grenzkurven auch Tabellen für die Fourierkoeffizienten 
der Mathieuschen Funktionen erster Art und für diese Funktionen selber. 


2. Aufgabenstellung. 


Die hier vorliegende Arbeit soll einen Beitrag darstellen zur Bereitstellung des wichtig- 
sten Zahlenmaterials. Es wurde der Einfluß eines dritten unabhängigen Parameters in der 
Hillschen Differentialgleichung auf die Stabilität der Lösungen untersucht. Die hier be. 
handelte spezielle Hillsche Differentialgleichung lautet 


rar, esty,sdyr0 . . . . 2.2... 


Berechnet wurden für die Werte 
y,=0; 05; 1,0 und 1,5 
und 
= +01; +02; +03; +04 und +05 
die Grenzkurven 
Cu Cm In, €, und S, 


(in der Bezeichnungsweise des Aufsatzes [1]). 

Trägt man die Parameter 4, y,, y, in einem dreiachsigen kartesischen Koordinatensystem 
auf, so erhält man den „Stabilitätskörper“. Zusammengehörige Werte der drei Parameter, 
die die Grenzen zwischen den Bereichen stabiler und instabiler Lösungen bezeichnen, bilden 
Flächen in diesem Stabilitätskörper. Wir wählen zur Darstellung Schnitte durch den Körper, 
die parallel zur A-y,-Ebene sind, also y, als Parameter haben; in diesem Sinne sprechen 
wir von Grenzkurven statt von Grenzflächen. Die bekannte „Struttsche Karte“, die 
über die Stabilität der Lösungen der Mathieuschen Differentialgleichung Auskunft gibt, ist 
also der Schnitt durch den Körper der 4-,y,-,y,-Werte für y,—=0. 

Über die Stabilitätsverhältnisse anderer Hillscher Differentialgleichungen von ähnlicher 
Bauart, z. B. der N 

+(+y,cos2+y,siu2x)y=0 
soll in weiteren FARENNEN, ee werden. 


3. Rechnungsgang. 

Da die Belegungsfunktion (A+y,cosx=-+y,cos2x) beschränkt ist, kann man die Hill- 

sche Lösungsmethode anwenden. Zu den Grenzkurven gehören periodische Fundamental- 
lösungen von der Form (n sei eine ganze, positive Zahl) 


an 


a) für die C„-Kurven: y=} „Ar cosvx 


on 


3) für die S,„-Kurven: y =2’B 


on 
y) für die Can +1-Kurven: y= As, + 1c0s 


‚= 0 9° 


5) 
2v+1 
5 x 


” 


6) für die Syn +1-Kurven: y= 2, Bart sin? eg 1 
a ‚—_( 





2 


. Die unendlichen Systeme algebraischer Gleichungen, die man zur Bestimmung ihrer Fourier- 
koeffizienten A und B erhält, stehen auf Formelblatt 1 (GIn. (3)). 
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Formelblatt 1. 
a) für die C„-Kurven 


1A,+ A ‚+34 
A a )a+%4 44% 4, 


5A+G-9A, +5 A, +34 


(Ak) Art (An + Ars) + (Art Ars) 


P) für die S„-Kurven 


” 


inerona dann 


3B,+@—-9B,+5B,+%B, 


(dk) B, +4 [Br + Br] + 2 [Be „+ Brk+:l 


y) für die Cy„ + |-Kurven 


2 


1 
142] ee 
1-7|4 +3 3 Ay 


2 A,+% Ay+la-7 Ay + GA +5 As, 


9 
en Fa 3 Bs,+ Bun 


25 ’, y 
: Bs,+ 2 B: + 3 B» 2 


+1V j 
\\Burı+% Ban + Bun 4a +3 [BuatBurs 
2 2 2 2 


3 [Au Hs +3 Aust Ass = 
2 2 


=0 
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Die gleich Null gesetzten Determinanten dieser Gleichungssysteme (3) geben die 


Gleichungen 
Ca (4, 71, Y)=0, Sul, yo Y)=0, 


Cyn+1(d, Yıy9)=0 und Syn+1 (Yo) = 0 


D} 


der Grenzkurven. Aus diesen Gleichungen erkennt man, daß (genau wie bei der Mathieu- 
schen Differentialgleichung) die Funktionen C„ und S„ gerade Funktionen hinsichtlich des 
Wertes y, sind, während die Funktionen C,,,;, und S;,„ + bei einem Wechsel des Vorzeichens 


z = 
von y, ihre Rollen vertauschen. Wir brauchen deshalb die Rechnungen nur für positive 
Werte von y, durchzuführen. 

Für die explizite Berechnung der zu den Grenzkurven gehörenden Wertetripel A, y,, y, 
wurden jedoch nicht die sich aus dem Nullsetzen der Systemdeterminanten ergebenden Glei- 
chungen unmittelbar herangezogen, sondern es wurde ein Annäherungsverfahren angewendet, 
wie es auch von F. Kluge [5, S. 135] benutzt wurde. Die Gleichungssysteme (3) wurden hierbei 
jeweils nach der fünften Gleichung abgebrochen. So entstehen Systeme von jeweils fünf homo- 
genen Gleichungen mit fünf Unbekannten. Eine der Unbekannten — hier stets die letzte 
wurde zu Eins normiert. Das Gleichungssystem wird dadurch zu einem inhomogenen System für 
vier Unbekannte; eine Gleichung ist überzählig; sie dient als „Fehlergleichung“. Im folgen- 
den wurde stets die erste Gleichung als Fehlergleichung benutzt. Für die Größen y, und y, 
wurden die in Abschnitt 2 genannten Wertepaare eingesetzt; für die’eigentlich gesuchte, bei 
diesem Annäherungsverfahren gleichfalls als Parameter behandelte Größe A wurden An- 
näherungswerte eingesetzt, die nach folgendem Gesichtspunkt anfgesucht wurden: Da die 
Werte A für die gewählten Werte von y, bei dem Wert y,=0 aus den Inceschen Rech- 
nungen [2] bekannt sind, so wählten wir in der Nähe dieser Werte A die Annäherungswerte 
für „= +0,11; in der Nähe der so erhaltenen Werte A wieder die Annäherungswerte für 
,=+023 usw. z 

Das folgende Beispiel zeigt die Bestimmung von A für das Wertepaar y,=1,5 und 

=0,5 der Grenzkurve S,—=0 nach diesem Verfahren. Das inhomogene Gleichungssystem 
lautet 
05B, +05B, =B3-4 


4+025B, +05B, +@-16)B,=--0,75 
025B,+075B, +@-9NB,+05B = " 
0,75B,+(A—4)B,+075B, +025B, = 0 
und die „Fehlergleichung“ 
BIER HBR ST . .: 5 222. 
Als angenäherte Werte Ä setzen wir hierin (da das zu y,=1,5 und y,=0,4 gehörige A den 


Wert 0,99 hatte, und mit wachsendem y, die A-Werte größer werden, und zwar um ungefähr 
0,04 für einen Zuwachs in y, von 0,1): 


4=1,03; 1,04 und 1,05 
Dazu gehören die folgenden Fehlerwerte f nach (4b) 
) % D 
5, 4 | 1,03 | 1,04 | 1,05 


Y3s, Bra | fr | 6,0076 | +6,9495 | + 19,8203 
(5 Gleichungen) | 








In Bild 1 ist die „Fehlergerade“ f(A) eingezeichnet; ihr Schnitt- 
punkt mit der 4-Achse liegt bei A = 1,0346. 

Um die Güte der so erzielten Annäherung zu erkennen, 
wurde für den höchsten y,-Wert und die ng dem Betrage 
nach höchsten y,-Werte, also für y„=15 und „=+0,4; +05 
der Wert von Ä auch aus einem System von us Gleichungen 
ermittelt. In der Zahlentafel 1 sind jeweils in der Spalte A, die- 
jenigen Werte von 4 aufgezeichnet, die sich bei Benutzung von 
fünf Gleichungen, in der Spalte A, jene, die sich bei Benutzung 
von sechs Gleichungen ergaben. Wie man aus der Zahlentafek1 
ersieht, ist die dritte Stelle nach dem Komma noch als sicher zu 
betrachten. Bezüglich der Rechengenauigkeit bemerken wir: Im 
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Zahlentafel 1. 
9 = — 0,4 


D) 





5, C, Co Cur2 Sur 5; C, 





0.8310 | — 0,82 7518 | 0,9934 | 1,4468 hs | — 0.6232 0,6065 | 0,4813 | 0,6326 | 1,5626 
— 0.8305 2 0,9933 | 1,4468 a ı — 0,6233 | — 0.6065 | 0,4815 | 0,6327 1.5624 


„,=1B; 


PA 





& S, c, R, ki ü, S, c, 





E 2,8660 | — 0,8580 | 0,7775 | 1,0345 | 1,4277 ls | — 0,6085 | — 0,5904 | 0,4411 | 0,5843 | 1,5740 
is — 0.8659 | — 0,8580 | 0,7777 | 1,0346 | 1,4277 hg \ — 0,6085 | — 0,5904 | 0,4411 | 0,5842 | 1,5740 


Verlauf der Rechnung (auf der Rechenmaschine) wurden sämtliche Zahlenwerte bis auf die 
sechste Stelle nach dem Komma mitgenommen. Das Resultat wurde auf die vierte Stelle 
nach dem Komma gekürzt. In den Fehlerdiagrammen (nach Bild 1) wurde der Maßstab für A 
so gewählt, daß der Wert 0,01 durch eine Strecke von mindestens 5 cm Länge dargestellt 
wurde, so daß stets die dritte Stelle nach dem Komma noch genau und die vierte Stelie 
geschätzt abgelesen werden konnte. 


4. Ergebnisse und ihre Erörterung. 


Die auf die oben beschriebene Weise erzielten Ergebnisse stehen in den Zahlentafeln 2 
und 3. Die Bilder 2a bis e zeigen Photographien nach einem Modell des Stabilitätskörpers. 





In Bild 3 und Bild 4 sind je zwei Schnitte 
durch diesen Körper parallel zur 4-y,- Ebene, 
und zwar für y,=0 (übliche Struttsche 
Karte) und y,=0,5 bzw. y„„=—0,5 überein- 
ander gezeichnet. Die Schnitte für die da- 
zwischenliegenden Werte von y, zeigen ein 
entsprechendes Verhalten, das sich mit größer 
werdenden |y,| immer deutlicher ausprägt. Be- 
merkenswert an Bild 3 ist die Überschneidung 
der Kurven C, und $, bei positiven Werten 
von Y,. Der Schnittpünkt der beiden Kurven 
wandert mit wachsendem y, von y,=0 (wo 
er für y„—=0 liegt) nach höheren Werten von y,. 
Bild 4 dagegen zeigt, daß für negatives y, 
überhaupt kein Schnittpunkt der Kurven C, 
und $, auftritt; diese Kurven rücken vielmehr 
Bild 2b. mit wachsendem |y,|immer weiter auseinander. 
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Zahlentafel 2. 4-Werte der Grenzkurven. 





01 | 02 | 0,3 0,4 05 





— 0,0013 | — 0,0051 | — 0,0112 | —0,019 | — 0,0811 
— 0,1224 | — 0,1338 | — 0,1480 | — 0,1650 | — 0,1840 
— 0,390 | —04172 | — 045 | — 0,4701 | — 0,5000 
| 0,7359 | — 0,7650 | —0,7960 | — 0,8310 | — 0,8660 





0,2432 | 02436 | 0,2354 0,2236 0,2088 
— 0,0398 | —0,0562 | — 0,0748 | —0,0%2 | — 0,1203 
— 0,3691 | — 0,3931 | — 0,4498 | — 04486 | — 0,4795 
— 0,1550 | —0785 | - 0,8217 | — 0,8580 











0,2436 | 0,2354 0,2236 0,2088 
0,4848 | 0,4896 0,4911 0,4892 
0,6411 | 0,6602 0,6763 | 0,6892 
0,6947 0,7242 0,7518 0,7775 





0,8987 | 0,8472 | 0,7951 | 07424 
1,0088 | 0,9662 | 0,9232 | 0,8802 
1,2369 | 1,209 | 1,1778 | 1,1588 
1,4814 | 1,4656 | 1,4468 | 1,4277 








1,098 | 1,472 1,1949 1,2419 
1,0767 | 1,1245 1,1716 1,2179 
1.010 | 1,0581 | 1,1081 | 1,1473 
0,9078 0,9509 | 0,9984 | 1,0846 











Zahlentafel 3. A-Werte der Grenzkurven. 


— 0,1 | —02 018 | —04 —05 








— 0,0013 | 0,019 | — 0,0311 


— 0,1076 — 0,1074 | — 0,1129 
— 0,3619 — 0,3286 | — 0,3235 
— 0,6637 — 0,6232 | — 0,6085 





0,2482 | 02236 | 0,2088 
— 0,0178 | — 0,0072 | — 0,0102 
— 0,3287 02885 | — 0,8 


| | 
— 0,6703 | 0,6065 | — 0,5904 





0,2482 02236 | 0,2088 
0,4500 08877 | 0,3613 
0,5678 0,4735 | 0,4383 
0,5950 | 04813 | 0,4411 





1,0497 | 1,1949 | 1,2419 
1,1345 1,2571 | 1,2972 
1,3204 1,4005 1,4267 
1,5251 1,5626 | 1,5740 





0,9497 0,7951 | 0,7424 
0,9295 0,7769 0,7243 
0,8700 0,7218 0,6711 
0,7737 0,6326 | 0,5843 
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Die Entscheidung darüber, 
welche Gebiete stabil und welche 
instabil sind, kann nach einem Satz 
von Haupt getroffen werden [6, S. 15, 
Satz I. Auf den Abbildungen sind 
die so ermittelten stabilen Gebiete, 
wie üblich, durch Schraffur kenntlich 
gemacht. 

Nach demselben Satz kann auch 
über Stabilität und Instabilität der 
Lösungen auf den Grenzkurven ent- 
schieden werden. Die Lösungen sind 
für diese Parameterwerte im allge- 
meinen instabil; nur zu den Para- 
meterwerten der Schnittpunkte zweier 
(gleichzahliger) Grenzkurven gehört 
eine stabile Lösung. So gehören zu 
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den Schnittpunkten der C,- und S,-Kurve (für positives y,) stabile Lösungen, ebenso zu den 
Schnittpunkten der C,,,- und $,,»- Kurve auf der A- Achse. 


Bild 5 zeigt noch den Schnitt y,=0 durch den Stabilitätskörper. Er stellt also die 
Stabilitätsverhältnisse der Differentialgleichung y"’—+(4+y,cos2x) y=0 dar. Die Kurven C,,, 
und S,,, liegen jetzt (zusammenfallend) ganz im (ersten) stabilen Bereich. 
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Der starre Kreiszylinder im isotropen elastischen Medium. 
Von F. Reutter in Karlsruhe. 


Aus dem elastischen Halbraum z>0 wird ein solcher Raumteil ausgespart, 
daß in ihn ein starrer Kreiszylinder eingepaßt werden kann. Die kleinen 
Verschiebungen, die die Punkte des elastischen Mediums durch das Ein 
drücken dieses Zylinders erfahren, werden bestimmt. Ihre Ermittlung 
führt auf zwei lineare Integralgleichungen erster Art, deren Lösung durch 
Zurückführung auf ein System von unendlichvielen linearen Gleichungen 
mit unendlichvielen Unbekannten gelingt. 


1. Einleitung und Formulierung des Problems. 


Zu den Problemen der räumlichen Elastizitätstheorie, die eine geschlossene Lösung ge- 
funden haben, gehört der starre Stempel auf einer ebenen elastischen Unterlage‘). Die vor- 
liegende Arbeit hat eine Erweiterung jener Untersuchung zum Ziel. Aus einem elastischen 
Medium — es sei der Halbraum 2 >0 gegeben — werde ein kreiszylindrischer Raumteil 
0=2=<k vom Radius c herausgenommen (Bild 1). 
In den ausgesparten Teil denken wir uns einen 
Zylinder vom Radius c so eingepaßt, daß der 
Zylindermantel am Rande des elastischen Mediums 
fest anliegt. Der Zylinder werde durch eine in 
Richtung der positiven z-Achse wirkende Kraft 
in das elastische Medium eingedrückt. Dabei soll 
der Zylinder als starr angesehen werden, d.h. sein 
Material sei gegenüber dem des elastischen Mediums 
so hart, daß seine Formänderung vernachlässigt 
werden kann. 














Unter dem Drucke des Zylinders erleiden alle an der Zylindergrundfläche anliegenden 
Punkte des elastischen Mediums eine konstante kleine Verschiebung in Richtung der positiven 
z-Achse. Wir nehmen ferner an, daß längs des Zylindermantels keine Verschiebungen in 
radialer Richtung erfolgen und es soll längs des Zylindermantels die Schubspannung gleich 
der Normalspannung multipliziert mit der Haftreibungsziffer sein. Aus diesen Annahmen 
ergeben sich sofort die Randbedingungen unseres Problems (Nr. 3). 


Von praktischer Bedeutung ist diese Untersuchung u. U. für die Baugrundforschung 
(z. B.: Einrammen eines Pfalles in den Baugrund. Das Material des Pfahles ist so hart zu 
denken, daß seine Formänderung gegenüber der des (elastischen) Baugrundes vernachlässigt 
werden kann). 


Die Verschiebungskomponenten genügen der Bipotentialgleichung. Deren Lösungen 
lassen sich mit Hilfe von Potentialfunktionen aufbauen. Man erhält so für die Verschiebungs- 
komponenten je von einem Parameter «a abhängige Partikulärlösungen, die aber die Rand- 
bedingungen noch nicht erfüllen. Durch Integration gewinnt man hieraus -Lösungen, die 
zwei unter dem Integralzeichen stehende Funktionen g(a) und h(a) enthalten. Diese sind 
gemäß den Randbedingungen zu bestimmen. Ihre Ermittlung führt auf lineare Integral- 
gleichungen erster Art. Während sich die Lösung h(a) aus bekannten Integralrelationen 
für Besselsche Funktionen sofort ablesen läßt, ergibt sich g(a) zunächst in der Form 

n 


g(a) = \ o(ß)Jır3(ac)a-?dß. Die „Verteilungsfunktion“ o (B) ist wiederum durch eine lineare 
1) 


Integralgleichung erster Art bestimmt. Indem man das uneigentliche Integral durch unend- 
liche Summen (bei halbzahligen /-Intervallen) ersetzt, erhält man g(a) als Summe zweier 
unendlicher Reihen, deren Glieder Produkte aus Besselfunktionen steigenden (ganz- bzw. 
halbzahligen) Parameters und negativen -Potenzen von a darstellen. Die Koeffizienten dieser 
Reihen lassen sich mittels eines Verfahrens, das für die Behandlung solcher linearer Integral- 
gleichungen erster Art von grundsätzlicher Bedeutung ist, durch ein unendliches Gleichungs- 
system bestimmen. Dessen Auflösung wird im ersten Schritt angegeben, und es wird hier- 
aus u. a. die Schubspannungsverteilung längs des Zylindermantels und längs des Deckels 
ermittelt und in Diagrammen dargestellt. 


1) E. Trefftz: Math. Grundlagen der Elastizitätstheorie in Frank-Mises, DifferentialgIn. d. Physik, Bd. II, 
Braunschweig 1925, S. 306 ff. 
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2. Die Grundgleichungen und ihre Partikulärlösungen. 


Zur analytischen Behandlung des in Nr. 1 formulierten Problems führen wir Zylinder- 
koordinaten z (Bild 1), r (Abstand von der z-Achse), ® (Winkel einer beliebigen Ebene durch 
die z-Achse gegen eine feste Ebene) ein. Da der Zylinder nur durch eine in Richtung der 
positiven z-Achse wirkende Kraft in das elastische Medium eingedrückt wird, sind die Ver- 
schiebungen vom Winkel d unabhängig. Daher ist der Verschiebungsvektor u für einen be- 
liebigen Punkt P:(r,z) des elastischen Mediums gegeben durch 


u=g .e+y:Pf. 
e ist ein Einheitsvektor in radialer Richtung, f der Einheitsvektor in der z-Richtung. gyundy 
genügen den Differentialgleichungen: 
1 1 u 
A4dyp- ziert Fe „419 rap 0, 14x=0 (4 Laplacescher Operator). 
Wir führen die biharmonischen Funktionen y und y auf harmonische zurück durch den 
Ansatz’) 


(1). 


o(r,z2)=n(r,2)+(2—k)y,(r, 2) \ 
] 


ı(r,2)= {(r,2) +(2— k)y;,(r,2 
Die Funktionen 7 (r, 2), £(r,2), y(r,z) genügen den Differentialgleichungen 


i i y —=( 
und es gilt auch 


Ay,=d, 
Außerdem liefern die Gleichgewichtsbedingungen der Elastizitätstheorie die Beziehungen: 


1 m 0 

L ) — iv — () 

ig „9 * 3, Av u . 
m 0 


- - dvu=0 
m —20z2 e 


Ax-+ 
m Querkontraktionszahl. Hieraus folgt unter Berücksichtigung der Gl. (2): 


Mm 


1 . 
HET — ı(- tr +2) 


Die beiden ersten Gln. (2) werden durch die Ansätze: 


n(r,2)=R(r) Z (2) S(r,2)=X(r) Y(2) 
gelöst. Dann gilt: 


1 | 1 
R'+.R4+(#—,)R=0 "+ X+BX=0 
zZ" — A'Z=0 Y'-BY=0. 


Da die Verschiebung für r=0 und 2>x nicht unendlich groß werden darf, erhält man zu- 
nächst die Partikulärlösungen 


n=aJ,(ar) e”*? PEN EEE Er 
a,b,a sind beliebig wählbare Größen. Daher wäre die allgemeinste Lösung: 
= \ ga) J,(ar)e"""da, {= \ h (a) J,(ar)e”*" da. 
0 0 
Die Funktionen g(a) und h(a) sind vorläufig unbekannt und später gemäß den Rand- 
bedingungen zu bestimmen. Wegen des Ansatzes (1) erhielte man jetzt für z=k: 


{ 2) Vgl. hierza Frank-Mises: Differentialgln. der Physik, Bd. I, Braunschweig 1930, S. 848, oder Grammel- 
Biezeno: Techn. Dynamik, Berlin 1939, S. 126. 





158 Reutter, Der starre Kreiszylinder im isotropen elastischen Medium 5; Bf; inth Moon. 





P)=|g(w) J,(ar)e”**da, 
N) 


x(k) )=|h(a )J,(ar)e-*kda. 
v 


Soll für die Punkte z=k, 0OS£r<c eine konstante Verschiebung in der z-Richtung statt- 
finden, so muß 

C(k)= M = konstans 
2Msin me. 
ZT a 


sein. Diese Bedingung ist erfüllbar für h(a)e-** = (s. Nr. 3, Gl. 7). Dann wird 


J,(ar)ek-2)da. 


SR JE! sinca 


ZT 
0 


Dies Integral ist aber für z2< k divergent und man gelangt nur dann .zu konvergenten ‚Inte- 
gralen, wenn man den Ansatz (4) der Partikulärlösungen von vornherein wie folgt abändert: 


n=aJ,(ar)e-®i:-kl, Z=bJ,(lar)e-ei:-kl, 


Damit besitzen zwar die Partikulärlösungen einen Sprung der ersten Ableitung an der 
Stelle z=k, doch läßt sich die Funktion g(a) so bestimmen, daß die durch Integration aus 
diesen Partikulärlösungen gewonnene endgültige (d. h. die Randbedingungen erfüllende) Lösung 
in den Ableitungen der Verschiebungskomponenten für 2=k keinen Sprung aufweist. Be- 
kanntlich „glättet“ die Integration. Mit dem neuen Ansatz: wird: 


n=|g(a) J, (ar) e-«!:-kIda, 
1} 


= nad, (ar)e-“*|2-klda. 
Aus Gleichung (3) folgt dann durch Integration naeh z: 
v2) [h (a) # g(a)] J,(ar)e-eie-kdatwir) . . »..:..60) 


mit Pe Dabei ist w(r) eine Potentialfunktion, die sich aber auf eine additive Kon- 
stante reduziert. Mit den aus (5) folgenden Werten für y,,y, erhält man schließlich: 


— (gla)J,(ar)e-e-kdatple—k ){ Ina) ) # g(a)]a J,(ar)e-«i:-kida 
1) 1 ö 


bo} 6) 


x= |'h(a) J,(ar)e-ei:-klda—p(z— k)|[g(a) F h(a)]a J,(ar)e-«!2-kida 
v v 


Das obere Vorzeichen gilt jeweils für z2>k, das untere für 2< k. 

Außer der durch den Produktansätz erhaltenen Partikulärlösung (4) für 4{=0 scheiden 
alle übrigen Partikulärlösungen dieser Gleichung von vornherein aus, da sie durchweg unzu- 
lässige Singularitäten besitzen, indem sie an bestimmten Stellen zu unendlichen Ver- 
schiebungen führen; 


3. Die Randbedingungen. 


Die Funktionen g(a) und kh(a) im Ansatz (6) sind jetzt gemäß den Randbedingungen 
zu bestimmen. 


a) Für 2=k, 0Sr<sc sollen alle Punkte der Deckfläche des Zylinders eine konstante 
kleine Verschiebung M in der z-Richtung erleiden. 


b) Für r=c, 0Sz2=<k soll keine Verschiebung in radialer Richtung erfolgen. 


e) Für r=c, 0<Sz<k soll die Schubspannung 1,, gleich der Normalspannung 0, 
multipliziert mit der Haftreibungsziffer f, sein. 
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Außerdem gehört für r >c die Ebene 2= 0 zur Berandung. Wir stellen für diese aber 
vorläufig noch keine Bedingung. Es wird sich später zeigen (Nr. 7), daß unsere im folgenden 
gewonnene Lösung für 2=0, r>c nahezu verschwindende Schub- und Normalspannungen r,,, 
o, liefert, was man eigentlich von vornherein als Randbedingung zu fordern hätte. 

Nachdem wir die Ebene 2=k, die die Deckfläche des Zylinders enthält, durch die 
Randbedingungen wie durch den Ansatz (6) ausgezeichnet haben, ist noch die weitere Be- 
dingung zu stellen: 

Für 2=k dürfen die Ableitungen der Verschiebungskomponenten keinen Sprung er- 
leiden. Denn diese sind lineare Funktionen der Spannungen. Wir schreiben also nur längs 
des Mantels zwei Bedingungen in der üblichen Form vor, längs des Deckels geben wir nur 
eine Verschiebungskomponente und fordern Stetigkeit der Spannungen in der Ebene des 


Deckels für r >c. Aus den Gln. (6) errechnet man aber leicht, daß sq .g von vornherein 


or’ da 
für 2=k stetig bleiben, daß aber =, 2 dann und nur dann keinen Sprung erfahren, 
wenn die Funktionen g (a) und h (a) die Bedingungen 


d) \ah(a)J,(ar)da=0 für r>e, 
N 


x 


e) \ag(a)J,(ar)da=0 für r>c 
0 


erfüllen. Die Randbedingung a) findet dagegen folgende mathematische Formulierung: 


a) |\h(a)J,(ar)da=M für O<r<ce. 
0 
Durch die Bedingungen a) und d) ist die Funktion h (a) vollständig bestimmt, man erhält®) 
RER... SUSE \ 


7 a 


4. Die Integralgleichungen für die Funktion 9 (ae). 


Wesentlich schwieriger ist die Ermittlung der Funktion g(a), die zugleich den Be- 
dingungen b), ce), e) genügen muß. Zunächst ist 


J,{ac)al-*” mit v»>1 


eine die Bedingung e) erfüllende Funktion, denn es gilt*): 


ao 


| J, (ar) J, (ac) a-"da=0 für r >c. 
0 


Das Integral ist konvergent für beliebige (rationale und irrationale) »>1. Also ist die all- 
gemeinste der Bedingung e) genügende Funktion 


se=lofhrsladatdB. ....:. 2.2.2.2... 
0 


mit o(0))—=0. Die Bestimmung von g(o) ist daher auf die Ermittlung einer Funktion o (ß) 
zurückgeführt. Wir ersetzen nun das uneigentliche Integral durch eine unendliche Reihe mit 
diskreten (halb- und ganzzahligen) $-Werten: 


a 


9 (a) = P2 a,J,,,(ac)a”” Bi: u, mlaa +, 


v—_ vw—1 


Der Bedingung o (0)=0 ist damit bereits Rechnung getragen. Alle Funktionen J,,,(ac)a”” 
und J,,,„(ac)a="**"* sind ungerade Funktionen, die an der Stelle a=0 eine Nullstelle 
erster Ordnung besitzen. Die Funktionen einer jeden der beiden Reihen 


3) Frank-Mises: Bd. I, S. 420, Bd. II, S. 3%. 
4) Nielsen: Handbuch d. Theorie der Zylinderfunktionen, Leipzig 1902, S. 199. 








. fi . f fi ‘ . 2. v. Mat ri 
160 Reutter, Der starre Kreiszylinder im isotropen elastischen Medium Ba.23 Nr.3- Tun de 





ı 


9. (a) = — a,Jy,,(ac)a”” 
—1ı 


(9) 


ao 
u; PR | 19 
9; la) = ji; b,J,+,2(ac)a 
s—| 
sind untereinander linear unabhängig. Zur (gleichmäßigen) Konvergenz dieser Reihen reicht 
es hin, wenn die Koeffizienten a,, b, beschränkt sind: 


la,I<N, b,I<N. 
ei J W-, 
Denn da N’ J,,,(ac)®) (gleichmäßig) konvergiert, konvergiert > 72,0 für |a| >1 
v—1i 
vl 


offensichtlich, für |a|< 1 aber wegen 


1 Jaf*: (5) 


hal< ri) ae 


_ ebenfalls und zwar gleichmäßig in a. Mit Hilfe dieser letzteren Beziehung zeigt man ebenso 


= 


leicht die Konvergenz von > Jr + 42 (a 4 Reihen des Typus (9) sind bisher wohl nicht 


vl 
bekannt. Bei Nielsen treten zwar Reihen mit einer ähnlichen Struktur auf, es läßt sich 
z. B. zeigen, daß jede ungerade Funktion in eine Reihe der Form®) 


a 


_ s—1—k 
fi) = > 4, (3) Jr+s(x) k fest, keine negative ganze Zahl 


s=1 
entwickelt werden kann. 


Wir suchen nun die Koeffizienten a,, b, im Ansatz (9) so zu bestimmen, daß die Rand- 
bedingungen b), ce) erfüllt sind. 


Die Bedingung b) verlangt 


b) | g(a)J, (ac) e-*!12—k da+p(2z—k) | [h(a) + g(a)] a J, (ac) e=-*2-klda=0 
f) ö 


für OSz<sk, 
entsprechend ce): 


0 d 





9,9%, 
Mn tZr 





für O<2<k, 
ge 


wo die Ableitung an der Stelle r=c aus den uln. (7) zu bilden ist. Die Gleichung c) enthält 
nach Ausrechnung die Funktion g(a) ebenfalls linear. Beide Gleichungen können auch als 


Bedingungen für die Funktion o(ß) (Gl. 8) geschrieben werden. Grundsätzlich ist dann 
jede der Randbedingungen vom Typus 


x 


| Jo(9-t(r,2,a)Jırs(ac)da dB=w(r,2). 
06 


o(ß) ist vollständig bestimmt, wenn für jeden Teil der Berandung (Bild 1) eine derartige 
Bedingungsgleichung vorgeschrieben wird. 

Wir geben nun den beiden Bedingungen b), c) noch eine Form, die zeigt, daß sie 
eigentlich nur eine Integralgleichung 1. Art für g(a) darstellen. 

Die Randbedingungen b), c) sind ja zunächst zwei Integralgleichungen 1. Art für ein 
und dieselbe Funktion g(a) von der Form 


B) 


|K,(a,2) g(a)da= 6, (e) 
3 O<z<k. 


= 


| K,(a,2) g (a) da= G,(e) 


v 


6) Frank-Mises: I, S. 447. — Graf-Gubbler: Besselfunktionen, H. 2, Bern 1900, S. 71. 
®6) Nielsen: Handbuch, S. 267. — Nyt Tiddskrift för Mathematik 9 B (1898), S. 83. 
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Das Grundgebiet ist aber 0O< 2< ©; über das Intervall k < 2< © des Grundgebietes ist gar 

9 nicht, über das Intervall O< 2< k doppelt verfügt. Nach Ausführung der Substitution 2’ = 2 
r k Jr ’ f k ’ ‚ r» . . 
mit R, (a, z)= (a, 2’), @,(2)= (2)=6 (2’) und Wiedereinführung der Bezeichnung 2 
an Stelle von 2’ genügt also g(a) den Gleichungen: 
ı reicht 
| K,(a,2)g()da= 6, (e) 0<z<k 
ö 
a|>1ı | Ki (a, 2) g(a) da= 6} (e) k<e<o. 
ö 
z _J$K, (a, 2) 0<2<k 
Kan=[0n k<2<o 
Oder endlich mit 0.) 
2 0<z<k 
“= ‚ k 

ebenso G, (2) k<2<o 

gilt die einzige Integralgleichung 1. Art mit dem für z=k unstetigen Kern K(a,z) und der 
nicht für 2=k unstetigen Störungsfunktion @ (2): 
st sich | K(a,2)g(a)da=G (e) 0O<2<w 
v 

Damit ist auch die Verträglichkeit der beiden Gleichungen b), c) für g(a) nachgewiesen. Auf 
die weitere Frage, ob die Bedingungen b), ec) zur eindeutigen Bestimmung von g(a) aus- 
reichen, kommen wir in Nr. 7 zurück. Entsprechendes gilt auch für die Integralgleichungen 
für o(ß). 

Rand. Unter Mitberücksichtigung der Sprungbedingung e) kann endlich der. Kern auch als 
Funktion von 3 Veränderlichen’), die Störungsfunktion als Funktion von 2 Veränderlichen 
geschrieben werden; es gilt die Integralgleichung 

MT 777 Vene 
N 
mit 
G, (2) für z<k,r=c 
G(r,2 =! 0 „ mh, r>e 
G/ (2) „sw>rhr=e 
. K,(a,2) fürz<k,r=c 
-_ K(r,2,a)={faJ,(ac) . serhr>e 
dann IK, (a, z) .„ s>h, r=e. 
Damit sind die Bedingungen b), c), e) zu einer einzigen Integralgleichung zusammengefaßt. 
5. Ersatz der Integralgleichungen durch ein unendliches Gleichungssystem. 
rtige Wir behandeln die Integralgleichungen b), c) mittels einer Modifikation eines bekannten, 
: aber praktisch wohl noch wenig angewandten Verfahrens®). 
sie 
Ist die Integralgleichung 1. Art 
ein l 
S K(x,2)f(e)de=g(x) 
vorgelegt, sind aber die Eigenfunktionen des Kerns K (x, 2) nicht bekannt, so bilde man aus 
der Integralgleichung durch Multiplikation mit den oo' Funktionen y, (x) eines abgeschlossenen 
Funktionensystems, das wir gleich als orthogonal annehmen wollen, und Integration nach x 











die ©o' Gleichungen 


”) Hellinger -Toeplitz: Integralgleichungen und Gleichungen mit unendlichvielen Unbekannten in. Enzy- 
klopädie der mathematischen Wissenschaften, Bd. II, III,, S. 1389. 

8) S. z.B. Hamel: Integralgleichungen, Berlin 1937, S. 93/94. 
11 
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2 2 
\\Kad)f@)y(@)dzede=\g@)y,(de=(, v=1,2.... 
00 1) 

l 
Oder mit \K(&,2)y,(«)de=u, (2): 
v 
I 
\@w(@)d2=(C, v=1,2,... 
vu 
Bilden die «,(z) bereits ein Orthogonalsystem, so ist 
fd)= PP} C„u,(2), 
v—1 

im anderen Falle sind die «,(z) erst zu orthogonalisieren. 

Im vorliegenden Fall läßt sich das Verfahren nicht unmittelbar anwenden, da f(2) ja 
teilweise schon durch die Bedingung e) festgelegt ist. Wir können aber jetzt f(z) als Reihe 
nach den ©" Funktionen y.u(2) — s. Ansatz (9) — ansetzen und erhalten: 

fi2) = Ya Xu (2). 
fu 

Mit 

l 
| Xu (2) Uy (2) dz= Bu v 
0 
ergeben sich daher die ©' Gleichungen 
2 Bu v Au = C, 

zur Bestimmüng der Koeffizienten Au. 

In dieser Weise behandeln wir die Integralgleichungen b) und c): Diese Gleichungen 
enthalten noch die Veränderliche 2. Es werden jetzt beide Gleichungen mit cos m (eh) 
multipliziert und über z von 0 bis k integriert, zugleich wird der Ansatz (9) in b) und ec) 
eingeführt. (Es würde zunächst naheliegen, an Stelle.des Funktionensystems cos (m - (— ») 
ein System von 0! Zylinderfunktionen zu verwenden, doch erweist sich dies als nicht durch- 
führbar.) 

Nun spaltet jede der Gleichungen b) und c) in je unendlichviele Gleichungen 
(n=0,1,2,...) auf und man erhält >"! Gleichungen für die ©o' unbekannten Größen a,, b,.: 

Setzt man 

0 
n a 
Cm (a) = \ e-“lk—2| cos (m 5 —r))d (k — 2) = Tr ((— 1)" e-k — 1), ’ 
e 7T 
k a@+ ee 
® n k n k k 
BE. —ak__ ch —_.a@a 
dm = | Herne 9-7 2m’ ei ( un. ı I Arc bene AR vol a ld A Be R 
®n n?n? 
(+) + | 


so gelangt man durch das angegebene Verfahren schließlich zu folgendem System von 2X &' 


Gleichungen: 
y nv nv a? p’ nv 
I ala a -pD+ola kp at a 

v—1i 


hd nv nv p?’ wi, Bi 11) 
+ > 2, lo p)+atal-p)—kp rt 2 on + ai) 


v1 








2 2 
= pt + KH) + en +) n=0,1,2,. 





. Mech, 
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(2) ja 
teihe 


Igen 


-h) 


l ce) 


h) 


rch- 


gen 
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ir nv v 4 n® v 1 
) 4, I 2f(m —1—p) (op + ®ore) u (m — 1) fs p loss + woı 
v‚—1 
En 2 (m — 1) A pP (oo E= k oo3 u wo) - ( 1 E= p) Im — 2) (wi + 112) 
+ Ai fü 9} nv ER nv 2 fo Oo nv k nv nv 
Er p (m— 2) (w134+ 185) — . plm 2) \wısı + k wiss + 0132) 


+6 — 


(m —2)(wj3ı + @132) + (fe +m—2)(oj4+ wi) 


-2p (+ m Yloli+kotstom)) 


ar nv nv 4 n? f nv nv 
+ > b, -2f, (m —1 p) (oo2ı + @u2) — pr (m 1) fh p log +w y 


v‚=1 
+2 (m —1)f,p (oysı + k on + wu) + (L+p) (m — : 2) (0151 + @Io) 


f "LP 


4 ” o nv nv nv nv nv 
na gr. pm-— 2) (ja + @is) — 7 m — 2) (0141 + k wis + @ia2) 


+ . (m — 2) (@i4ı + 142) + u. nr (fs — m — 2) (734 4 0185) 
-2p ll +m-Ylolı +kola+ ob) 


[_,n n\ 4 
—=2f,(1—- pP) +%0o) + c P pfa(lm— 2) (08 + x04) 


—2pfolm 2 (al, + ka + u) (pP — lm — 2) (xı + #35) 
Bi. n? > f Ä 
\ 


pe 4 
zu m 9° 


x»14+ x) +2p (m 232 Le (et + kai + #10) 


2 


47 
u (m 2) pl + u) +2plm May +krnt x %) 


=01,2,. 


Die Größen w},., %,, sind Integrale von der Form 


A| 





025) 





ao 6 
le a Tl eye 2 et 
v 1) 
= [+] 
all = | ah.” ug ae erde 
a + k? \ a +- k? 
v0 V 
on ao 
Pr SE J; (ac) ru SR BR LU LU ERTT I O Do 
|" a + 5 A a +- pe 
0 0 
© © BSR 
a, = Ber Si 14 = lade „or da 
| +) | (+) 
0 0 














(12). 


(13) 
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Br RL 4 Be. a-e-kda et J,laejsinae as da 
® - ——— ” a’+ nn“ 
a [5 p 
0 v 
(14) 
2M J, (ac) sinac ’ 2 2J (ac) sin ac 
“= 7 PrELH e-kda = | ge e-“k da 
+) Ne+) 
0 v 





usf.; sie sind z. T. noch mit: Vorzeichenfaktoren (— 1)" zu versehen, doch wird hier von ihrer 
ausführlichen Aufzählung abgesehen. 


6. Auflösung des unendlichen Gleichungssystems. 


Wır behandeln das unendliche Gleichungssystem (11), (12), indem wir es sukzessive 
durch ein endliches System mit einer immer größeren Anzahl von Gleichungen und Unbe- 
kannten ersetzen. Dazu denken wir uns die Gln. (11) in der abgekürzten Form 


= 0} 

% 

z ara, + Bd, cn, „=1,2,... 
vw—] v—1 


die GIn. (12) in der Form 


ao 


z Yv»a,+ B; ö.,b,=d,, #»=1,2... 


v—L ul 


geschrieben. Von den =" Lösungen a,, b, ist zu verlangen, daß die mit ihnen gebildeten 
Reihen .(9) konvergieren. Dazu genügt es nach Nr. 4, daß die a,,b, beschränkt sind: 


la,I<N, Ib,|<N. 


Unter dieser Annahme ist eine Theorie der unendlichen linearen Gleichungssysteme von 
H. v. Koch?) entwickelt worden, die sich des Apparates der unendlichen Determinanten 
bedient. Für unsere Zwecke ist aber trotz der andersgearteten Anforderungen an die 
Lösungen der Anschluß an die Theorie von E. Schmidt‘) vorteilhafter, deren Konvergenz- 
bedingungen hier leichter zu überblieken sind. Notwendige Bedingung für die Existenz 
der Lösungen x, eines Gleichungssystems von der Form 


2 Auv, = Cu 
unter der Voraussetzung, daß )' x’ konvergiert, ist nach E. Schmidt die Konvergenz von 
+ AR, für jedes u sowie die Konvergenz von 2» C}. Diese Bedingungen sind für unser 
Gleichungssystem (11), (12) erfüllt. Denn in der Zeile (variabler Index », feste Nr. n) stehen 
von links nach rechts, wie man aus den Gln. (13) erkennt, Besselfunktionen steigender Ord- 
nung unter dem Integralzeichen. Die Absolutbeträge der Integrale lassen sich aber ins- 
gesamt unter Verwendung von nur sechs verschiedenen Integralen majorisieren. Es ist z. B. 


a 


| - )J, i 2? -v k? » 
ol = \* (ac) Fr u e-kdel Err) . ana 


ge m 1 
u 7 167) er —y! 


wo II(z) die Gaußsche I/-Funktion darstellt. Sämtliche übrigen Integrale lassen sich in 
derselben Weise mit Hilfe der fünf folgenden abschätzen: 


%) H. v. Koch: Acta mathematica, Bd. 16 (1892), S. 217 bis 295. 
10) E. Schmidt: Rend. di Palermo, Bd. 25 (1908), S. 53 bis 77. 
11) N Sonine: Math. Annalen, Bd. 16 (1880), S. 39. 
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, u. da= 
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. n(® 5 =) 


n(5) | 9, (5) 


1-v 
I4.. (ac) a'-*’ da= ;- 


, n(® = 1 


L) 


- 1 . 9512 + 
932» (5) N | >; ’da— er ae 2 
J,,,,(ac)a’’"da= Bs2-» 1 , r 1 = h 
ce? 1Iw—1) > 
0 0 
© Ze 
\ 972-v II (5) 
Jr u (ac) a" "dam u 


ü ce? 1IW —2)' 
v0 
. . . . > Cu . r » » 
Alle wei nehmen daher ab wie eine Reihe mit Gliedern wei Die Koeffizienten an», An» der 
Gleichungensysteme’ (11), (12) sind aber lineare Kombinationen der Größen w,,., somit ist 
die Konvergenz von 
a} y " 2 Y 2 ry 2 
P: Ay» + Sr Pi v und 2 Yn ‚+ = 0 
für jedes n nachgewiesen. Es konvergiert sogar schon die Reihe der absoluten Beträge 


P3 Anv + 3» Bnv und 2 Ynv +) Önr . 


zz 


Noch einfacher läßt sich die Konvergenz der Quadratsumme der rechten Seiten der 


1 
r abnehmen. 
l+n? 
Deren Konvergenz ist nach E. Schmidt'?) zur Existenz eines Lösungssystems x, mit kon- 
4 R . Be . 
vergenter I’x, auch hinreichend, wenn nach ÖOrthogonalisierung des unendlichen 


u 


Gleichungssysteme nachweisen, da diese nach den GlIn. (14) mindestens wie 


Gleichungssystems die Reihe der Quadratsummen der neuen rechten Seiten }) 9% ebenfalls 
konvergiert. Die g„ sind lineare Kombinationen der c, von der Form 


n 


u rn 
In= — YnvCv; 


v—1 


die r„„ durch den Orthogonalisierungsprozeß bestimmbare Konstante. Bei der komplizierten 
Form der Koeffizienten unseres unendlichen Gleichungssystems ist diese Orthogonalisierung 
nicht olıne weiteres durchführbar, die hinreichende Bedingung also zunächst nicht nachzu- 
prüfen. Wir dürfen aber wohl um so eher annehmen, daß sie erfüllt ist, nachdem nicht nur 
die Quadratsummen der Koeffizienten, sondern schon die Reihen ihrer absoluten Beträge 
konvergieren. Zudem sind bei E. Schmidt alle Bedingungen unter der Forderung der Kon- 
vergenz der Quadratsumme der Lösungen aufgestellt, während es hier sogar ausreicht, wenn 
die Lösungen beschränkt sind. 

Um endlich die Eindeutigkeit der Lösungen zu gewährleisten, müßte nach E. Schmidt 
die Quadratsumme der Spaltenquadrate der linken Seiten den Wert 1 besitzen‘’). Dann ist 
das vorgelegte Gleichungssystem „vollständig“. Wir können nun zwar leicht die Konvergenz 
der Spaltenquadrate 


Y 2 


> a}, „+ pP Yn v und pP: Bar + pP ö2 v 
nachweisen, da innerhalb einer jeden Spalte nach den Gln. (13) die Glieder mindestens wie 
1 \ ’ u. A ’ 
ir" abnehmen. Der Nachweis, daß diese Summen jeweils den Wert 1 besitzen, gelingt 


aber natürlich nicht ohne weiteres. Amdererseits ist die Existenz einer eindeutigen Lösung 
des Randwertproblems durch die Existenzsätze der Theorie der Randwertprobleme ge- 
währleistet. D.h. die Eindeutigkeit der Lösung folgt daraus, daß unsere Reihen g, (a), 9, (a) 
nur die Funktion o (ß) ersetzen, die wiederum durch die Randbedingungen eindeutig bestimmt 
ist, allerdings unter der Voraussetzung, daß auf allen Rändern je zwei Spannungs- bzw. 
Verschiebungskomponenten vorgeschrieben sind. Unser Ansatz genügt aber dieser Forderung 
formal zunächst nur für die Mantelfläche des Zylinders. Wir lassen nun die Frage nach der 
Eindeutigkeit noch offen und kommen darauf in Nr. 7 zurück. 





12) E. Schmidt: ]. e., S. 73. 
13) E. Schmidt: ]. c., S. 69. 
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Der erste Schritt der numerischen Auflösung wurde durchgeführt und dabei die 
Koeffizienten a,, b, bestimmt. . Die auftretenden Integrale &},,, #,, wurden durchweg nach 
der Simpsonschen Regel auf 5 Stellen berechnet'*). Die Haftreibungsziffer f, wurde zu 
0,5 angenommen, die Querkontraktionszahl m = 4 gesetzt. Ferner wurde angenommen 
k=c=1. Zur praktischen Rechnung wird in den durch die GIn. (13), (14) definierten Inte- 
gralen die Substitution ae=w ausgeführt, so daß insbesondere 
ge ae ei.) 7 A)? u Tan 
u er zur“ k? 

c o 6) 9 
—u’—+n?’n? 

j e 
1] 


’ a 
usf. Dadurch wird erreicht, daß nur noch Abhängigkeit vom Verhältnis „ besteht und daß 


die auftretenden Besselfunktionen von k und c ganz unabhängig werden. Daher kann die 
Rechnung auch für verschiedene Werte k, e rasch durchgeführt werden, wenn die Integrale 
einmal berechnet sind. 


Es ergaben sich so die Werte '*) 


2M 2M 
a, = — 1,30535:”—, b,= — 0,2694 -——. 
IT ZT 
Die numerische Auflösung bis zu weiteren Koeffizienten und die eingehende mechanische 
Diskussion soll einer späteren Mitteilung vorbehalten bleiben. 


7. Die freie O'b.itiache z=0. — Eindeutigkeit der Lösung. 


Die freie Oberfläche des elastischen Mediums, z=0, müßte an und für sich von 
Normal- wie Schubspannungen frei sein. Wir haben das bei Aufstellung der Randbedingungen 
nicht gefordert und können daher nicht erwarten, daß unsere Lösung dies liefert. Immerhin 

erhält man mit dem bisher numerisch durchge- 
Tu führten ersten Schritt die in Bild 2 wiedergegebene 
” Verteilung der Schubspannungen längs z=0,r >e. 


i i N M 
Als Ordinaten sind die Werte r,.: 2 aufgetragen. 
- Tr 


1 Dabei ist zu beachten, daß die Einsenkung M eine 
kleine Größe ist. In einer Entfernung von 5 Zylinder- 
radien ist r,, praktisch abgeklungen, nimmt aber 














\ M 
ER En ” für 2 Zylinderradien im Maße 7g och endliche 
Bild 2. Werte an. Die Verteilung der Normalspannungen 


o, wurde zwar noch nicht numerisch errechnet, 
doch sind die im Ansatz für o, auftretenden Integrale vom selben Typus wie bei r,,, so 
daß o, in ähnlicher Weise abklingt wie r,,. Wir können daher sagen, daß unsere Lösung 
zwar in der Umgebung des Deckels, d. h. in genügender Entfernung von der freien 
Oberfläche, nahezu streng gilt. An der Stelle, da der Zylinder in das elastische Medium 
„eintaucht“, werden aber die wirklichen Verhältnisse u. U. noch nicht genügend genau 
wiedergegeben. Eine exakte Berücksichtigung der Randbedingungen für die freie Ober- 
fläche ist aber in folgender Weise möglich: 


Das Verschwinden von (r,-),., kann von vornherein gefordert werden. Man gelangt 
dann zu der Bedingung: 


\ 


Y 
2ky.,,.—y, Hm +24 2)-24=0, 


und zwar kann man zeigen, daß diese Gleichung nicht nur für 2=0, sondern für alle 
Werte r,z gilt. Setzt man für y, seinen Wert aus Gl. (3) ein, so folgt nach zweimaliger 
Integration nach 2 


14) In den Werken von Watson (A treatise of the theory of Besselfunetions, Cambridge 1922) und Nielsen 
(Handbuch, s. Anm. 3) sind nur einige wenige dieser Integrale enthalten, so daß die numerische Rechnung unum- 
gänglich war. 

15) Für die Durchführung dieser recht umfangreichen numerischen Rechnungen bin ich Herru eand. math. 
Robert Claß zu Dauk verpflichtet. 








Z.a 


mi 


we 


fol 


0 














































Bar" a nt 1985 Reutter, Der starre Kreiszylinder im isotropen elastischer Medium 167 
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J,tar) 
nach yiz) = \ [A-g(a) # B-h(a)] ge BET te 
> zu : Ö 
men mit o a 
Inte- fr SERRERE m — 1 Dan m. 
k3m—4’ k(öm-—4) 
Die Bedingung {r,.),.„0 läßt sich also verwirklichen, ohne daß der Funktion g(a) eine 
weitere Bedingung auferlegt wird. Dagegen führt die Forderung (o,),.„=0 für r>e auf 
folgende Bedingung für g (a): 
\ a ga) J,(ar)e-"*da= \sincaJ,(ar) e-kda Beeren ., 2, Bm 
daß Re I Pe; 
Oder als Gleichung für o(ß) geschrieben: 
die 
nn 0 nm 
rale i RE. ; ge 
I \ J\o(MJIı+slac)J,(ar)-ae="*kdadß=\sineaJ,(ar)e-"*da, r>ec. 
0 0 0 
Es bleibt zu prüfen, ob diese Bedingung noch mit den beiden (bzw. drei) anderen Bedingungen 
für g(a) bzw. o(ß) — Nr. 4, s. auch Gl. (10) — verträglich ist. Doch kann dies geschlossen 
werden, da damit erst über alle Ränder verfügt ist. An Stelle der halbzahligen f-Intervalle 
f wählt man jetzt zweckmäßig #-Intervalle von der Größe w Die Funktion h (a) ist aber nach 
‚che \ i i , > i en R - 
wie vor schon durch die beiden Bedingungen a), d) bestimmt. Bezüglich der Eindeutigkeit 
der in Nr. 2 bis 6 gewonnenen Lösung ist daher zu sagen: 
Die Lösung ist die einzige, die den in Nr. 3 formulierten Bedingungen a) bis e) genügt 
und für r,, (und o,) die oben in Nr. 7 angegebenen Verteilungen längs 2=0, r > c liefert. 
Wir hatten — wie das auch manchmal bei der Behandlung ebener Randwertprobleme (mittels 
von konformer Abbildung) geschieht —, die Bedingungen für 2=0 noch freigelassen und können 
gen jetzt nachträglich die mit der gewonnenen Lösung sich einstellenden Randwerte (1,.).=o, 
hin (0,)... als vorgeschrieben denken. Da beide Werte nur in unmittelbarer Umgebung von 
1ge- r=c,z=0( merklich von Null verschieden sind, wird wegen der stetigen Abhängigkeit von 
ene den Randbedingungen unsere Lösung nicht sehr von derjenigen abweichen, für die (t,.).-o 
> c. und (0,),-. exakt Null sind. Diese exakte Lösung könnten wir aber unter Berücksichtigung 
ei der Gln. (15), (16) ebenfalls mittels unseres Verfahrens gewinnen. Wir werden hierauf 


zusammen mit der weiteren numerischen Auswertung in einer späteren Mitteilung zurück- 
ine kommen. 














ler- 8. Weitere Diskussion. 
ber Bild 3 zeigt die mit unserer nullten Näherung errechnete Schubspannungsverteilung — 
3 . Tr z ah . . ; . 
che d. h. wiederum 2: M — längs des Zylindermantels, die wegen der Randbedingung ce) zugleich 
J T 
zen — bis auf den Faktor f, — die Normalspannungsverteilung wiedergibt. Das Bild zeigt, daß 
1et, 
so i 
en en 023 050 a1 M > 
Ing | 
ien 1 00 05 WM_, | | 
um | T | | | | 
au 7} 7 | 
‚er- 2 
) 
ar 7 
gt . l 
Bild 3. | 
"frau 


Bild 4 (rechts). 


1 die Normalspannung durchweg von einerlei Vorzeichen, ein Ablösen des elastischen Mediums 
or vom Zylinderrand also nicht zu erwarten und daher die Randbedingung ce) wirklich sinnvoll 


zer R 

ist. Bild 4 endlich zeigt die Verteilung der Schubspannung, d.h. von a M am Deckel des 
ur Zylinders. Für r=c ist r,, unendlich groß. (Der .Meridianschnitt des elastischen Mediums 
ım- 





hat an dieser Stelle den Innenwinkel ed s. Bild 1) Dabei wird r,, aber von niedrigerer 


als erster Ordnung unendlich; denn mit 0 (ß) geschrieben ist 
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= \o(M)J,(ar)J,; ;(ac)a'-?dadß. 


00 


Nun gilt für jedes feste $ und c=1: 


- BER 
\ J (ar) J, ,s(ac)a'-? da=777 U >: 
ö 


Ein Pol liegt also vor für O<ß<1l. Er tritt auch auf, wenn die Forderung (r,,.);.„=0 — 
Nr. 7, Gl. (15) — gestellt wird; in diesem Falle erhält r,, ein Zusatzglied, das für r >e 
’ RER 2 \ 3 
singulär wird wie die hypergeometrische Funktion FS.5,2+B: 
Von Interesse ist ferner der Widerstand W, den der Zylinder beim Eindrücken in 
das Medium in Richtung der 2-Achse erfährt. Er ist gegeben durch den Ausdruck 


e k 
W=2anlo,rdr+2n| 1,,dz, 
0 0 


da längs des Deckels nur die Normalspannung, längs des Mantels nur die Schubspannung 
einen Beitrag zu W liefert. Unter Benützung der für o, und r,, sich ergebenden Werte — 
es ist 


0p7 


4 1 DL 
= ale + 2), =26% ss; 1 +52) 


ds "Ir . 'a-30r 'v’3 


r 


erhält man schließlich für ud : 


no 


9() J,(ac)da+Bp+1) | (e-“* —1)g(a) J,(ac)da 


„’ 


0 


G el-plm-2) 


W \ 
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x 


r 


+2pk \ ae-"k g(a)J,(ac)da 
ö 


+ sm pn +°R A| sin ca J,(ac)da 





k wo 
\ sincaJ, (ac) e-*k da! 
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Die Größe M (die Strecke, um die der Zylinder in das elastische Medium eingedrückt wird), 
ist auch in allen Gliedern des ersten Teils der rechten Seite von Gl. (17) als linearer Faktor 
enthalten. Denn nach Ersatz der Funktion g(a) durch die darstellenden Reihen (9) haben 
alle Koeffizienten a,,b, den Faktor M.. Die rechten Seiten des Gleichungssystems (11), (12) 


zur Bestimmung der a,,b, sind nämlich durchweg Integrale der Form % — Gl. (14) —, die 
alle den Faktor M enthalten. Daher kann jetzt mit Hilfe der Gl. (15) M durch den Ge- 
samtdruck W des Stempels ausgedrückt werden: 


M= A nd W, 
wo A ein durch die bestimmten Integrale und konstanten Faktoren festgelegter Ausdruck ist. 
Endlich erhält man für die Radialverschiebung o längs des Deckels: 
e= Prem = | gl) J,(arnda. 


Man zeigt leicht, daß o negativ wird, das elastische Medium wird also unter dem Stempel 
nach innen zusammengedrückt. 
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9. Die allgemeine Rotationsfläche. 
Das eingeschlagene Verfahren läßt sich grundsätzlich auch für eine beliebige Rotations- 
fläche verwenden. Man geht wiederum aus von den unter Nr. 2 aufgestellten Partikulär- 
lösungen und baut jetzt die allgemeine Lösung hieraus etwa in der Form: 


n n 


= | 9(a)-J,(ar)e=-*da— pz\|[g(a)—h(a)]aJ,(ar)e-*:da, 
VÜ 0 


x=\hla) J,(ar)e-da— pz\[p(a)—h(a)]a J,(ar)e-":da 
0 0 
auf. Für die Berandung 2 = f(r) der Rotationsfläche und längs der Oberfläche z2=0 des 
elastischen Mediums sind irgendwelche Randbedingungen vorzuschreiben, durch die die 
Funktionen g(a) und h(a) wiederum mittels linearer Integralgleichungen erster Art fest- 
gelegt werden. Wenn nur eine Berandung 2=f(r) mit geschlossener Darstellung vorliegt, 
erscheinen die Randbedingungen z. B. in folgender Form: 


PM = A == B 


(Hr); =0 (9,)2- = 0 


.. (18). 


Für g(a) und h(a) sind jetzt möglichst Reihen nach den Funktionen y„(a) eines solchen ab- 
geschlossenen Funktionensystems anzusetzen, daß zwei von den vier Bedingungen (18) schon 
erfüllt werden. Die beiden anderen Gleichungen sind dann in geeigneter Weise mit den 
Funktionen y,„(r) eines vollständigen Orthogonalsystems zu multiplizieren, es ist nach r zu 
integrieren und man erhält auf diese Weise 2X 0' Gleichungen für die 2X oo! unbekannten 
Koeffizienten der Reihen für g(«a) und h(a). Ohne Berücksichtigung der Bedingungen für 
z=Ü ist es jetzt u. U. möglich, für x„(a) und y,„(r) die Funktionen ein und desselben 
Systems zu wählen und damit mittels des in Nr. 5 geschilderten Verfahrens — s. Anm. 9 — 
unmittelbar die Fourierkoeffizienten von g(a),h(a) in bezug auf die y„(a) zu bestimmen, ohne 
auf das unendliche Gleichungssystem zurückgreifen zu müssen. 

Den Herren Prof. Dr. Haack und Prof. Dr. Sonntag danke ich für das der Arbeit 
entgegengebrachte Interesse. 450 


Zur Stabilität der elastischen Schalen II” 


Von Ernst-August Deuker in Hannover. 


VII. Die Stabilität einer geschlossenen Kugelschale unter gleichförmigem Außendruck. 


Die Stabilität einer geschlossenen Kugelschale, die durch einen gleichmäßigen äußeren 
Normaldruck beansprucht wird, hat zuerst R. Zoelly untersucht‘). R. Zoelly sowohl wie 
auch E. Schwerin'*) beschränken sich jedoch auf rotationssymmetrische Knickfiguren. Von 
dieser Beschränkung hat sich erst A. van der Neut freigemacht'*). Die in dieser Arbeit 
entwickelten Methoden zeigen auf, daß dieses Stabilitätsproblem einige hervorragende Eigen- 
schaften besitzt, die es zu einer der einfachsten Aufgaben der Stabilitätstheorie dünner 
Schalen machen. 

Es lassen sich drei Gründe angeben, die in ihrer Gesamtheit die analytische Behand- 
lung der Aufgabe wesentlich erleichtern: 

1. Eine besondere Eigenschaft der Kugel. 

2. Die Symmetrie der Grundlösung, die durch i. bedingt ist. 

3. Das Fehlen von Randbedingungen, die durch die Forderung zu ersetzen sind, daß 
die Lösungen auf der ganzen Kugel regulär bleiben. 

Als körperfeste Koordinaten auf der unverformten Mittelfläche mit dem Radius a wählen 
wir die Poldistanz q, und die Länge q,. Das Koordinatensystem ist rechtwinklig, und die 
Komponenten des Maß- bzw. Krümmungstensors haben die Werte: 


0 (n) 0) - 
m= a; =0 ge= sing, 
2 8.1). 
(0) __ B (0) _n. m ie ( ) 
z1=-4;, m=0; %9=—a sin’g, 


Die besondere Eigenschaft der Kugel besteht darin, daß ihr Krümmungstensor dem Ma&- 
tensor proportional ist. Es gelten nämlich die Beziehungen: 


*) Die Abschnitte I bis VII sind in dieser Zeitschrift Bd. 23 (1943), S. 81 bis 100 erschienen. 
14) Vgl.{11]. :5) Vgl. [9]. 1°) Vgl. [7]. 
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(0), 1 


(0) (0), 3 u % 
Kon — 000 = —-—h A 


er a 


Unter dem Einfluß des konstanten Außendruckes p [kg em” ?] geht die Mittelfläche in eine Kugel 

von etwas kleinerem Radius über. Für den Grundzustand bleiben deshalb die Beziehungen 

(8.2) erhalten, nur daß statt des Radius a der etwas kleinere Radius zu nehmen ist. Die 

Änderung des ursprünglichen Radius ist jedoch so gering, daß sie in den Gleichgewichts- 

bedingungen (3.7) bei der Berechnung der Grundspannungen unberücksichtigt bleiben kann. 
Die Gleichgewichtsbedingungen erhalten dann die Gestalt: 

(od Ti 

d g“ ; 

VOM 1 \2 


h 
9 u u Wr BE 
er +2 sm pl! +3,) =t 


—=o0 


Sie haben die spezielle Lösung: 


(0) (0) 


ne —- FE u BE... ern ii 


"Zur Bestimmung von T und M reichen die Gleichgewichtsbedingungen jedoch nicht aus, 
es ist notwendig, den Zusammenhang von T und M mit der Normalverschiebung zu be- 


2 


ae ö 
trachten. Man findet dann, daß sich M und T verhält, wie 1 Es 1. M ist also sehr klein 


gegenüber T. Man sieht hier, daß das Nullsetzen der Momente bei der Membranlösung mit 
der Verformung zwar nicht im Einklang steht, der Fehler aber unbedeutend ist, so daß wir 
M gleich Null setzen können. Für die Konstante T erhält man dann den Wert: 


r-"2(i4z.) te ee er VD 


2a 


wo, wie beim Kreiszylinder, der Klammerfaktor noch gleich 1 gesetzt werden kann. 

Die Gl. (8.4) zeigen die große Symmetrie der Grundlösung. Wie auf der Kugelfläche 
die Krümmung in einem Punkt für jede Schnittrichtung gleich ist, so auch im Falle kon- 

0 0) 

stanter Normalbelastung die Spannung. Die Tatsache, daß die Komponenten er und 7° in 
jedem Koordinatensystem auf der Kugel die gleichen Werte annehmen, gestattet es, die 
folgenden Rechnungen ohne Benutzung eines speziellen Koordinatensystems durchzuführen, 
und-ist der innere Grund dafür, daß bei der Betrachtung nichtrotationssymmetrischer Knick- 
figuren keine Änderung der Knicklasten gegenüber dem rotationssymmetrischen Fall auftritt, 
was zuerst Verwunderung erregt hat, durch die Invarianz der Tensorgleichungen (8.2) und 
(8.4) aber unmittelbar erklärt wird. Wir behandeln das Stabilitätsproblem zunächst nach 
Ziffer V an Hand der Gleichgewichtsbedingungen. 

Unter ‘Berücksichtigung von (8.2) erhält man aus (6.7) und (6.11) für die Variationen 
der Komponenten des bzw. Maß- und Krümmungstensors der verformten Mittelfläche die 
nachstehenden Ausdrücke, wenn man statt der Normalverschiebung öv, noch die dimensions- 


6%, . 
lose Größe döw— — einführt. 


Mo , Mo 
de Tip 
oM d2 öw 1 7 

“= ag) ieta ö, dw . 
Da in den folgenden Formeln die kovariante Ableitung stets auf das unverformten System 
bezogen ist, lassen wir den Index ( von jetzt ab fort und schreiben d statt (Wd. 

Für die Ableitungen der in Frage kommenden Größen hach t ‚ genommen für t=0, 
gelten dieselben Betrachtungen, auch hinsichtlich der Bezeichnungen, wie sie beim Kreis- 
zylinder angestellt worden sind. 


Die Berechnung der ö*.; führt hier zu den Ausdrücken: 


1 ( „„dDöv Döw u 
3 (6 su +% za rad dm 


uo 86 u 3 
urn atom! -%) 


Ö ge = 290 dw. 
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= 
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Bei der Konstruktion des Längstkrafttensors und des Momententensors ist es jetzt aber 
zweckmäßig, zunächst den Schubmodul @ und die Querdehnungszahl m als elastische Kon- 
stante zu wählen. Weiter setzen wir wieder 

ur; 
12a 
und führen die Invariante 


sowie den invarianten Differentialoperator 
D: 


L= a?g', dr dp +2 


ein, der in engem Zusammenhang mit dem Differentialoperator von Beltrami steht. Man 
erhält so: 


u 


ste anlarnl|gie dd 429% or am+2 _e+2sm] 
( Aue nn + )I go dg:  Dg* r=o9,on “m ZREFERD 


ME 1, Ö, \ 
2k\arge ddp +5, w+ een L(ö m) | 


Br u.dor, , dr" ; Ö; 
sm! =2Ghak| di zn + Y +25 dw+2,. (de+2öm) 





a: ug Bm ur 0; I 
a Io) dg‘ dr 1 on m 1 Liör)|j 


Die hochgradige Symmetrie der Kugel macht sich in dem vollkommen gleichartigen Aufbau 
von Längskraft und Momententensor geltend. 
Wir sind nun in der Lage, die Stabilitätsgleichungen aufzustellen. Es zeigt sich, daß 
die Grundlösung (8.5) in die ersten beiden Gleichungen gar nicht eingeht, da die Differenz 
joul 2 ji ul u 
IS : T 
FE Be 7 


für T; = — Tö,, identisch verschwindet. Die Spur des variierten Krümmungstensors ist gleich 
ee L(öw), so daß die Stabilitätsgleichungen (5.2) die einfache Form annehmen: 


„u nd 
DT; 


R a ren an Ve 
ad +3 öM%- _r L(ön)=0 

Wir führen in die ersten beiden Gleichungen die Ausdrücke (8.9) ein und erhalten: 

DT, 

d y“ 


— J alte d’öv 2 dv „9 öw 2 d ’ | 
=@hjttk) KO dr de ddr tr" di ur“ Fri ai 
Mi d D? öw dw 1 Ri) L(ö w) ' | . ” E 
—- dk 2 % ’ Du 
= k ß Io) d gq“ (da 5 } ( 
- u u ; a dör: 
Es empfiehlt sich in die Gl. (8.11), außer der bereits auftretenden Divergenz de = de’ 


N) g* v m—1 Ö q’ 
den 


Drehtensor: 
(8.12) 


1 
u 3 dgq“ dq 


5 ze ar) 
High d gy’ 
einzuführen. Wir geben zu diesem Zweck den Gl. (8.11) die Ferm: 

„j2m 2 TRONEIR® ou Pövi: D’öv“ 
A+Kk) m i er t dgl" de dt dig“ 


9296 m dv’ dw“ $- d? ö v" \ EP UN... Ö L(öw) 
dgl dg“dg* di dgef m—|1 og’ 


“es d? döw d? RE:ET) döw 
BE E24 VERRBRÄRRESER, | Ncdlioröiiuchl: BERRRRERESTS SURNENEEN Zune nein | HERGe- Banclbnec ET 
2rla Io) kr dgq® ) Trier] | ( 


( ö ", Ri) Ö ) 
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Die Differenzen der zweiten kovarianten Ableitungen eines kovarianten ‚Vektors sind 
nach Gl. (6.14) zu bestimmen. Danach gilt: 


d? dv“ ie d’ dot 
dDq"dgt Ddgtdg“ 
d? (® ö 2) 3 d’ ( öw dw 


2 hi — 1- _ —— || = a? 4? _ 
a Alrzr. dgl dg* dg“ dq ) a Io) , Au dq” 


=gn (0) we 10% 


Die Komponenten des Riemannschen Krümmungstensors werden speziell für die Kugel dar- 
gestellt durch: 


(0) 962). 


OR n! I (99 g0 u— gen Io 


Die betreffenden Differenzen von kovarianten Ableitungen gehen damit über in: 


d’ dv“ dv? öv“ 1 4% ‚(3 ö e\- d’ $ ö "= döm 


dg“ dgi _ dgl dge —— a? Ö v2 a? go) \d g“ dg* gi | dg* dg" d gi d gi ' 


Ferner gilt: 
ou D’öva Döut „ufpdröv dv’ övuı 


Ko) du dge Dgrdg“ = Ko) dgudge dgide| 
Br ou d° 0) ® 1 d? Ö Vu \ & j (0) v E a | 

190 de dg“ —pgdgi rt ORien R. ua) ö vor 

FF: ai A 1 

I) d g? d gq“ d g* a? 


ud OWDur 1 döw. 2 
=2g6, dg +— övi d.gqe 


dv 


==8 -+Zz6n. 


a? 


Die Gl. (8.11a) erhalten unter Berücksichtigung obiger Gleichungen eine überraschend 
9 


; EL 2 m 3 IE 
einfache Form. Nach Division durch ee... und mit es ergeben sie sich zu: 


jpdp(de+2 2öW) 


142 -0... (813) 


zovi 


dgq® r7 urr; 


dom, 1 6 :] „2 L(öm) 


®° 
+atn)| dg. 


DöM; N 
Um nun die dritte Stabilitätsgleichung aufzustellen, braucht man nur — d 


(8.11) und (8.13) abzulesen: 


aus den G1.(8.9), 


u 


DM; _Ehak|, dlde+2dm) | dv, 1 er) en 


Fe ey us Hi +07, er 


und erhält unter Berücksichtigung von (8.13) zunächst: 


döM; 7 :3Llöw 
En, en VA 
d ge I-#itk dgQ 
} DOM". ) 
Die Invariante dgidge ist nun nach Gl. (8.14) leicht zu berechnen, wenn man sich an die 
Definition des Operators L erinnert. Die mehrfache Anwendung des Operators bezeichnen 
wir durch die entsprechende Potenz von L und können dafin schreiben: 


PöM Ehk 1-—k i 
dgyrdg“ = - r da ig 2] + k | L? (ö w) — 2 L (ö „))  . a . 2 > (8.15). 


Die in der dritten Stabilitätsgleichung auftretenden ersten Invarianten 6 T}; (des Längskraft- 


tensors) und öM7; (des Momententensors) ergeben sich zu: 
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9L 
STI=GHl+M 2(0e+23m)+,,_de+2öm) m Lim) 


— ER + ui+hlde+2m)—kL(öm)) 


De LTE 2(de+25W) — L(ön). 


Setzt man noch zur Abkürzung 
T1—")_ 


ir en. 
so lautet die dritte Stabilitätsgleichung (8.10) 


ke’ (öm) —2 Lion» ++) +böe+25m)+e1 Liön)=0 . . . (8.16). 


Wir kehren nun zu den Gl. (8.13) zurück. Die linken Seiten stellen die kovarianten Kom- 
ponenten eines Nullvektors dar. Die Divergenz und der Drehtensor dieses Vektors ver- 
schwinden natürlich, d. h. aus (8.13) folgen die Gleichungen: 


(1+k){L(de+2ödmw) — 2(de+ 2ömw) — (1 — »)[L (dw) — (de +2 öm)]} | 


8.17), 
— k[L? (öw) — 2 L(öm)] — 0 ae 
von‘, 6 W. u 2 om öweu 2 , 

= ÖWui = 9: baw. Io) = de ae dgl == +-zd5m=0 . . (818). 


dd dgide "a a 
Das Bemerkenswerte an diesen Gleichungen ist, daß (8.17) wie (8.16) nur noch de und öw 
enthalten, während (8.18) eine Differentialgleichung für den Betrag des Rotors allein ist. 
Durch eine Umrechnung, die wir hier unterdrücken, kann man (8.18) auf die Form: 


I ) us 
a’ rd 0 Eee (8.18a) 


bringen. Die Gl. (8.18a) ist aber nicht etwa L(öw,,), denn öw.ı nimmt zwar bis auf das 
Vorzeichen nur einen Wert an, döw,=—öw,,, ist aber in (8.18a) als schiefsymmetrischer 
Tensor zu behandeln. Um nun aus der obigen Gleichung eine Invariante zu bilden, multi- 
plizieren wir sie skalar mit einem geeigneten schiefsymmetrischen Tensor mit den Kom- 
ponenten: 


1 
lt, == ZT lt, 


Er hat wie der Maßtensor die besondere Eigenschaft, daß seine kovarianten Ableitungen 
verschwinden. Wir können deshalb schreiben: 


en . uk 
re +20. es ee Se 


eine Gleichung, die bis auf den Faktor 2 gleichbedeutend ist mit 


TE a 


ar 
wobei do = ——? der absolute Betrag des Rotors ist. 


0 
Sind q, und g, Poldistanz und Länge eines Punktes auf der Kugel, so hat in diesen 
Koordinaten der Operator L die Gestalt: 
Ee 
sin?q,09 


BE: KR 
uhr TH, 7 Yin 


Es kommen für dw nur die Lösungen der Differentialgleichung L=0 in Frage, die auf der 
ganzen Kugel regulär sind. Außer der trivialen Nullösung ist nur noch die Kugelflächen- 
funktion 1. Ordnung: 


WB a 


Y,=4Acosq,+sing,(asing,+tbcosq.). :. . 2.2.2... (821) 


eine derartige Lösung von L=0. 
Durch den zugehörigen Drehtensor: 


Ben Pe er ae er RR 
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ist aber keine Verformung, sondern nur eine starre Drehung der Kugel bestimmt. Wir 
wollen etwas näher darauf eingehen. Aus der allgemein gültigen Definitionsgleichung der 
starren Verschiebung’): 

dDövı  dDövu 


dg“ er 


läßt sich die Gleichung: 
DÖWwur 


BE 5 ‚OD. = 94) 
dr AR 5 5 MENGE En 


ableiten. Durch Divergenzbildung erhalten wir: 


dDöw‘; 
de 
Die Gl. (8.25) ist bei zwei Veränderlichen mit (8.24) gleichwertig, da diese in dem Falle 
auch nur zwei verschiedene Gleichungen umfaßt. Nun ist — g’“ R’,„a=K 6; (K=Gaußsches 
Krümmungsmaß der Fläche), so daß (8.25) übergeht in: 
d ö [077 N 


a t+Kin=0.....2.2222... @ 
9° 


RE te een 


; 2 i Döwr* TER D?Ome* , ! 
Bildet man nun die Divergenz von - = + Koöv* und beachtet, daß di = identisch ver- 


{ Die _. ; d A . 
schwindet, während de= 5 a gleich Null ist, weil es sich um eine starre Verschiebung 
handeln soll, so folgt als weitere Bedingungsgleichung: 


oK,, 


die bei der Kugel als Fläche mit konstanter Krümmung, K-oı identisch erfüllt ist und 


z. B. für eine beliebige Drehfläche sofort erkennen läßt, daß nur die Drehungen der Fläche 
als Ganzes um die Rotationsachse in Frage kommen. Andererseits liefert der Rotor von 
(8.26) die Differentialgleichung (8.18), von deren Lösungen nur die eine starre Drehung be- 
deuten, die den aus (8.23) und (8.25) folgenden Bedingungen: 


WERTEN“ 
Ar Fmu 7 


genügen. Man überzeugt sich leicht, daß die Lösung (8.22) diese Bedingungen erfüllt. 

Wir sind auf die starren Drehungen etwas ausführlich eingegangen wegen der be- 
merkenswerten Folgerungen, die sich bei zwei Veränderlichen aus der Gl. (8.24) ziehen lassen. 
Eine starre Drehung ist ohne Einfluß auf die Formänderungen und kann als nicht vorhanden 
angesehen werden. Wir dürfen deshalb den Drehtensor bei der geschlossenen Kugelschale 
gleich Null setzen. Es bleiben also nur noch die beiden Gl. (8.16) und (8.17) für de und dm. 
Diese lassen sich nun in besonderer Weise behandeln; denn wegen des Verschwindens des 
Drehtensors ist es möglich, den Vektor öv; als Gradient einer Skalarfunktion einzusetzen: 

OU 

A. van der Neut führt die Integration seiner Stabilitätsgleichungen mit einem dem obigen 
gleichwertigen Ansatz, der als vereinfachende Annahme eingeführt wird, durch und zeigt 
nachträglich, daß der allgemeinere Ansatz: 

U oV 
v, Ig® öv,= dq? 


der auch in der Biegetheorie der Kugelschale Eingang gefunden hat, sich hier auf V=U 
reduziert. Die Brauchbarkeit des Ansatzes (8.28a) beruht darauf, daß es bei der Kugel mög- 
lich ist, eine Gleichung allein für öw, d.h. für V— U aufzustellen. 

Für unser Stabilitätsproblem genügt es und ist es auch bequemer, mit de und öw zu 
rechnen. Durch Addition der beiden Gl. (8.16) und (8.17) erhält man eine einfache Beziehung 
zwischen de und öw, die es uns ermöglicht, de zu eliminieren. Diese Beziehung lautet: 


17) Vgl. [2], [3]. 


0 ER 


ö (8.288), 
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+) Lie +20m) (1 - +) L(öw) +; ° Lö) = ee 


Wir wenden nun auf die Gl. (8.16) nach einmal den Operator L an und erhalten unter Be- 
rücksichtigung von (8.29) eine Differentialgleichung, die nur noch öw enthält: 


1+k 
k£{L’ (dw) —2 LP’ (ön)y ++ k)(l—v*) L(öw)+o; Tr (dw) —(1+») L(öw)}—=0 . (8.30). 
Von den Lösungen dieser Differentialgleichung kommen wieder nur solche in Betracht, .. die 
auf der ganzen Kugel regulär sind. Nun läßt sich die partielle Differentialgleichung 6. Ord- 
nung für dw in bekannter Weise durch den Ansatz: 


en 


in drei Differentialgleichungen von der 2. Ordnung aufspalten, wobei der Parameter A der 
kubischen Gleichung: 


„4; „‚„i+rk l+k 
4 +24+ k ( N A  ü+l+n =(0 . . . . . (8.32) 


zu genügen hat. Die drei Wurzeln dieser kubischen Gleichung sind: 


R ü l+ko l+k o : 1—»? 
RR jagt) f at) - k 


Aus der Theorie der Kugelfunktionen ist aber bekannt, daß die Differentialgleichung 
nr rt 


nur dann eine auf der ganzen Kugel reguläre Lösung besitzt, wenn der Parameter A die 
Bedingung: 
iA=n(n+1l)—2 ae 


erfüllt. Die Wurzel 4,=0 kommt hier nicht weiter in Betracht. Die kritischen Werte 
von 0, für welche die Grundlösung (8.5) instabil wird, hat man also aus der Gleichung: 
i+ko I+k o 1 --»? j a 
3-18 Vi Hr) u =n(n+1l)—2. . . . . (835) 
zu bestimmen. Nun ist die rechte Seite von (8.35) reell, ja sogar ganzzahlig. Offenbar muß 
dann stets die Beziehung: 
es 0 )= ei 1—»? 


I-E2R +r k (8.36) 


gelten. Das Gleichheitszeichen liefert eine untere Schranke für den kleinsten Eigenwert, 
nämlich: e 

A — k) 8.37 
e=277, WR v’) —v te 


< 


Für diesen Wert von o ist die linke Seite von (8. 36) von der Größenordnung ”. also bei 
kleinem k eine große Zahl, die nur wenig von einer ganzen Zahl abweicht, so daß praktisch 
die untere Schranke (8.37) gleichzeitig den kleinsten Eigenwert darstellt. Der kleinste kritische 
Außendruck ist also gegeben durch: 


r0- —k TEE 
oder bei Vernachlässigung von Gliedern höherer Ordnung: 


(8.38). 


Wir wollen nun dasselbe Problem noch einmal mit Hilfe des Ansatzes (4.1) für die Form- 
änderungsarbeit behandeln. Es zeigt sich, daß wie beim Kreiszylinder auch hier die Stabilitäts- 
gleichungen komplizierter sind. Es genügt wieder, die ’; bis zu den Gliedern zweiter 
Ordnung zu berechnen. Diese sind für die Kugelschale gegeben durch die Ausdrücke: 
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ER u: REN s : 
heb,+tu+ DI W+tu+tsin’,+40,0,+0)) 
1 (4 ” . ’ . 
=zW, +sin’q,öd,+2v,v,+2sin?’gq,®,v, 


+2sin q, c08 q, 0,9, +%,%, +sin? q, 0,9%, +0, 0, 


1 v 

4 

7} on v‚+v, +ctg 4 v,+ 57 v) + sin? q r v; 
Pi 


+ tr (ctg? 9,7 -1) v; + 4 ctg Yı dd; + 4 d, v; + 4ctg gı v, v, 





Die Berechnung der skalaren Produkte ’j; j, und ’ee bis zu den Gliedern dritter Ordnung 
geschieht nach Gleichung (7.17) und führt zu den Formeln: 


ei 1/( a Dr 
ee tl | 
v,+sin?’q,t, 


2 . 2 . 2 ‚ G 
sın u Pe vv 
Ai u q, 9. + sin’ q,6,) sin’g, [id 


+ sin? I v, d,+ (v, — sin? 9ı ö,) eig q,dı — ö, v,] 
eenteten ei, ++ eig qıd,) + +v,0, 


or -3[; a 


249 7. sin’a 2 \sin? gg, 


- — sin? )) 
sin?q, sin?’gq, 4. v3 


ee=(, +, +27, +egq, dv ld, +9, +20, +ctgg, 2|ö 


en 
+4 +24 +2V, 0, + 2 q, %T sin?g, 





Der Ausdruck für die Formänderungsarbeit der Kugelschale läßt sich jetzt angeben, da die 
Krümmungsänderungen bereits berechnet und aus den Gl. (6.11) zu entnehmen sind. 
Man erhält: 


n 


n 2 
a=6n| \ke. +0.’+3- -w + sin? q, d,)? 
v6 


2 sin? q 


[4 1 . [4 
+(W,+v,+ctg 4 fe @#,+%,+2v,+ etg q, v,)” 


.. ' ’ 1 (4 . . R 
+k (%,— ö,) + erg 3 (ö v,— ctg 4,d% — 2 (v, + sin? gı ö,) 


„ 


Ha + ctg 9 a v, er ctg 9ı v,) 
1 v ; , » 
+ tet) tree) | 


(ö, +]; + +0,1,-0% +5 3 la q, — sin? 1) | 


e,+sin’q,® r R) 
ra ’v v,%,+ sin? q, Ö, ®%,+(v, — sin? g, d,)ctg 94, %, —d; v,) 
er v) ” 
-, tu, +dtg,uv)ua tu +00, + sin? qı sin? qı 





1 2 EEE 
gl q, — sin’ 9: ö) Ui | @,+® + 2», + ctg gı v,) 
E) . 


” 4 [3 . V® . Un 
tät +2 tn = rn Ne sin g,dg,dg, 
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Die Berechnung der Arbeit, die der konstante Normaldruck mit Jen Komponenten p; =— p£i; 
bei einer virtuellen Verschiebung leistet, geschieht wie beim Kreiszylinder. Die mit yg 
multiplizierten Komponenten der Flächennormale erhalten hier die Werte: 


Vsu;=alttrPsing to), —hl+r)] | 

v9; =a(l+v,)sin,+v)ß,%,—v,(l+0,))] 

vs; =e®(l+r’sin, +eY)[lito)d+e)—e,d,] | 
so daß die virtuelle Arbeit durch die Formel: 

ö Ay er | | d+ v,)” sin (4: + v,) [tv, v; A Ö; (1 + v,)} öv, 
00 
+, —v,(l+ö01ö0, + +0) hv)— v0 ör,]a’dg,dg, 

gegeben ist, die als die erste Variation des Integrals: 


n 2n 
. 


1 ' ’ ". ' 
4A,= ua; ap\ \ A+rP’ıl+H)l+Vr)—v,d}atsin(g, +er)da,dg, . (8.42a) 
v6 | 
angesehen werden kann, d.h. der Belastung läßt sich ein Potential zuordnen. 
Der Grundlösung (8.5) entspricht die auf den ursprünglichen Radius a bezogene Radial- 
verschiebung: 


.(0) ap Eh \, ] ©_ 
“= "r5pür ml= >); Pet Pt .... . Eu). 


Die Eulerschen Gleichungen des Variationsproblems 
en Vu 
haben die Grundlösung: 
vi =—b; vs” =0; % N (8.44). 


Vernachlässigt man bei der Berechnung der konstanten Radialverschiebung wie beim Kreis- 
zylinder höhere Potenzen von v,, so erhält die Konstante den Wert: 


vu) — ö ap 

’ 2D(l-+») 
Wir berechnen nun die zweite Variation der potentiellen zen U=A—A,, die der Grund- 
lösung (8.44) zugeordnet ist, bezeichnen die Variationen von v v, mit %,v,w und erhalten 
folgenden Ausdruck: 


(8.44). 


19 Us, 


2n 


ö° u=6n\ Ni + m)’ -(# + sin? q, d)? 
PR 


+3 sin? q, 


+(® + m +.ctg q, u)? 4 (“+9 +2m+ctgq, u)’ 


2 ‚ 1 we .\ 
+k (& — u’ + sintg (" — ctggqg, w — 3 u — 5 sin’q,6) 
1 


„ 2 
+ q, -+egg,m—v—ctgg, u u) 


w 


2 
u —9— cig a) |a’sin q,dg,dg, 


+ lt eig uote qı 
ap| ER : 
+ 4 w+R+2u OF zineg, 


1? 
sin? q, 


1a’ sin q,dq, dq, 


+6(& +0 +ctgq,u)w + 10 w? — w 





OR ee 1 
+5 uv—uo+Wt+r)cg,u— Zu 


Die der zweiten Variation zugeordneten Jacobischen Differentialgleichungen lauten: 
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dl. Bi a | 
ERLERNEN 


—k ti Pier lir+ eig q,w+ - u v9 —ctgq,u u) sin 4) 


sin? q, 


v— egq,w— ZU 


HERE 3 0X SB TERN 
+sing, 3a, rer (# + sin? q,®) sin’q, (ni 


sin? q,5)} - Nitwtelgqut —Wr+rv v+2mw+ctgg,u) 


w 
- kI-—— +c w— v — ctggq,u 
sin?g, tg q, 8 4ı 


e w : 
+ an (+ etgq,w+ sintg, u—v—cigg, «) cos q, 


eilig lernt er een) 


\ 


i an = 
+sing, 36 ui is sin? q, -+3ctg q, w 


9 
++ oetgqu—ul sin q1=0 
1 z ? = D:. Re 
b) Y | FALL + sin? g,ö) — ki — ctggq, mw — 3 u — 7 sin? q, ö) sın 4 
Hunse tw teigg, u + \ ı +v+2mw+ctgg, u) 


—k sintg, -t+etgy,w— vu —ctgq, u 


2 
+3 ‚(® +ctg,w + — U —r —ctggq, «) 


m — sin? q, 


Bi RN 
tor rn ( — 3 i) sin al En sın q, 3 ® + 3w 


‘ 


+3 lite; |} = 0 


2 ap 5:2, w 
c) Ic 97% m u . —u—v—ctgq, u «) sın a 


dd, Isinq\ eig q, w -g4-5sin’g, ) 


1 w 


0 %g slaor: qı +egqm—r—cetggu 


PA 


+ fü +cgq,w+- — (ü+v+ctgg, )) 


sin?q, 


d w 
3 (lie sims q, tes, w—v— ctggq,u 


PrA 


1 
ten ‚(ö+ AB + si; 2, — U —-d —ctgq,u u) cos 4, 


a a ae 
sin?q, dq, w— eg, w—- gü— sing, 2 


+ atom ++ gu + te +2 mw-+.ctgg,w)|sing, 





+ an let ent -+10m+3 (di +ö+ tg q,u)} sin q, = 0 





. Mech. Math. Mech. SR BER 
ani 1943 5 Ba Nr. 3 Juni 1943 ‚Kleine Mitteilungen 179 
u... 





Die Gl. (8.46) sind in der vorstehenden Form noch recht unübersichtlich. Unter Ein- 
führung der Größen e, » und des Operators L formen wir sie so um, daß sie möglichst den 
Gl. (8.13) und (8.16) entsprechen. Man erhält dann das System: 

° s 
%q,' 


a) (1+k)(e+2w) —kL(w)+o(e+3m)} 


IHR) (i ERROR 2 - 


; [| ® 
in \sin q, 


A) 
b) . 


dq (1+k)(e+2w) —kL(w)+o(e+3m)} 


+(1+%k)(1—r) {osing, +singe—w}—=0 
c) k£I? (w) — L(e+2w) — (1+v) L(w)} +1 + r)(1+k)(e+2w) 
+o{L(mw) +8w+3e} =0 
Die Grundlösung geht in diese Gleichungen etwas komplizierter ein als in die ent- 
sprechenden Gleichungen vorher. Insbesondere erscheint sie auch in den ersten beiden 
Gleichungen, aber in so einfacher symmetrischer Form, daß auch hier die Gl. (8.19) gilt, wie 
man sich sofort überzeugt. Die Gleichung, die der Formel (8.17) entspricht, lautet nun: 
(I+k){L(e+2w) — (1 —») L(w) -—(1+r)(e+2m)} | 
— k{L’(w) — 2 L(w)} + ot{L(e+3m) — 2(e+3m)}= 
Aus (8.48) und (8.47c) folgt durch Addition die Gleichung: 
L(e+2w) —(1—v) L(w)+o{2L(w)+L(e+2mw)+(e+2mw}=0 . . . (8.49). 
Diese Gleichung kann wie (8.29) dazu dienen, (e+2w) aus einer der Gl. (8.48) oder 
(847c) zu eliminieren, denn in ihr ist zumindest o (e+2w) zu vernachlässigen. Die weiteren 
Überlegungen gehen denen der ersten Stabilitätsuntersuchung parallel. Die in Frage 
kommenden Wurzeln der kubischen Gleichung erhalten die Werte: 


_1+k@-ı [1+k8—») ’ 1— . 
2 R 1+] | 3% +) - — "za ).(3—r). . (8.50). 





(8.48). 


Die Unterschiede gegenüber den früheren Werten sind bedeutungslos, da k (3 — r) gegen- 
‚2 


1 
über 1, und auch noch za -r)(3— vr) gegenüber ( —— vernachlässigt werden kann. 
416 
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Ergänzende Bemerkungzuder Arbeitvon 
D. Küchemann: „St Örungsbewegungen 
in einer Gasströmung mit Grenzschicht‘‘. 

Die Methode der kleinen Schwingungen hat bei 
der Untersuchung der Stabilität laminarer Grenz- 
schichten wertvolle Ergebnisse gezeitigt, welche 
im Hinblick auf das Problem des laminar-turbu- 
lenten Umschlags die Frage der Turbulenzentste- 
hung durch kleine Störungen bis zu einem ge- 
wissen Grade geklärt haben. Erstmalig hat nun 
D. Küchemann (198) die Frage nach dem d 

A 








Einfluß der Kompressibilität einer Gasströmung 
auf die Stabilität der laminaren Reibungsschicht in 
Angriff genommen [1}*). Aus Gründen der rechne- 
rischen Einfachheit legte er seinen Untersuchungen 
nach dem Vorgang von Lord Rayleigh [2] ein EIN 
aus zwei Geradenstücken zusammengesetztes Ge- Bild 1. Geschwindigkeitsprofil der Grundströmung. 
schwindiekeits 1 si i o ». (Das 
mag er 7 st u auf den Standpunkt stellen, daß in der Küche- 
herigen Stabilitätsrechnungen bei inkompressiblen Mäannschen Arbeit als erster orientierender 
Grenzschichten bedenklich stimmen, da sich die rich-. Schritt ein Vergleich des Stabilitätsverhaltens 
tige Berücksichtigung der Krümmung des wirklichen I"gendeines — mathematisch möglichst einfach zu 
Profils für das Resultat der Rechnung als sehr behandelnden -— Geschwindigkeitsprofils mit und 
wesentlich erwiesen hat. Allein. man kann sich ohne Berücksichtigung der Kompressibilität voll- 
2 n zogen werden soll, und daß es daher für diese 
1) Die Zahlen in eckigen Klammern verweisen auf das Rechnung gleichgültig sein darf, ob das Profil ein 
Schrifttumverzeichnis am Schluß der Arbeit. wirklich vorkommendes Grenzschichtprofil so gut 
12° 
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approximiert, daß die Resultate auf dieses an- 
wendbar sind, oder nicht.) Für die Störungsrech- 
nung wurde die Reibung vernachlässigt. Aber 
selbst bei diesem so vereinfachten mathematischen 
Modell setzte die rechnerische Kompliziertheit der 
Entwicklungen den Untersuchungen vorläufig enge 
Grenzen. Küchemann beschränkte sich auf die 
Bestimmung der möglichen indifferenten 
Schwingungen. Die besonders interessierende 
Frage nach der Berechnung oder auch nur nach 
der Existenz angefachter Störungen blieb 
offen. Man weiß, daß die inkompressible 
Strömung mit demselben Geschwindigkeitsprofil 
stabil ist. + 

Es soll nun im folgenden lediglich auf eine Kon- 
sequenz hingewiesen werden, die man bei Hin- 
nahme der Küchemannschen Vereinfachungen 
und Vernachlässigungen aus seinen Ergebnissen 
fast unmittelbar folgern kann: Aus den nachge- 
wiesenen indifferenten . Schwingungen lassen sich 
durch einen |Grenzübergang gewisse angefachte 
Wellen herleiten, die zwar nicht zu der Klasse der 
von Küchemann angesetzten Partialstörungen 
gehören, aber als Lösungen eines allgemeiner an- 
gesetzten Störungsproblems ein lineares An- 
wachsen der Schwingungsamplituden mit der Zeit 
aufweisen ?). Inwieweit diese Lösungen physi- 
kalisches Interesse besitzen, muß vorläufig 
dahingestellt bleiben. Die folgende Herleitung der- 
selben ist so abgefaßt, daß ein gewisses Bedenken 
in dieser Hinsicht sich aufdrängt. Ich komme an 
gegebener Stelle ‘darauf zurück. Da dieser Lö- 
sungstypus auch bei der Untersuchung anderer 
Geschwindigkeitsprofile auftreten wird, vielleicht 
jedoch ohne so deutlich sein Wesen zu offenbaren 
wie in der hier gegebenen Herleitung, dürfte die 
folgende Mitteilung erwünscht sein. 

Wir legen unseren Betrachtungen einen gegen- 
über dem Küchemannschen Ansatz fortschrei- 
tender Störungswellen allgemeineren Störungs- 
ansatz, in welchem die Gestalt der Zeitabhängig- 
keit dieser Störungen noch offen gelassen wird, 
zugrunde. Im übrigen stimmen die im folgenden 
gemachten vereinfachenden Voraussetzungen mit 
jenen der Küchemannschen Theorie überein. 

Die laminare ebene Grundströmung besitze die 
Geschwindigkeit U, welche parallel zur z-Achse 
gerichtet sein und, wie bei allen Störungsrechnun- 
gen, von denen einleitend die Rede war, aus Grün- 
den der rechnerischen Durchführbarkeit allein von 
der dazu senkrechten %-Koordinate abhängen 
möge. Dieser Grundströmung werde eine kleine 
Störung überlagert derart, daß die Gesamtströmung 
durch ihre Geschwindigkeitskomponenten x und ® 
in z- bzw. y-Richtung und ihre Dichte o in folgen- 
der Gestalt gegeben sei: 


ur, yd)=U(y)+u(y,t)eiex 
v(z,y,t)= uly,t)eiex 2: ME 


eo, y,)=ant1+o(yt)eiex, 


Dabei bezeichnet o, die als konstant angenommene 
Dichte der ungestörten Grundströmung. Mit reel- 
lem @a—2x/A bedeutet A die Wellenlänge der an- 
gesetzten Partialstörung. 

Die Eulerschen Gleichungen, in welchen der 
Druck p durch die als konstant angenommene 
Schallgeschwindigkeit a vermöge a?=dp/de er- 
setzt wurde: 


du Du d%u_ ade 
Bd Tau Ver Er e 
dv dv dv_ ade 
FLäMle 7 Ad Peer E77 


2) Diese Bemerkung entstammt einigen Rechnungen, 
die ich im Herbst 1937 auf Anregung von Herrn W 
Tollmien durchführte. 








und die Kontinuitätsgleichung 


do do d0 fe u, dv 
-+u U — =—o|-- 

set + dy 7; Er (8) 

liefern, wenn wir nur die in 4,© und ® und deren 

Ableitungen nach % und ? linearen Glieder be- 

rücksichtigen, das Gleichungssystem 


I ac : 24 

5 tie lurl v=—iaa?d 

DE?) “ v5 

: jaU%  ). &- i 

3, riealv | *.y „... (die 
ES EEE EN... 

34 riale a 7 


zZ 


Hieraus gewinnen wir durch Elimination von 
und v die folgende Störungsgleichung für o: 
0°0 

oy? 


do 00 ee 
un "zarte u tel 
2 „00 r 80 
2 ze ei D 2 Tr 
+ a? (a nt a 7° 
+ ia! (a -_UNUO—=0 


Für die von D. Küchemann angesetzte 
Störung in Gestalt einer in z-Richtung wandern- 
den ebenen Welle ist insbesondere 


u(y,t)=u*r(y)e-iPt | 
vyu=urWW)e-ift * ,, (6). 
0 (y,t) = o* (y) e-ißt 


Mit diesem Ansatz geht unsere allgemeine Stö- 
rungsgleichung (5) in die von Küchemann an- 
gegebene gewöhnliche -Differentialgleichung 2. Ord- 
nung für o*(y) über, in der fla—=c gesetzt ist: 


m U’ r 3 BER 2 
A ra +] a er=0 . (m. 


Der Realteil f, des Parameters ß bedeutet die 
Kreisfrequenz der Störung, f; das logarithmische 
Inkrement derselben. ; ist für angefachte Schwin- 
gungen positiv, für gedämpfte negativ. Der Real- 
teil c, von c liefert die Phasengeschwindigkeit, 
mit der sich die ‘Welle in z-Richtung fortpflanzt. 
Für den in (7) auftretenden Differentialoperator 
führen wir die Kurzbezeichnung D ein: 


2 ’ E58 2 
RE. OB » A. 3 + 1) . I 


Das von D. Küchemann zugrundegelegte Ge- 
gen hat folgende Gestalt (siehe 
ild 1): 


uo=fl) für An 


U=ybfür y>b 


(r=coonst). Nach Einführung der neuen dimen- 
sionslosen unabhängigen Veränderlichen £ = y/b 
—- c/U, und der Machschen Zahl M = U,/a geht der 
Operator D (nach ‚Multiplikation mit 5°) in den 
speziellen Ausdruck 


(9) 


d? 2 d : i 
Derazare® {M? 2? — 1}. (10i) 
innerhalb der Grenzschicht und 
d?. ey 
== —__ \2 2 rn a 
Da= + (ab) {m (1 7) 1} . (10a) 


außerhalb der Grenzschicht über. 


Die Randbedingungen des Problems er- 
geben sich aus den Forderungen, daß 
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1.v—0 an der Wand y=0, 
2.o am Profilknick y = 5b stetig und 
3. auch ® in y—=b stetig 


ist. 

Stellen o, una 2, ein Fundamentalsystem der 
Differentialgleichung D;o* — 0 dar, so lautet die 
allgemeine Lösung in der inneren Schicht: 

An. tan. .- - -»-» (IN. 


Die Stelle {=0 (die jenem kritischen Wand- 
abstand 4 entspricht, in welchem die Phasen- 
geschwindigkeit der Wellen mit der Grundströ- 
mungsgeschwindigkeit übereinstimmt) erweist sich 


als Nebenpunkt, die Lösungen verhalten sich also 


in der ganzen Grenzschicht regulär. Dasselbe gilt 
für die Geschwindigkeitskomponenten [1] S. 209 
bis 210. 
Außerhalb der Grenzschicht ergibt sich aus 
D, 0 * — 0: 
PB ii Bei... . 2. 


Am. 
N - m. AU: 
0 


Verschwindet der Realteil ZL, von ZL, so führen 
die Randbedingungen auf drei Gleichungen für die 
vier Integrationskonstanten A, As, B,, B, und 
physikalisch ergeben sich Störungen beliebig vor- 
gebbarer Wellenlänge und Phasengeschwindigkeit. 
Es handelt sich um von außen kommende Störun- 
gen. Das eigentliche Stabilitätsproblem ergibt sich 
für L.$0. Es werde dann ohne Einschränkung 
der Allgemeinheit Z, >0 angenommen. Dann muß 
notwendie Bz — 0 sein und man erhält Störungen 
in Wandnähe. Für die verbleibenden 3 Konstanten 
Ay Ay, Bı =B liefern die Randbedingungen das 
homogene Gleichungssystem 


4,0 + 4,0,’ =0 


mit 
L=avfı- (1 


für = — c/Uv (Wand) | 

A,0ı + 4,0, =BerL: | 

für t=t4=1-—ciV, | (14) 
’ ’ -1r schicht- 

Aı0, + 4,0, =—LBe Fr | rand) 

für t=ık=1-—clUs 


und die Lösbarkeitsbedingung ergibt zu jedem 
festen U, Eigenwertkurven a=a(ß). Die Natur 
und der Gültigkeitsbereich der von Küche- 
mann hergestellten Lösungen setzte dem Ver- 
such der Berechnung komplexer Eigenwerte ß sehr 
große Schwierigkeiten entgegen, so daß Küche- 
mann sich im weiteren Verlauf seiner Rech- 
nungen auf die Bestimmung der möglichen in- 
differenten Schwingungen (ß; — 0) beschränkte. 

Für die reellen Eigenwerte ergeben sich die in 
Bild 2 wiedergegebenen Kurven. Ich habe dabei 
auf Grund des mir von Herrn Küchemann zur 
Verfügung gestellten Zahlenmaterials eine andere 
Darstellungsweise gewählt als in der Kiüche- 
mar nschen Veröffentlichung. Der Grenzfall 
M=0 zeigt den guten Anschluß an die von 
0. Tietjens für das Geschwindigkeitsprofil (9) 
berechnete Eigenwertkurve des Falles der volum- 
beständigen Strömung, in welchem sich Stabilität 
herausgestellt hatte [3] °). 

Als wesentlich neues Merkmal im Verlauf der 
Kurven für M>0 gegenüber der Tietjensschen 
Eigenwertkurve des inkompressiblen Falles zeigt 
sich ein Umbiegen derselben nach Erreichen eines 
maximalen a-Wertes (minimale Wellenlänge A) und 


Rückkehr zum Ursprung. Zwar gestatten die 


3) Die Gleichung der Eigenwertkurve lautet dort: 
Bb 
Uo 


=ab-2(12e-2eb), 


12 
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Bild 2. Eigenwertkurven der indifferenten Schwingurgen 
für verschiedene Machsche Zahlen. 


Küchemannschen Lösungen der Störungs- 
gleichung nicht, die Eigenwertkurven für die prak- 
tisch in Frage kommenden Machschen Zahlen M 
in den Bereich größerer a-Werte zu verfolgen — 
es gibt dort auch weitere Äste der Eigenwert- 
kurven —, jedoch wird die eben hervorgehobene 
Eigentümlichkeit der Kurven sichergestellt durch 
die mittels asymptotischer Lösungen für große 
Machsche Zahlen berechenbaren Eigenwertkurven, 
von denen allein aus diesen grundsätzlichen Er- 
wägungen heraus jene für die große Machsche 
Zahl M — 10 berechnet wurde ®). 

Gerade das Vorhandensein dieser Maxima der 
obigen Eigenwertkurven a—a(ß) — nennen wir 
sie a—=a, — gibt uns nun die Möglichkeit, aus 
den zugehörigen indifferenten Lösungen angefachte 
Lösungen mathematisch zu konstruieren. Während 
für eine zu einem vorgegebenen M-Wert gehörige 
Eigenwertkurve zu jedem a<.a, zwei Kreisfre- 
quenzen ß, also zwei indifferente Wellen möglich 
sind, rücken diese, für @a>a, zusammen. Über- 
lagert man für @a<ao zunächst die beiden mög- 
lichen indifferenten Wellen: 


E (y,t) eiex— Ao*(y;d,) eilex—Pıt) 
+ B o*(y; Ba) eilax— Pat 
und wählt man B = — A, so ist 


e (y,t) eiex— A t0* (y;B,) — o* (y; Ba)} eilax—Pıt) 
— A 0o* (y; ß,) fei(ex— Pt) — eilax—Pıl)), 





Geht man nun zur Grenze a>, also auch 
ßı— ßa>0 (fi, f2> Po) über und hält man dabei 
A (fı— Be) =o als von -Null verschiedene end- 
liche Zahl fest), so wird in der Grenze 


& . a8r- . : 
(y.8) ein 00, 9 itor) a) FR RR 


4, Eine Neuberechnung nnsererseits ergab kleine, hier 
unwesentliche Abweichungen gegenüber der von Küche- 
mann angegebenen Kurve. 

5) Diese Annahme: o=const für den Grenzübergang ist 
physikalisch nicht evident. Es läge näher, statt o die 
Größe | A| als konstant anzusehen. (Anders liegen natür- 
lich die Verhältnisse bei erzwungenen Schwingungen, 
denen dauernd Energie zugeführt wird.) Auf diese Stelle 
der Rechnung bezieht sich die einleitend gemachte Be- 
merkung zur physikalischen Deutbarkeit der Lösung. 
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Dabei sind also f, und a, = a (ß,) die Parameter- 
werte des Extremums der Eigenwertkurve a — a (ß) 
der indifferenten Eigenschwingung (deren Ampli- 
tude 0 * (y; a,, Po) ist) bei der betreffenden Mach- 
schen Zahl M, und es bezeichnet o eine Konstante 
mit der Dimension (Zeit) ' 

Da der vollzogene Grenzprozeß legitim ist er 
daıs) 
ds 
existiert, stetig ist und in = ß, verschwindet, 
-— stellt (15) eine Lösung des Störungsproblems 
von dem angekündigten Anfachungscharakter dar. 
Man kann diese Lösung selbsverständlich auch 

direkt durch den Versuchsansatz 


setzt die Kenntnis der Tatsache voraus, daß 


o(y,t) =otR,(y) -itR,(y)}e-ißt (16) 
einer linear mit der Zeit angefachten Welle ge- 
winnen. (Der Faktor i bei R, ist hier nur aus 
Zweckmäßigrkeitsgründen gewählt.) Geht man 
nämlich mit diesem Ansatz in die allgemeine Stö- 
rungsgleichung (5) ein, so erhält man 


RY 
(U-—-c) | 


wa ar, I 
ii - u = IR, 
a! Pr ER? ben it DR, 


wo D der Operator (8) ist. Da diese Beziehung 
identisch in ? gelten soll, müssen linke und rechte 
Seite für sich verschwinden. Es ist also einmal 
DR.,—=0 und mit Rücksicht auf die Randbedin- 
gungen muß also 


BEN. 


DR, - 
‘ 17) 


sein. Damit sind auch bereits die möglichen Para- 
meterwerte a, ß auf jene der Eigenwertkurve 
a—.a(ß) von o* beschränkt. Die inhomogene Dif- 
ferentialgleichung für R, (nullgesetzte linke Seite 
von (17)) hat, wie man dann leicht bestätigt, die 
d o* 
GB --, 
bezeichnet ®). Die additive indifferente Eigenlösung 
mit der Amplitude k o* bleibt natürlich frei, da 
sie stets jeder Lösung der linearen Differential- 
gleichung (5) überlagert werden kann, und wir 
können ohne Einschränkung für unsere Betrachtun- 
gen k — 0 setzen. Die Randbedingungen, die man 
für R, direkt auf Grund der oben im Text ange- 
gebenen drei Forderungen an vo und v in = £y 
bzw. in = {x herleitet, stimmen mit den Aus- 
drücken, die man durch totale Ableitung der Be- 
dingungen (14) nach £ erhält”), dann überein, 


a a 
— 0 ist, 


Lösung R, +ko*, wo k eine Konstante 


und man gelangt damit gerade 


zu der Lösung (15) für die Parameterwerte a — a,, 
B = Bo. 

Es ist nur noch die Regularität der Lösung zu 
prüfen. Innerhalb der Grenzschicht ist, da ja nach 
der ersten Bedinzunz (11) o*(£; a, $) bis auf 
einen freibleibenden Normierungsfaktor in der 
Form 


+ (lu) RW). (19) 


geschrieben werden kann und dabei 


= 
ey 


6 
u -2Vv -2Vv ” fm 
= 2 a,L? » bay£? (20) 


v—( ö y 


=’) 
=0 


BR Bir Er 
6), Die allgemeine Lösung lautet A, = 3 3 +Cı 01 +C30s, 
aber da man R' (im) —=0 zu fordern hat (siehe unten), 
® 


ist wegen Fr ({y)=0 zu verlangen, daß C, 0; (pr) 
+03%(öw)=P, also Cı:C3=Aı:4a nach der ersten 
Gl. (11) ist. 

?) Das darf man ja, da diese Beziehungen gerade längs 
a=a(ß) identisch in 9 erfüllt sind. 


(!1], Gleiehunz (13): dort sind auch die Rekursi- 
onsformeln für dıe Entwicklungskoeffizienten an- 
gereben) ist, 


d 0o* 


R, = dd 3=3% 


’ “ “ = \ 
TO \sW)2a\s 


Dabei ist die Tatsache benutzt worden, daß 
dE : dp 1 da 


B wegen 
dB 3=3 dd |3=% 7 4B|8=3 


BlaU, und 
Außerhalb 
wird entsprechend, wegen 


- =—=( 
do I 0 


ist und ferner, daß nur &—y/b 

Sy ßlaU, in (19) von f abhängen. 
der Grenzschicht 
daL Ei 
dß 3 Fa 
ab? M? 
U,L 


ab?M?/ c 


U, L U, 


I=)’g 


(1 „Be LZ+-Ce u 


aus den Randbedingungen bestimmt 
do* 


R,= f 


wo sich C 
der Kenntnis, daß R, = ist, 
fa, 


dd 3 71 
BL] 
( me I — I } rec at 8 37 
lad Tadıls-a )erechnet®). Es 


bzw. mit 
sorleich zu 


ist somit R, und also auch © (y. f) innerhalb und 
außerhalb der Grenzschicht in y regulär (und am 
Knick ©= {x auf Grund der Randbedingungen 
stetig und stetig differenzierbar). 


Für die Geschwindigkeiten ist entsprechend (15) 
anzusetzen: 


u(yt)=o{F(y) --itur(y)yemidt \ 


| 


Wenn man zur Auffindung der Lösung zunächst 
versuchsweise den allgemeinen Ansatz (16), wie 
dies oben geschildert wurde. macht. braucht man 
für die Randbedingungen, die an AR, zu stellen 
sind, die Kenntnis von v. also von G (y), denn es 
soll ja v® an der Wand verschwinden und im Pro- 
filknick stetig sein. Nun gilt ([1], Gl. (6)) all- 
cemein 


(23). 


v(y,t)=0{G(y)—itv*(y)}e 


> 


ia? 


- ad: MEET. 
a(U—c) 


1 * (y)= 


. “ R u . P . 
Damit errechnet man G (y) unter Berücksichtigung 
der zweiten Gl. (4) zu 


G(y)= s {v* (y)+ia? R; (y)} (25). 
z 


1 
a(U 


(Außerhalb der Grenzschicht ist insbesondere 
U U, zu setzen.) Es wäre demnach von Rı zu 
fordern: 1) Rı =0 an der Wand. 2) und 3) R,_und 
R', stetir am Profilknick. Dies aber nur nebenbei. 

Hier interessiert uns nur die Frage der Regu- 
larität der zur Lösung (15) gehörigen Geschwindig- 
keitskomponenten» Es ist dabei nur die kritische 
Stelle 0 besonders zu prüfen. Bei Berücksich- 
t'eunz von (9), (19), (21) ist 


5 ia? j o*"”(£) o*’(£) 
(( = = - 
aU, 


) adbU.Eel 


aU,f 
| 
\ 
J 


+9 (ww) (dd) - er (Ew)o: (£ 


AR . a x 
8) | 3 läßt sieh demnach, nebenbei bemerkt, direkt aus 


{ 
den Randbedingungen für R, ins={xK berechnen. 
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Setzt man in (19) die Potenzentwicklungen (20) 
ein, um sie dann in (26) einzuführen, so heben sich 
hier gerade die negativen Potenzen von { heraus, 
denn nach steigenden Potenzen angeordnet ist 

ih. 


‚2 


Damit ist G@ (2) und also auch v (y,t) in {=0 
regulär. Aus der dritten Gl. (4) folgt dann un- 
mittelbar auch die Regularität von uıy,®). 

Es ist somit die angefachte Störung (15) eine 
Lösung des Störungsproblems, die allen aus 


physikalischen Gründen zu stellenden Regularifäts- 
bedingungen genügt. 


Schrifttumverzeichnis. 

{1} D. Küchemann: Störungsbewegungen in 
einer Gasströmung mit Grenzschicht. Z. an- 
gew. Math. Mech. Bd. 18 (1938), S. 207—222. 

| Lord Rayleigh: On the stability or in- 
stability of certain fluid motions. ‚Papers 3, 17. 

[3] 0. Tietjens: Beiträge zur Entstehung der 
Turbulenz. Z. angew. Math. Mech. Bd. 5 
(1925), S. 200— 217. 


Göttingen. H. Görtler. 452 


BUCHBESPRECHUNGEN 


Dozent Dr.-Ing. habil. ERNST ECKERT, Die Be- 
rechnung des Wärmeübergangs in 
der laminaren Grenzschicht um- 
strömter Körper. (VDI-Forschungsheft 416, 
Beilage zu „Forschung auf dem Gebiete des In- 
genieurwesens“, Ausgabe B, Bd. 13, Sept./Okt. 
1942.) 24 S. m. 29 Bildern, 3 Arbeitsblättern und 
12 Zahlentaf. Berlin 1942, VDI-Verlag G.m.b.H. 
Preis brosch. 5 M. 

Da bei Wärmetauschern vielfach laminare Grenz- 
schichten auftreten, stellt ihre theoretische Be- 
handlung ein technisch wichtiges Problem dar. Für 
nieht zu große Strömungsgeschwindigkeiten und 
Temperaturunterschiede kann man dabei Dichte, 
Zähigkeit und Temperaturleitfähigkeit als kon- 
stante Stoffwerte ansehen und die Wärmeerzeu- 
gung durch. innere Reibung vernachlässigen. Die 
laminare Strömungsgrenzschicht wird dann von 
den Prandtlschen Grenzschichtgleichungen (Bewe- 
gungs- und Kontinuitätsgleichung) beschrieben, 
während man die zugehörige Temperaturgrenz- 
schicht mit Hilfe der Energiegleichung berechnen 
kann. Die Grenzschichtgleichungen lassen sich 
für den Fall keilförmiger Körper, bei denen die 
Geschwindigkeit in der äußeren Potentialströmung 
nach einem Potenzgesetz von der Bogenlänge 
längs der Wand abhängt, streng lösen. In geeig- 
neter dimensionsloser Darstellung erhält man zu 
jedem Keilkörper ein bestimmtes Grenzschicht- 
geschwindigkeitsprofil. Diese bereits von Falkner 
und Skan und Hartree berechneten Profile werden 
nochmals dargestellt. Entsprechend wird vom 
Verf. jetzt auch die Energiegleichung behandelt. 
Die zu den jeweiligen Geschwindigkeitsprofilen ge- 
hörigen Temperaturprofile werden für den Fall 
konstant vorgegebener Wandtemperatur berechnet. 
Aus dem Temperaturgradienten an der Wand er- 
gibt sich die jeweilige Wärmeübergangszahl. Die 
so gewonnenen exakten Lösungen werden dann 
als Grundlage eines Näherungsverfahrens genom- 
men, mit dessen Hilfe sich die Geschwindigkeits- 
und Temperaturgrenzschicht an beliebigen Kör- 
performen berechnen läßt. Ausgehend von der be- 
kannten Lösung im Staupunkt wird die Änderung 
der Grenzschicht unter der Annahme berechnet, 
daß das Anwachsen der Grenzschichtdicke an 
einer Stelle des gegebenen Körpers ebenso groß 
ist wie bei einem Keilkörper an der Stelle, wo 
die Grenzschichtdicke und das Geschwindigkeits- 
gefälle mit den entsprechenden Werten an dem 
aan Körper übereinstimmen. Zur bequemen 
ösung der dabei auftretenden Differentialglei- 
chung sind einige Rechenblätter beigegeben. Ge- 
genüber den bisher bekannten Rechnungsverfah- 
ren kommt Verf. mit einem sehr viel kleineren 
Zeitaufwand aus. Zum Schluß wird das Verfahren 
an einer Reihe von Beispielen erläutert und mit 


Messungen und. den bisher bekannten strengen 
J,ösungen und anderen Näherungsverfahren ver- 
glichen, wobei sich eine befriedigende Überein- 
stimmung ergibt. 
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Die Tätigkeit der Baltischen Geo- 
dätischen Rommission iin den Jahren 
1938—1941. Berichtet von dem Präsidium. 113 S. 
Helsinki .1942, Osakeyhtio Weilin & Göös, Aktie- 
bolag. 


Diese Veröffentlichung enthält außer einem all- 
gemeinen Tätigkeitsbericht einen Landesbericht für 
Dänemark von N. E. Nörlund, ginen Landes- 
bericht für Finnland von J.Bonsdorff wad einen 
Landesbericht für Schweden von G. A. Rune. 
In Dänemark wurden ein hydrostatisches Nivelle- 
ment mit Hilfe eines mit Wasser gefüllten Blei- 
rohrs, eine Breitenbestimmung der trigonometri- 
schen Station Kopenhagen, Schwerebestimmungen 
mit einem dynamischen Schweremesser auf sechs ' 
Stationen und solche mit zwei verschiedenen sta- 
tischen Schweremessern auf zahlreichen Stationen, 
vier Basismessungen sowie ein fast 10000 km 
langes Feinnivellement ausgeführt. In Finnland 
wurden die Triangulation, die astronomischen. 
Arbeiten, das Feinnivellement und die Ausglei- 
chungsrechnungen fortgesetzt; ferner wurden ver- 
schiedene, auf die Invardrähte und die Basismes- 
sungen sich beziehende Untersuchungen ausgeführt. 
In Schweden wurden Triangulationen, Grundlinien- 
messungen, astronomische Messunsen sowie 
Schwerkraftsbestimmungen durchgeführt. 

Außer diesen Tätigkeitsberichten enthält die Ver- 
öffentlichung noch Mitteilungen von G. Ner- 
gaard über „Un gravimötre nouveau et des me- 
sures & l’isle de Bornholm“, B. Wideland über 
relative Schwerebestimmungen in Schweden im 
Jahre 1941, K. Ros&n über zwei Sehnendreiecks- 
formeln und den Satz von Legendre, V. R. Ölan- 
der über die Ausgleichung von unsymmetrisch 
angeordneten Richtungsbeobachtungen und 1. 
Bonsdorff über zufällige und systematische 
Fehler bei den Winkelmessungen der Triangulation 
1. Ordnung des Finnischen Geodätischen Instituts. 
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Dr. K. SCHWIDEFSKY, wissenschaftl. Mitarbei- 
ter d. Opt. Werke Zeiß, Jena, Einführung in 
die Luft- und Erdbildmessung. 3. er- 
weiterte u. verbess. Aufl. VI+ 176 S. m. 85 Abb,, 
3 schwarzen und 2 farbigen Tafeln im Text, 
1 schwarzen Tafel, 1 farbigen Brille und 2 Stereo- 
bildern im Anhang. Leipzig und Berlin 1942, Ver- 
lag B. G. Teubner. Preis geb. 7,80 M. 
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Zuschriften 
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ZUSCHRIFTEN AN DEN HERAUSGEBER 


Zur Näherungskonstruktion 
der logarithmischen Spirale. 


Die in dieser Zeitschrift Bd. 22 (1942) S. 168 bis 
169 veröffentlichte Kleine Mitteilung des Herrn 
H. Meincke, „Annäherung der logarithmischen 
Spirale durch Kreisbogen“, ist eine nur gering- 
fügig veränderte und erweiterte Wiedergabe einer 
unter Herrw Meinckes und meinem Namen ver- 
öffentlichten Arbeit „Eine Näherungskonstruktion 
für die logarithmische Spirale‘‘ (Deutsche Math. 
Bd. 3 (1938) S. 269 bis 272). Da ich für die Rich- 
tigkeit des Inhalts und der Formulierung dieser 
Arbeit die alleinige Verantwortung trage, da ferner 
Herr Meincke sie ohne Nennung meines Namens 
zitiert und schließlich in die Mitteilung Unklar- 
heiten und Fehler hineinbringt, die in meiner Dar- 
stellung nicht vorhanden waren, fühle ich mich 
zu folgender Richtigstellung veranlaßt: 


Es handelt sich um eine Näherungsk onstruk- 
tion der logarithmischen Spirale, für die Be- 
rechnung der Koordinaten einzelner Punkte 
bringen die Näherungsformeln keinen praktischen 
Vorteil. Auch bei Berechnung von Bogenlänge 
und Flächeninhalt wird man kaum Näherungs- 
formeln benutzen, da die durch Integration ge- 
wonnenen exakten Formeln 

r? — a2 


r—a e 
gs ——, F= —— 
co8T 4ctgrt 


hinreichend einfach sind. 


Falsch ist Formel (4) der Mitteilung; sie gilt bei 
kleinem 4 (s. Formel (8)) näherungsweise. Daher ist 
auch die Bemerkung im vorletzten Absatz unrich- 
tig, daß die in der Deutschen Math. veröffentlichte 
Zahlentafel entbehrlich sei. Bei gegebenem r und « 
bestimmt sich ® aus einer transzendenten Glei- 
chung. Die Tafel soll dem Konstrukteur die Mühe 
der Auflösung ersparen. Sie zeigt ferner, daß man 
sich in allen praktischen Fällen auf die Werte 
a=15°, 30°, 45° und 60° beschränken kann. 
Ferner ist = r+ = (vor Formel (5)) unrichtig; 


- 


a . 
es muß » = r—Z——ı nd oder zumindest on 


—5— tr heißen. In (10) bis (15) der Mitteilung 


muß das Argument der Sinusfunktionen überall 


a+t w statt r— > und o statt + = heißen. Wenn 


diese Winkel auch bei kleinem 4 näherungsweise 
Supplementwinkel sind, so darf hier das 4 nicht 
vernachlässigt werden, da der Gegenstand des Be- 
weises gerade die Güte der Annäherung ist. Die 
ohne Beweis mitgeteilte Grenzwertformel (6) ist 
nach Setzen von Klammern zwar richtig, ist aber 
nach dem Gesagten für Konstruktion und Rech- 
nung ohne Bedeutung. 


Berlin. Günther Schulz. 


Diese Zuschrift wurde Herrn Meincke zugesandt. 
W. 470 


Bemerkung zu H. Görtler, Berechnung 
von Aufgaben der freien Turbulenz auf Grund 
eines neuen Näherungsansatzes. Z. angew. Math, 
Mech. Bd. 22 (1942) S. 244 bis 254. 


Ich weise auf zwei Druckfehler hin, die leider 
nicht ausgemerzt wurden: 1. Seite 246, Zeile 19, 
ließ F’(0) =1 statt F’ (0) = 0; 2. Seite 250, erste 
Gleichung (2.7), rechte Seite, ließ y x/s statt Ys/x. 


Bei dieser Gelegenheit möge folgende ergän- 
zende Bemerkung zum ersten der durchgerechneten 
Beispiele (Vermischungszone zweier Strahlen) an- 
gefügt werden. Die Rechnung lieferte in den aus 
geführten ersten vier Schritten der Iteration außer- 
ordentlich rasch klein werdende Korrekturen der 
Ausgangsnäherung im praktisch interessierenden 
&-Intervall (etwa |Z!<2) und führte damit zu der 
erstrebten guten Näherungslösung für die Ver- 
mischungszone. Zu der in der Arbeit nicht näher 
untersuchten Frage jedoch, ob sich die exakte 
Lösung in Gestalt der dort nur formal angeschrie- 
benen unendlichen Reihe (1.9) für alle & 
in x SESoo wird darstellen lassen, ist fol- 
gende Warnung angebracht. Nach orientierenden 
Rechnungen, die ‚Herr W. Mangler kürzlich 
durchführte, ist Konvergenz des Verfahrens in 
ö|= oo nicht wahrscheinlich. Das leuchtet auch 
ein, denn das Iterationsverfahren wird wohl bis 
auf einen gemeinsamen für !E|>  asymptotisch 
exponentiell verschwindenden Faktor aller Reihen- 
glieder für sehr große |ä| im wesentlichen auf 
eine Entwicklung nach steigenden Potenzen von & 
führen, die dann natürlich in || — oo (wesentlich 
singuläre Stelle der- Lösung) nicht konvergieren 
kann. Im übrigen kommt dieser Frage im Hinblick 
auf die schon am Rande des praktisch interessie- 
renden #-Intervalls nicht mehr zutreffende physi- 
kalische Annahme der Theorie nur ein unter- 
geordnetes Interesse zu. 
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Nachtrag zu der Kleinen Mitteilung: „Eine 
Schieberanordnung für die Schlüssel- 
gleichung 


NE +YA) +R WM) = Ta hun), 


wobeinurdasEndresultatabzulesen 
ist.“ Z. angew. Math. Mech. Bd. 22 (1942) S. 169 
bis 170. 


In Ergänzung der obigen Mitteilung wird emp- 
fohlen, den Rechenstab mit der fı-Kurve (Schie- 
ber 1) durch einen Kreisrechenschieber zu ersetzen, 
dessen äußerer Ring die fı-Kurve als Rand hat. 
(Polarkoordinatendarstellung.) Dieser Kreisschie- 
ber wird je nach Verwendungszweck auf dem 
festen Teil oder auch auf der Zunge des Schie- 
bers 2 befestigt. Gerechnet wird mit dieser An- 
ordnung, wie s. Zt. für die zwei Stäbe beschrieben. 


Durch Einführung des Kreisrechenschiebers ist 
die Anordnung wesentlich leichter und einfacher 
zu handhaben. 


Böblingen. H. Schminke. 4% 
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